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INTRODUCTION 

Actuarial science is peculiarly dependent upon the Theory of 
Probabilities, the solution of many of its problems is best effected 
by resort to the Differential and Integral Calculus and in practical 
work the Calculus of Finite Differences is almost indispensable. 
Excellent text-books on these subjects are, of course, available but 
none of them has been written with the special requirements of 
the actuary in view. In beginning his training the student is, 
therefore, confronted by the difficulty of judicious selection and in 
the circumstances it has appeared to the Council of the Institute 
of Actuaries that a mathematical text-book sufficiently compre- 
hensive, with the standard works on Higher Algebra, to provide the 
ground- work of an actuarial education would be of great value. At 
the request of the Council, Mr Alfred Henry has undertaken the 
preparation of such a work and the resulting volume is issued in 
the confident expectation that it will materially lighten the toil of 
those who essay to qualify themselves for an actuarial career. 

A. W. W. 

May 1922. 
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AUTHOR'S PREFACE 

Actuarial science is essentially practical in that, whilst it is based 
on the processes of pure mathematics, the object of the worker 
must be to produce a numerical result. 

For this reason it is necessary for considerable prominence to be 
given, in the curriculum of the actuarial student, to the subject of 
Finite Differences, and it thus becomes convenient, in the study of 
those subjects not included under the heading of Algebra, to deal 
with this part of the syllabus first and, notwithstanding certain 
theoretical objections, to treat the fundamental propositions of the 
Differential and the Integral Calculus as being, substantially, 
special cases of similar propositions in Finite Differences. The 
subjects enumerated cover so wide a field that it has been necessary 
to exercise considerable compression and to include only such 
problems as are requisite for a proper knowledge of the subjects 
within the syllabus. 

In the chapter on Probability it will be seen that the numerical 
or "frequency" theory of probability has been adopted. Having 
regard to the practical nature of the actuary's work, it is thought 
that strict adherence to this aspect of the subject is necessary if 
the student is to acquire sound views from the outset. The subject 
of Inverse Probability has been excluded from the examination 
syllabus in recent years and for this reason it is not introduced into 
the present work. 

In conclusion the author would wish to tender his best thanks 
to many colleagues and other members of the Institute of Actuaries 
for their kind assistance and useful criticisms. In this connection 
he is particularly indebted to Mr G. J. Lids tone, who was good 
enough to read the chapters relating to Finite Differences and 
made many valuable suggestions. 

A. EL 
August 1922. 
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NOTE AS TO SYMBOLS 

In a few cases certain mathematical symbols have been employed, 
the use of which in modern mathematical work may be unfamiliar 
to students. 

! Used instead of | to represent a. factorial. Thus 

nl=\n==n(n — l)(n-2) 1, 

where n is positive or integral. 

== Denotes an identity which is true for all values of the variable. 

[ ) Used instead of n C rt (n, r) or — r? r , . 

-*. Tends to the value ; has the limiting value, e.g. 
As#-»-oo, (l + -)*-»- e n and e~ x -^0. 
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CHAPTER I 

FUNCTIONS. DEFINITION OF CERTAIN TERMS. 
GRAPHICAL REPRESENTATION 

1. When the value of a certain quantity y depends upon, or 
bears a fixed relation to that of another quantity, x, y is said to be 
a function of x, and the relationship is written as y =f(x). 

[Other notations used are u X) v x , <f> (x), etc.] 

Thus, we may have y = x 2 , y=a x , y = sinx, y=* log x. The 
quantities forming the right-hand sides of the equations are all 
functions of x. 

When expressed in this way the relationship of y to x is said to 
be explicit. But if, for example, ax 2 + 2bxy + cy 2 = 0, it is clear that, 
whilst the value of y depends upon that oix, it cannot be determined 
in any case by direct substitution of the value of x. The relation- 
ship in such circumstances is said to be implicit In this particular 
example the relationship can, of course, be made explicit by solving 
the equation for y in terms of x. 

The quantities x and y are called variables. 

The quantity x is called the independent variable since it may 
assume any value, whereas y is called the dependent variable since 
its value depends upon that of x. 

The independent variable is sometimes referred to as the argu- 
ment. 

2. Functions may, of course, involve more than one variable. 

For example we may have y = aPz* or u = — j an( J a ll these 

z 

variables may be involved implicitly. A function of more than one 

variable would be denoted by f(x,y,z), u x . y , etc. 

If we have a function u, such that 

u = aPy^z? + aP'yVz ' + ... 

and a + 6 + c = a' + 6' + c / =...=n, 

then u is said to be a homogeneous function ofn dimensions. 

H. T.B.I. 1 
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••'•*«• QecUtitiefor Co-ordinates. 

3. A function of two variables may be represented graphically 
according to the scheme of the figure shown. 

Let two straight lines of indefinite length OX, OF be drawn at 
right angles to each other. The point is 
called the origin, OX is called the axis of x 
and Y the axis of y. 

Then if ON is measured along the axis of 
x, equal in value to x, and at that point a 
perpendicular line PN is drawn equal in 
value to y =f(x), ON is called the abscissa 
and PN the ordinate of the point P. The 
convention is taken that measurements of x in the direction OX are 
considered to be positive and those in the contrary direction nega- 
tive. Similarly measurements of y in the direction OY are treated 
as positive and those in the contrary direction as negative. The 
point P, written for convenience as (x, y), is thus completely deter- 
mined from given values of x and y. 

If for every value of x the corresponding value of y were plotted 
on a diagram such as the above a continuous curve would be 
obtained which would be the graphical representation of the equa- 
tion y =f(x). It is, of course, impossible in practice to plot every 
value of the function, but generally a few values can be filled in so 
as to enable the curve to be drawn by sight. 

It should be noted that the lines OX and OY, and consequently 
the origin 0, can be chosen quite arbitrarily and that the position 
of the point P can be fixed, in the manner indicated, with reference 
to any suitable axes of co-ordinates. 

4. The point P can be fixed with reference to its distance from 
two straight lines not at right angles to each other, the distance 
PN— y being measured along a line parallel to the axis of y. 

The values x and y corresponding to a given point P are called 
the rectilinear co-ordinates (or simply the co-ordinates) of the point 
P. Where the axes are at right angles to each other the system 
can be distinguished, if necessary, by referring to x and y as the 
rectangular co-ordinates of P. 

As a rule rectangular co-ordinates are the more convenient to use 
in practice and lead to simpler results. Unless otherwise expressed 
it is to be understood that rectangular co-ordinates are implied. 
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5. The following examples give simple cases of the graphical 
representation of explicit functions. t 

(i) The equation x = a clearly represents a straight line parallel 
to the axis of y and at distance a from it; for the value of x at any 
point is constant and equal to a. 

(ii) The equation y = mx represents a straight line passing 
through the origin and making an angle 6 with the axis of x } where 
tan 8 = m; since at any point the ratio y : x is constant and equal 
to tan 9. 

(iii) 
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Let AB be any straight line cutting the axes of x and y respec- 
tively at the points A and B, so that OA = a and OB = b. 

Let P be any point on the line AB, of which the co-ordinates are 
(#, y). Then, if perpendiculars PN and PM be dropped upon the 
axes of x and y, MP = x and NP = y. 



Also 



and 



x = ON = PB 

a OA AB 

y_NP_AP 

6" OB" AB' 

x y PB + AP 
a + b~~ AB 

= 1. 



#2/ 



Hence - + t = 1 is the equation of the straight line AB. 
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4 FUNCTIONS 

6. An implicit function can be similarly represented. For ex- 
ample, it is obvious from the ordinary properties of the circle that 
the implicit relationship x 2 + y* = a 2 represents a circle of radius a 
with its centre at the origin. 

Note. The function y = a + bx + ex 2 + dx* + . . . is sometimes 
called a parabolic function, since the equation y = a + bx + ex* 
is represented graphically by a curve which is known as a 
parabola. 

7. It does not follow that for every value of x there will always 
be a real value of y. 

Thus, consider the function y 2 = (x — a) (x — b) (x — c), where 
c > b > a. If a? is negative, the right-hand side of the equation is 
negative and y can have no real value. If x is positive and < a, the 
position is the same. If, however, x > a and < 6, then the right- 
hand side is positive and y has a real value; but when x > b and 
< c, y is again unreal and remains so until x > c when a real value 
of y results for each value of x. 

The form of the curve is shown below, where OA = a, OB = b 
and OC — c. 

In circumstances such as these, 
where one or more parts of a 
curve are isolated from the 
others, the function and the 
curve representing it are said to 
be discontirmous. 

8. It is convenient here to 
introduce the conception of the 

limiting value of a function, or simply a limit 

If y =f(x) and y continuously tends towards a certain value and 
can be made to differ by as little as we please from that value, by 
assigning a suitable value to x, say a, then f(a) is said to be the 
limiting value off(x) when x tends to the value a. 

A convenient notation is as follows : 

y~*~f( a ) when x-*~a. 
Also /(a) would be expressed as Lt f(x). 
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GRAPHICAL REPRESENTATION 5 

Thus let y = . By writing y in the form 1 — we see that 

x x 

by making x indefinitely great, we can make the value of y differ 

from unity by as little as we please. 



Thus 



1*2=5-1. 



9. We will now give an example of another form of discontinuity, 



\x — al 



shown 



and for this purpose we will take the curve y- 
below. 

Here it will be seen that as x -+■ a, the value of y -*. oo . Similarly 
if x -*- oo , y -*■ 1. 

Thus if we draw two lines, one y 

PN parallel to the axis of y and at 

distance a from it, and the other tL 

QM parallel to the axis of x and at ^X 

unit distance from it, the curve will N 

continuously approach these lines 

but will not actually touch them 

except at an infinite distance from the origin. 

Such lines are called asymptotes to the curve. 

In general, actuarial functions are finite and continuous ; but in 
mathematical work, as will be seen later, attention to these points 
is necessary in the consideration of certain problems. 
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10. Explicit functions involving two variables, and implicit 
functions of three variables, can be expressed in a diagram of three 
dimensions. Thus if we have z=f(x,y) or f(x, y> z) — we may 
measure x, y and z by reference to their perpendicular distance, not 
from two lines, but from three planes at right angles to each 
other. 

A useful simile is that of the floor and two adjacent walls of 
a rectangular room of indefinite size. Any point in space is fixed 
by reference to its perpendicular distance from the two walls and 
the floor, x. and y being the respective distances from the walls, 
and z that from the floor. The conventions as to positive and 
negative values are similar to those previously explained. 
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11. As an example, if z = xPy 2 , the values of the function for unit 
intervals in the values of x and y are shown in the following table : 



Value 
of y 


Value of x 





1 


2 


3 


4 





1 
2 
3 

4 










1 

4 

9 

16 




4 

16 

36 

64 




9 

36 

81 

144 




16 

64 

144 

256 


... 



We will conceive the floor as being covered with a linoleum of 
chess-board pattern, the sides of the squares being at unit distance 
apart. The perpendicular distance from the floor of all points in 
these squares is nil and, therefore, z has the value 0. The corners 
of these squares will thus represent the various points (0, 0, 0), 
(0, 1, 0), (1, 1, 0), (1, 0, 0), (2, 1, 0), etc. If we were to erect at each 
corner a peg of height equal to the appropriate figure taken from 
the above table, the tops of the pegs would give points on the surface 
representing the equation z — a?y 2 . If the value of z were plotted 
for every possible combination of values of x and y> we should, 
of course, obtain the continuous surface corresponding to z = x*y\ 

Polar Co-ordinates. 
12. An alternative method of defining a point in a plane surface 
is as follows. 

Let an origin be taken and a fixed line OX be drawn from it ; 
then the position of any point P is known if the distance OP and 
the angle XOP are given. 

<-P 




Thus if OP = r and Z XOP = 0, r and are known as the polar 
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co-ordinates of P, and the point P can be written as (r, 0). The 
distance OP is called the radius vector and the angle XOP is called 
the vectorial angle. 

The convention adopted is that the angle XOP is reckoned 
positive if measured from OX in a direction contrary to that in 
which the hands of a clock revolve, and negative if measured in the 
reverse direction. 

Further, the radius vector is considered positive if measured 
from along a line bounding the vectorial angle, and negative if 
measured in the opposite direction. To illustrate this system, let 
PO be produced to a point Q such that 0Q = OP = r. Then the 
point Q may be written alternatively as (r, it + 0) or (— r, 0). 

13. The relation between rectangular and polar co-ordinates can 
be easily established. For if OX be taken as the axis of x then 
OY, the axis of y, is perpendicular to it. Also let PNbe drawn 
from the point P perpendicular to the axis of x. 




Then, clearly, if x, y be the rectangular co-ordinates of P, 
x = ON = OP cos = r cos 0, 
y=PN=OPsm0 = r sin 0. 

Any equation in rectangular co-ordinates can therefore be trans- 
formed into an equation in polar co-ordinates by the above substi- 
tutions. 

14. Three simple examples of the graphical representation of an 
equation in polar co-ordinates are now given. 

(i) The polar equation r = a clearly represents a circle of radius 
a with its centre at the origin; since the radius vector is constant 
and equal to a. 
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(ii) In the diagram shown in § 13, if ON = a and PN be produced 
indefinitely in either direction, then the polar equation of the 
straight line so obtained will be 

r cos 6 = a, 

since if P be any point in the line 

OP cos = r cos - ON = a. 

(iii) Let OA be a diameter of a circle OP A of radius a. Then 
if P be any point on the circle such that OP = r and £AOP = d, 
OP = 0A cos 6 = 2acosd. 




The polar equation of the circle, if be the origin, is therefore 

r = 2a cos 0. 



CHAPTER II 

FINITE DIFFERENCES. DEFINITIONS 

1. The subject or calculus of Finite Differences deals with the 
changes in the values of a function (the dependent variable) arising 
from finite changes in the value of the independent variable (see 
Chapter I, § 1). 

Many questions arise which can be dealt with on systematic lines, 
but probably the most important problems which require to be 
solved in actual practice, and' with which we are concerned at this 
stage of the subject, are the summation of series, and the insertion 
of missing terms in a series of which only certain terms are given. 

It will be convenient to proceed in the first place to some ele- 
mentary conceptions and definitions. 

2. If we have a series consisting of a number of values of a 
function, corresponding to equidistant values of the independent 
variable, and from each term of the series we subtract the algebraic 
value of the immediately preceding term, we shall obtain a further 
series of equidistant terms. The process is known as differencing 
the terms of the series, and the terms of the new series are known 
as the first differences of the original terms. By repeating the pro- 
cess with the terms forming the first differences, we shall obtain 
a further series forming the second differences of the original function, 
and so on. Thus if we have f(x) for the first term of the series 
and/(# + A) for the second term, the first difference of f(x) is 
f(x + h)-f(x) and is designated A/(#). The second difference of 
f(x) is A/ (x + h) - A/ (x) and is designated A 2 /(#). This may be 
set out as in the following scheme : 



Function 


First Differences 


Second Differences 


Third Differences 


/(*+A) 

/(*+3AJ 
/(tf+4*J 


/(*+2A)-/(*+A) 

'/(#+3A)-/(# + 2A) 
/(a?+4A)-/(*+3A) 


/(a?+2A)-2/(*+A) 
/(*+3A)-2/(#+2A) 

+/(*+2A) 


f(x+3h)-3f(x+2h) 
+ 3/<*+A)-/<*) 

/(*+4A)-3/(#+3£) 
+ 3/(*+2A)-/(#+A) 
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FINITE DIFFERENCES 



The first term of the series is known as the leading term and the 
terms in the top line of differences are known as the leading dif- 
ferences of the series. 

It must be clearly understood at the outset that A is merely a 
symbol representing the operation of differencing f(x) once ; it is 
in no sense a coefficient by which f(x) is multiplied. This point is 
dealt with again in § 5. 

3. An examination of the character of the series which ultimately 
results from the process of differencing repeatedly, leads to the 
development of certain important theorems. Before proceeding 
further, it will be helpful to give a practical example. 

Example 1. Obtain the differences of the series given by f(x) = a?, 
where x has all integral values from 1 to 6. 





/« 


Firpt 


Second 


Third 


Fourth 


X 


Differences 


Differences 


Differences 


Differences 


1 


1 


7 


12 


6 





2 


8 


' 19 


18 


6 





3 


27 


37 


24 


6 




4 


64 


61 


30 






5 


125 


91 








6 


216 











It will be observed in the above example that the fourth and, 
therefore, all higher differences are zero ; it will be seen later that 
this would equally have been the case had more terms of the series 
been taken. We can therefore construct all the remaining terms 
of the series by a process of continuous addition. 



, 4. Although most functions with which the actuary has to deal 
are not of the simple character of that shown above, yet it will 
usually be found that the differences of the function for which 
further values are required tend to the value zero and are sus- 
ceptible to treatment by methods which will be developed subse- 
quently. 

The student should obtain confirmation of this fact and insight 
into the character of certain series by taking out the differences of 
tabulated functions such as logarithms, annuity-values, etc. 
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5. Before proceeding to the consideration of the various problems 
which arise, it is necessary to develop certain fundamental formulas. 

In § 2 A has already been defined as the symbol of the operation 
by means of which the value of f(% + h)—f(x) is obtained. 
Similarly, it is customary to use the symbol E as representing 
the operation by which the value of f(x) is changed to the value 
/(# + h), so that 

ty<*> -/(* + h ) -/(«) + A/O). 
It must be carefully remembered that these symbols represent 
operations only and must be interpreted accordingly. Thus EPx 2 is 
clearly not the equivalent of (Ex) 2 ; the former expresses the result 
of operating twice upon the function x 2 in the manner indicated 
above, giving a value (x + 2h) 2 , whereas in the latter case the opera- 
tion is applied once to the function x and the resulting term (x + h) 
is squared. 

6. If, then, these symbols are found to obey the ordinary alge- 
braical laws, they can be dealt with algebraically provided always 
that the results are interpreted symbolically in relation to the 
function which is the subject of the operation. This principle is 
known as that of Separation of Symbols or Calculus of Operations. 

The algebraic laws referred to above comprise : 

(1) The Law of Distribution. 

(2) The Law of Indices. 

(3) The Law of Commutation. 
Taking these laws in succession : 

(1) The symbol A is distributive in its operation, for 

-[/i(«)+/.W+/.« + -] 

-[/i(*+*)-/i(*)]+[/^+*)-^wa 
+[/■(*+*)-/•(*)]+••• 

-A/ 1 (*) + A/ 1 (*) + A/,(«)+.~. 

Similarly the symbol E is distributive, for 

-^/i (*) + */■(*) + */•(■) + •••• 

(2) The symbol A obeys the law of indices, for in the case of 
positive integers the symbol A m /(#) represents the operation, 
repeated m times, of differencing f(x). 
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Thus A w /(#) = (A A A ... m times) /(#), 

.-. A n A m /(#) = (AAA ...n times) (AAA ... ra times)/(#) 
= (AAA ...(m + ri) times)/(#) 

= A n+m /0*0- 
Similarly it may be shown that the symbol E obeys the law of 

indices. 

(3) The symbol A is commutative in its operation as regards 

constants, for, if c be a constant, 

A[o/(«)]-o/(« + A)-o/(«) 
-c[/(« + A) -/(«)] 
-cA/(*). . 
The like result can be deduced as regards i?. 

7. It follows that, since 

Ef(x) = (l + A)f(x) f 
therefore 2? = 1+A 

and A = #-l. 

The two operators are thus connected by a simple relation, which 
will be found later to lead to important results. 

8. As an example of the manner in which the relationship 
between the operations represented by E and A can be utilised in 
the solution of problems, we may take the following : 

Example 2. Prove that 

/(0) + */(l) + g/(2) + g/(3) + ... 



:*[/(0) + *A/(0) + ^A«/(0)+...] 



Since 
/(l) = Ef(0) = (1 + A)/(0); /(2) = tf ■/(<>) = (1 + A)»/(0), etc., 

we have /(0) + */(l) + |/(2) + g/(3) + ... 



= [l+*(l + A) + g(l + A)» + ...]/(0) 



= [V (1+A) ]/(0) 



= e* l+«A + *^+...l/(0) 

/(0)+tfA/(0) + gA°/(0) + 



= 6* 



CHAPTER III 

FINITE DIFFERENCES. GENERAL FORMULAS 
AND SPECIAL CASES 

1. Starting from the relationship proved in Chapter II, § 7, it is 
now possible to develop two formulas of the utmost importance. 

2. To express f(x + mh) in terms of f(x) and its leading dif- 
ferences. 

By definition fix + mh) = E m f{x) 

= (l + A)-/(#), 
since by Chapter II, § 7, E = 1 + A. 

Expanding the above expression by the Binomial Theorem we 
have 

/( a; + mA) = [l + mA + (^)^ + (™)A'+... + (^)A»»]/(, c ) 

+ C)*'V(*) (D- 

Bearing in mind that the symbol A obeys the ordinary algebraic 
laws and that the Binomial Theorem holds for all values of the 
index, positive, fractional or negative, it will be realised that the 
above proof is perfectly general. 

It is instructive, however, to show how the formula may be 
deduced for fractional and for negative values of m. 

3. Fractional value. Let — be a positive proper fraction, and 
let the interval between the given terms of the series be h. It is 
required to find the value of / f x H — h). 

Also let f(x + h) -f(x) = A/(#), 

and /(* + £)-/(<> = $/(*). 



14 FINITE DIFFERENCES 

In the second case, the unit of differencing has been altered to 

- and, bearing in mind that n is a positive integer, we may write 
n 

at once from the theorem in the preceding article 

/(x + h) = (i + sy/(x), 

therefore (1 + A)f(x) = (1 + S) n f(x), 

and (l + A)»/(*) = (! + «)/(*)• 

Since m is also an integer, it follows that 

(i + sr/(*) -/(« + ™ *) - a + a)-/(«) 

-/(«)+5A/(«)+(S)Ay(«H— 

4. Negative value. It is desired to find the value o{f(x — mh). • 
Now, from the preceding theorems, it is clear that 

(l + A)~/(«-mfc) =/(«), 

therefore 

/(*-mA)-=(l + A)-/<«) 

= [l + (-m)A + (-;)AH...]/W 

-/<«) + (-m)A/(«) + (-*)A»/(«)+.... 

The above proofs show that the theorem holds universally and 
illustrate how the principle of Separation of Symbols can be applied 
whenever the symbols of operation obey the ordinary laws of algebra. 

5. To express A m f(x) in terms of f(x) and its successive values. 

= [^-m^ n - 1 +(2)^ w - 2 -...+(-l) w *]/(^) 

=f(x -r mh) - mf(x + m — Ih) + ( ™ jf(x + m — 2k) — . . . 

+ <-l)-/W (2). 

Alternatively both the above formulas can be easily proved by 
the ordinary methods of induction. 

Formula (1) also follows directly from the ordinary formula of 
Divided Differences (see Chapter VIII). This method has the 
advantage of showing directly the application of formula (1) to cases 
where m has a fractional or a negative value. 
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6. The above formulas are expressed in a form which applies in 
the most general way, i.e. when the interval of differencing is h and 
the leading term is/(#). It is clear, however, that by altering the 
unit of measurement the formula will be simplified although the 
result is not affected. Similarly by changing the leading term 
(which process corresponds to shifting the "origin") so that the 
leading term is expressed as /(0) a further simplification in form 
is made. 

If, therefore, the interval of differencing becomes unity and the 
leading term can be represented by /(0), the first formula can be 
written 

/(»)-/(0) + nA/(0) + g)A»/(0) + (3). 

An example will make this clear. 

Having given the values of/(10), /(15),/(20), etc. it is desired 
to express /(17) in terms of /(10) and its leading differences. 

The original formula (1) gives the value of /(10 + 1'4A), where 
h = 5, and therefore we write 

/(10 + V*h) =/(10) 4- 1-4 A/(10) + (1 ' 4 > ( ' 4) A*/(10) + ■ . ■ , 

where A, A 2 , ... are taken over the interval h. 

But the same result is secured if the unit of measurement is 
changed from 1 to 5 and if at the same time /(10) is made the 
initial term of the series, for then /(10), /(15), /(20), ... can be 
written as -F(O), F(l) t F(2), ... and the required value, viz./(l7), 
becomes jF(1/4) which by formula (3) is equal to 

J(0)+1'4AF(0) + (1 ' 4 ^' 4) A 2 ^(0) + .... 

7. The above formulas are of general application if sufficient 
terms of the series are known, but it is convenient at this stage to 
consider the particular forms taken by the differences of certain 
special functions. 

8. f{x)-a%\ 

A/(^) = a(^ + l) n -o^ n = a[^ n - 1 +(2)« n - 3 +...+l]. 

The result of differencing has been, therefore, to change the term 
involving the highest power of x from ax n to anx n ~ l (thus reducing 
its degree in x). 
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Similarly a further process of differencing will reduce the degree 
of x to n — 2 and the coefficient of the highest power of x will be 
an (n — 1). By repeating the process we arrive at the result that 
the nth difference of ax n is independent of # and is equal to a.nl. 
The (n + l)th difference is therefore zero. 

Corollary. It follows that the nth difference of 
ax n + foe 71 " 1 + cx n ~ 2 + ...+& 
is constant and equal to a.nl. 

9. f( x ) = x(x-l)(x-2)...(x-m+ 1). 

This expression is usually denoted by xl m \ 

&f(x) = (x+l)x(x-l)...(x-in + 2)--x(x-l)(x-2)...(x-m+l) 

= nix (x — 1) . . . (x — m + 2) 

= 7nx^ m " l \ 

Similarly A a /(#) = m (m - 1) ^ m "" 2 ). 

By repeating the process we arrive at the result 

A m /(a?) = m!, 

which is otherwise obvious from the preceding article since f(x) is 
of the wth degree in x. 



10. f(x) x ( £€ + l )( cc + 2)...(x + m-I)' 

Corresponding to the notation already used, this can be denoted 
by x^~ m) . 

A ^^ = (x+l)(x + 2)...(x + m) " x(x+l)(x + 2) ...(x + rn-1) 
— m 



~" x (x + 1) (x + 2) . . . (x + m) 


= — mx^' 


-m+i) 




Similarly 




A 2 /(#) = m (m + 1) ®h™*# 


Ld so on. 






11. 




Af(x) = a^-a x 
= a*(a-l). 


Whence 




A n f(x) = a* (a - l) n . 
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12. For many purposes it is convenient to have a table of the 
leading differences of the powers of the natural numbers. These 
can be represented as the differences of [a^a^o a ^d are sometimes 
known as the "Differences of 0." 

The following table gives a number of values of the first term 
and leading differences of [& m x n ] x=0> which, for convenience, can 
be written as A m n : 



n 


/(0) 


A 


A2 


A 3 


A* 


A6 


A 6 


1 


















2 







2 





— 


— 





3 







6 


6 





— 





4 







14 


36 


24 








5 







30 


150 


240 


120 






13. A working formula for constructing a table such as the 
above by a continuous process may be obtained as follows: 

& n f(x)=f(x + n)-nf(x + n-l) + (j^f(x + n-2)-.... 

[See formula (2).] 
Therefore, when f(x) = x m , 

A n x m = (x + n) m -n(x + n-l) m +(^(x + n~2) m ... > 
whence, putting x = 0, 

A n m = n m -n(n-l) m +(^(n-2) m -... 

= n [n^-in- l)(n - l)™" 1 + ("g 1 ) (»- 2)™" 1 - ...] 

= n [(1 + n-l)™- 1 - (n - 1) (1 +1^2)™-* 

^[A*- 1 ^- 1 ]^ (4). 

But /(l)=/(0) + A/(0), 

and A n " 1 /(1) = A M /(0) + A*/(0). 

Therefore 

[A*- 1 ***- 1 ]*.! = [A*- 1 ^- 1 ]^ + [A***"- 1 1-0 

_ ^n-iQm— l i A 71 © 771 " 1 . 
H. T.B.I. 3 
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Hence A n O m = n [A^O™- 1 + A n m " 1 ] (5). 

It follows that the differences of \aF%=* can be constructed from 
those of \x m ' r ] x=a0> and so on. 

To take an example from the table given above, 
A 4 5 = 4[A 8 4 + A 4 4 ] 
= 4 [36 + 24] = 240. 

14. By using the result given in § 9, it is possible to expand 
f(x) in terms of x^\ x^\ x&\ ... . 

Let f(x) = A + A x xU + A 2 ot* + A 3 x® + . . . . 

Then, putting x = 0, we* see that 

f(0) = A . 
Differencing both sides of the equation, we get 

Af(x) = A X + 2A 2 xW + SAtofl + . . . . 
Again putting x = 0, we find 

A/(0) = ^. 
By repeating the above processes, we obtain successively 

A'/(0) = 2U 2 , A>/(0) = 3!^„... A*f(0) = n\A n , 
whence 

A„ m . VWl A- A -^),... A.-&P. 
and 

/(«)-/(0) + «A/(0)4-^4V , (0)+ ... + £a»/(0)+ ... (6). 



CHAPTER IV 

FINITE DIFFERENCES. INTERPOLATION 

1. The subject of interpolation is one of the most important in 
Finite Differences and may be enunciated as follows. 

It frequently happens that we have given a number of values of 
f(x) corresponding to different values of #, and we wish to find 
a value of the function for some other value of x. If the form 
of the function is known or can be deduced from the given values, 
the problem is, of course, simple, although in many cases it is more 
convenient to proceed by the methods of Finite Differences. But it 
is frequently the case, especially in actuarial work, that the function 
cannot be expressed, algebraically or otherwise, in any simple form, 
and resort must be had to other devices. 

2. Looked at from the point of view of a problem in graphs, we 
may regard thegiven values of the function as representinga number 
of isolated points on a curve, and it is desired to plot a further point 
corresponding to a given value of the abscissa. 

It follows that if the form of the function (i.e. the equation of 
the curve) is unknown, some assumption must be made as to the 
relationship between the different values. The formulas of finite 
differences assume that this relationship can be expressed in the form 

y = a + bx + cx 2 + dx z + ... +kx n ~\ 
This assumes (see Chapter III, § 8) that all orders of differences 
higher than the (n — l)th vanish, but, as pointed out in Chapter II, 
§ 4, this assumption can be made without introducing important 
errors in practically all cases where actuarial functions are involved. 

3. The above equation contains n constants, and therefore n 
values of the function must be known if the values of the constants 
are to be determined. Conversely, if n values only are known and 
the methods of finite differences are to be applied, it must be 
assumed implicitly that all orders of differences higher than the 
(n — l)th vanish. 

4. The most obvious method of procedure is to obtain the n 
equations given by the n values of the function and to find the 
values of the constants therefrom. The assumed form of the function 

2-2 
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is then completely determined and the value corresponding to any 
value of x can be obtained. 

In the majority of cases, however, this is not the most simple 
method of working, for other devices can be adopted which will 
materially shorten the arithmetical work. It is important to note, 
however, that alternative formulas, in which the same values of the 
function are used, lead to identical results. 

In some cases there is scope for the exercise of the ingenuity of 
the solver, but usually the problems fall into the main categories 
which are illustrated in the following examples. 

5. Example 1. When n equidistant values of a function are given 
and it is required to find the value of some intermediate term or terms. 

This can be done readily by the application of formula (1) of 
Chapter III, or by the simpler formula (3). From the given values 
the successive orders of differences are calculated, and the result is 
obtained by direct substitution. 

Thus, taking the numbers living by the H M table at ages 45, 
50, 55, 60 and 65, it is required to find the value for age 57. 

In conformance with formula (3) the given values can be denoted 
by/(0),/(l), ..., so that the required value is/(2*4). Then 

2*4 x 1-4 



/(2-4)=/(0) + 2'4A/(0) + ' 



A'/(0) 



+ !±xM)^ At ^ 
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The working is as follows: 



X 


/<•) 


Af(x) 


A 3 /(s) 


A»/W 


AVW 





77918 


-5123 


-1106 


-389 


+75 


1 


72795 


-6229 


-1495 


-314 




2 


66566 


-7724 


-1809 






3 


58842 


-9533 








4 


49309 











From above: 

/(2-4) = 77918 + 2-4 (- 5123) + 1-68 (- 1106) 

+ -224 (- 389) - -0336 (+ 75) 
= 77918 - 12295-2 - 18581 - 871 - 2-5 
- 636751. 
The value given by the table is 63677. 
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The difference between the interpolated value and the true value 
is due to the fact that the interpolation curve, which is based on 
the assumption that all differences of higher order than the fourth 
vanish, represents only approximately the true function. 

6. Example 2. When the values given and the value sought 
constitute a series of equidistant terms. 

If there are n terms given of which n — 1 are known, then, as 
explained in § 3, it must be assumed that the {n — l)th order of 
differences is zero. 

Thus, using formula (2) of Chapter III, we have 

A- 1 /(0) = 0=/(«-l)-(n-l)/(r l -2)+(^ 1 )/(n-3)-... 

+<-i)"- , /(0). 

In this equation there is only one unknown quantity and its 
value can, therefore, be readily obtained. 
For example, if 

/(0) = log 3*50 = '54407, 
/(I) -log 3-51- -54531, 
/(2) = log 3*52 = -54654, 
/(4) = log 3'54 = -54900, 
and it is required to find log 3*53, i.e./(3). 
From above: 

A«/(0) - =/(4) - 4/(3) + 6/(2) - 4/(1) +/(0), 

whence /(3 )JJMMM) 

= -54777, 
which agrees with the true value to five decimal places. 

7. Example 3. If more than one term is missing from the com- 
plete series, a somewhat similar process may be followed. Thus, if 
two terms are missing, only (n— 2) terms are known and the (n— 2)th 
order of differences must be assumed to vanish. It is then possible 
to construct two equations: 

A^/(0) -/(»- 2) - (» - 2)/(n - 8) 4- . . . 4- (- 1)^/(0) = 0, 
A«/(D-/(»-l)-(»-2)/(n-2) + ...+(-ir»/(l)-0. 

From these equations, the values of the two unknowns can be 

calculated. 

Similarly if a larger number of terms is missing, the method 

can be extended. 
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8. Example 4. If several equidistant values are given, together 
with one isolated term. 

For instance, if three values /(0), /(l) and /(2) are given, to- 
gether with a further value f(h). Having four values of the function 
it must be assumed that the fourth order of differences is zero and 
it remains to find the values of the other three leading differences. 
The first two leading differences are obtained at once by differencing 
the first three terms of the series, and the value of the third dif- 
ference is then given by the equation 

/(*)«/(<>) + AA/(0) + (J) A*/(0) + (J) A'/(0). 

For example, taking the numbers living by the H M table at ages 
45, 46, 47 and 50, it is required to find values for ages 48 and 49. 

A A» 

/(0)= 77918 -954 -32 

/(1) = 76964 -986 

/(2) = 75978 

/(5) = 72795 
/(5) =/(0) + 5A/(0) + 10A'/(0) + 10A°/(0), 
A . /(0 ) = /(S)-[/(0) + gA/(0) + 10A*/(0)] 

72795 - [77918 - 4770 - 320] 



10 



= -33. 



The table is then completed by addition. Thus: 



Age 

X 


/w 


A 


A 2 


A3 


45 


77918 


- 954 


-32 


-3-3 


46 


76964 


- 986 


-35-3 


-33 


47 


75978 


-1021-3 


-38*6 


-33 


48 


74956-7 


- 10599 


-41-9 




49 


73896*8 


-1101*8 






50 


72795*0 









The work is checked by the reproduction of the value for age 50. 

The tabular values for ages 48 and 49 are 74957 and 73896 
respectively. 
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Values precisely the same as those obtained above would have 
been given if the two missing terms had been inserted by the 
method described in Example 3. It is instructive to confirm this 
by actual calculation and to compare the two methods of procedure. 

9. Example 5. Subdivision of Intervals. 

This problem arises when a series of equidistant terms of a series 
is given (usually every fifth term or every tenth term) and it is 
desired to find by interpolation the values of all the intermediate 
terms. 

The simplest method of procedure is to calculate from the given 
values the differences corresponding to the individual terms of the 
series (the subdivided differences) and thence to construct the table 
by summation. The calculation is checked by the reproduction of 
the values of the original terms. 

Thus assume that the given terms are /(0), /(l), .../(5) and it 
is desired to complete the series /(0), /(£), /(f), etc. It is con- 
venient to adopt the notation 

/(l)-/(0) = A/(0), 
and /(*)-/(<>) = 8/(0). 

The problem then becomes to express 8/(0), &f(0), ... in terms 
ofA/(0),A«/(0) 

Writing /(l),/ (2), ... in terms of the subdivided differences, 

/(0) = /(0) 

/(l)=/(0)+ 55/(0) + 105 2 /(0) + 1053/(0)+ 554/(0) + S 6 /(0) 

/(2) = /(0) + 105/(0) + 458 2 /(0) + 1205 3 /(0) + 210S 4 /(0) + 2528 6 /(0) 
/(3)==/(0) + 155/(0) + 1055 2 /(0) + 4555 3 /(0) + 136554/(0) + 3003^/(0) 
/(4) = /(0) + 205/(0) + 1905 2 /(0) + 11405 3 /(0) + 48455 4 /(0) + 155045 6 /(0) 
/(5) = /(0) + 255/(0)+300S 2 /(0) + 23008 3 /(0) + 12650S 4 /(0) + 53130SV r (0) 

Differencing successively both sides of the equation, we have 

A/(0)- 53/(0)+ 105 2 /(0) + 105 3 /(0) + 53*/ (0)+ 3 6 /(0) 

A/(l) = 5S/(0)+ 358 2 /(0)+ 1103 3 /(0)+ 2055 4 /(0) + 2513 6 /(0) 

A/(2) = 53/(0)+ 605 2 /(0)+ 335S 3 /(0) + 115534/(0)+ 27515 6 /(0) 

A/(3)- 5S/(0)+ 853 2 /(0)+ 6858 3 / (0) + 348034/(0) + ! 250 13 6 /(0) 

A /(4) = 53/(0) + 1 10S 2 / (0) + 1 1608 3 / (0) + 78053 4 / (0) + 37626S 6 /(0) 

A 2 /(0) = 253 2 /(0)+ 1003 3 /(0)+ 20034/(0)+ 250S 6 /(0) 

A 2 /(l) = ?53 2 /(0)+ 2258 3 /(0)+ 9508 4 /(0)+ 25005^(0) 

A 2 /(2) = 25S 2 /(0) + 3503 3 /(0) + 232534/(0)+ 97505^/(0) 

A 2 /(3)= 25S 2 /(0)+ 4758 3 /(0) + 4325S4*/(0) + 251253 5 /(0) 

A 3 /(0) = 1253 3 /(0)+ 75054/(0)+ 22503 6 /(0) 
A 3 /(l) = 1255 3 /(0) + 137554/(0)+ 72505 5 /(0) 
A 3 /(2) = 1255 3 /(0) + 200054/(0) + 153755 6 /(0) 

A 4 /(0) = 6255 4 /(0) + 50005 6 /(0) 
A4/(l)= 6255 4 /(0)+81255 6 /(0) 

A 6 /(0) = 31255 6 /(0) 
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Whence the values of S/(0), S 3 /(0), ... 8 5 /(0) can readily be 
obtained. 

10. Alternatively the formulas for S, S 2 , . . . can easily be written 
down by using the method of Separation of Symbols. 

Thus (1 + hyf (x) = (1 + A)/ (x). 

Therefore (1 + S)/(*) = (1 + A)V(#), 
and S/(*) = [(! + *)* -!]/(#) 

= [-2 A - -08 A 8 + -048 A 8 . . .]/ (#). 
Hence #*/(#) = ['2 A - '08 A 2 + '048 A 3 . . .] 3 / (x) 

• = [/04A 2 --032A 3 + -0256A*. ..]/(*) 
and so on. 

For convenience the coefficients of A, A 3 , . . . occurring in the 
values of S, S 2 , . . . are given, for the intervals 5 and 10, in the fol- 
lowing tables. 

Subdivision into 5 intervals 









Coefficient of 




Value of 










A 


A* 


A 3 


A* 


A 8 


a 


+ •2 


-•08 


+ •048 


- -0336 


+ -025536 


a 2 




+ •04 


-•032 


+ -0256 


-•02112 


a 3 






+ •008 


-•0096 


+ -00960 


a* 








+ -0016 


- '00256 


a 6 










+ -00032 



Subdivision into 10 intervals 



Value of 


Coefficient of 

i 


A 


A 2 


A 3 


A* 


A» 


a 
a 2 
a 3 
a* 

a 6 


+ •1 


-•045 
+ •01 


+ -0S85 
-•009 
+ •001 


- -0206625 
+ -007725 

- -00135 
+ •0001 


+ •01611675 
- -0066975 
+ -0014625 
-•00018 
+ -00001 
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11. The following example gives an illustration of the method 
of working. 

Given the present values, at 3 per cent, interest, of an annuity of 
1 per annum for 20, 25, . .. 45 years, it is required to find the inter- 
vening values. 

We have 



X 


/w 


A 


A 2 


A3 


A 4 


A 5 


20 


14-8775 


2-5356 


- -3483 


•0478 


-•0065 


•0007 


25 


17-4131 


2-1873 


--3005 


•0413 


-•0058 




30 


19-6004 


1-8868 


- -2592 


•0355 






35 


21-4872 


1-6276 


- -2237 








40 


23*1148 


1*4039 










45 


24-5187 













Applying the factors given in the former of the above tables, we 
have 

s /( 20 ) = [' 2A - ' 08 A 2 + -048A 3 - -0336A 4 + -025536A 8 ]/(20) 

= -2 x 2-5356 4- '08 x -3483 + '048 x -0478 + '0336 x -0065 

+ -025536 x -0007 
= •5375146752. 



Similarly 

S 2 m - -015642784, 
S 3 = -000451520, 



#=- -000012192, 
»_ -000000224. 



The table is then constructed by addition from these leading 
differences as shown below. 

It is necessary to consider how many decimal places should be 
retained in the working. This depends, in the first place, on the 
degree of accuracy desired in the result. Thus if four decimal 
places are required in /(#), our result must be given to at least 
one place more. Further, the range of the formula is 25 terms, and 
in the final term the values of S, S 2 , . . . S 5 are multiplied respectively 

ty \l) > \%)> •••{&)> which are respectively equal to 25, 300, 
2300, 12650, 53130. In view of the magnitude of these coefficients, 
we shall need to retain at least six decimal places more in the value 
of S 5 than are required in the final value of /(#). Since the co- 
efficients of the lower orders of differences are smaller than 53130, 



26 



FINITE DIFFERENCES 



a correspondingly smaller number of decimal places can be retained. 
In practice, however, there is little to be gained by cutting down, 
unless only a rough result is required. 

The working of the first five terms in the example is shown 
below. It will be noticed that the accuracy of the work up to this 
stage is checked by the exact reproduction of the value of /(25). 
The interpolated values agree exactly with the true values to four 
places of decimals. 



X 


/(*) 


8 


5 2 


5 3 


3* 


5 5 


20 


14-8775 


•5375147 


-01564278 


•000451520 


-•000012192 


•000000224 


21 


15-415015 


•5218719 


1519126 


439328 


11968 


224 


22 


15-936887 


•5066806 


1475193 


427360 


11744 


224 


23 


16443568 


•4919287 


1432457 


415616 


11520 


224 


24 


16*935497 


•4776041 


1390895 


404096 


11296 


224 


25 


17-413101 


•4636951 


1350485 


392800 


11072 


224 



12. Example 6. Lagrange's Theorem. 

We have now to consider the construction of an interpolation 
formula which will apply when the given terms of the series are 
not equidistant. 

Let n values of the function be given, namely 

f(a), XV, Ac),. ..f(n). 
Then the function must be assumed to be a parabolic function of 
x of degree (n — 1). (See § 3.) Assume therefore that the function 
can be represented as 

f(x) = A (x — b) (x — c) . . . (x — n) 
+ B (x — a)(x — c) ...(x — ri) 
+ C(x — a) (x — b) . . . (x — n) 
+ etc., 
there being n terms in all, each composed of (n — 1) factors multi- 
plied by a constant, the values of the constants having yet to be 
determined. It is clear that the right-hand side of the equation 
is of degree (n — 1). 

To find the values of the constants we proceed as follows: 
Put x = a. 

Then f(a) = A (a -6)(a - c) ... (a - n). 



Therefore A = 
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(a — b)(a — c) . . . (a — n) ' 



Similarly J - fl-a)ft- < oj...(ft-») ' 

and so on. 

Substituting these values of A, B,... in the original equation 
we have 

fM - f(n\ («-6)(«-c)-(«-«) , f(h s («-o)(g-o)...(g-n) 
/W-/W( _6)( _e)...(a-n) +/w (6-a)(6-c)...(6-n) 

+ ... + / W »- a )<*-g< <> -^" ...a> 

^ v / (n — a)(n--o)(n-c)... v ' 

By an obvious transformation, the formula can be put in a some- 
what simpler form for calculations, namely 

/(«) /(*) 

(% — a)(x-b)...(x-ri) (x — a)(a-b)(a — c)...(a—n) 

+ (b-aXx-b)...{b-n) + '" + (n-a)(n-6)(n-c)...(«-n) 

(2). 

In memorising the formula it should be noted that the denomi- 
nators are made up of the product of the algebraic differences of 
the values of the variable, the term {a — a) being replaced by (a — a) 
and so on. 

13. The formula is somewhat laborious to apply, and careful 
attention to signs is required, but it is convenient to use where only 
one or two unknown values of the function are required. Since the 
assumptions underlying it are precisely similar to those previously 
explained, its use in any particular case will give identical results 
with those which can be obtained by the use of the ordinary 
methods of Finite Differences where a sufficiently high order of 
differences has been taken into account. 

To illustrate this point and to provide an example of the use 
of the formula, we will calculate by Lagrange's formula the value 
for age 49 in Example 4. In this case we have 

/(0) = 77918, /(1) = 76964, /(2) = 75978, /(5) = 72795, 

and it is required to find the value of/(4). 
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We have accordingly 

m yxo) 

(4_0)(4-l)(4-2)(4-5)""(4-0)(0-l)(0-2)(0-5) 

/(l) /(2) 

+ (l-0)(4-l)(l-2)(l-5) + (2-0)(2-l)(4-2)(2-5) 

^(5-0)(5-l)(5-2)(4-5)' 
or, - £/(4) = - A/(0) + A/(l) - A/(2) - *V/(5). 

Whence we find/(4) = 73896'8, as before. 

It should be noted, as a check on the formula, that the sum of 
the coefficients of the terms on the right-hand side of the equation 
must equal the coefficient of the term on the left-hand side of the 
equation, 

14. Problems of interpolation between terms at unequal inter- 
vals can also be dealt with in a simple way by the formulas of 
Divided Differences (see Chapter VIII). 



CHAPTER V 

FINITE DIFFERENCES. CENTRAL DIFFERENCES 

1. It has already been stated that in interpolating between given 
values of a function the form of the expression connecting these 
values is assumed to be parabolic, and that this assumption is usually 
only an approximation to the truth. It remains therefore to be 
considered by what methods the best result can be obtained by the 
processes of Finite Differences. 

2. In developing the formulas of this chapter, it will be assumed 
that a number of equidistant values of the function are given. 
Let us assume further that it is desired to interpolate a value f(x) 
intermediate between /(0) and /(l). It is clear that our knowledge 
of the shape of the curve on which the points lie is increased if we 
are given values of the function lying on both sides of /(0), and 
that generally the best value of f(x) will be obtained, if a limited 
number of terms is to be used, when the required value occupies 
as nearly as possible a central position in regard to the terms used 
in the interpolation. 

The formulas of Central Differences are designed to give effect 
to these considerations. 

3. The more familiar formulas of Central Differences are as 
follows : 

Stirling's : 

. «(*-l) A'/(-l) + A'/(-2) ,,<(!^l) A4 ^_ ? x 

x (a? - 1) (a? - 4) A'/(- 2) + A'/(- 3) 

+ 5! 2 

+ *<*-l)(*-4) A . /( _ 8) + (1) . 
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Bessel's : 

, »(-!) A«/(-l) + A'/(0) («-j)«(«-l) 

2! 2 3! ^ ' 

(« + 1) « (« - 1) (« - 2) A«/(-l) + A</(- 2) 
+ 4! 2 

+ (g-j)(a» + l)«(«-l)(«-2) A „^ ( _ 2) 

, (g + 2)(g + l)«(g-l)(a?-2)(g-3) A«/(-2) + A«/(-3) 
+ 6! 2 "•"■■■ 

(2)- 

Gauss': 

/<•) -/(0) + • A/(0) + f^L> A'/(-l) + fe+WSzl) A »/(- 1) 

t («+ !)«(«-!)(«- 2) A ^ ( 

TO. 

+ (> + 2)(« + y-i)(,-2) Ay( _ 2)+ (3) 

Everett's : 

+ ^(o) + ^)A»/(-i) + ^^- 1 6 y- 4 ) Ay ( -2) 

+ y ^-l)( y ' 7 -4)(^-9) A8/( _ 3) + (4) 

[where y = 1 — x\. 

4. These formulas can be obtained in various ways from the 
ordinary formulas of advancing differences. Once, however, the 
scheme of differences entering into a formula is settled, the co- 
efficients can readily be calculated by the method of Separation 
of Symbols. An example may be given of the demonstration of 
Gauss' formula by this method. 
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Example 1. 
To express /(#) in terms of/(0), A/(0), A a /(- 1), A 8 /(- 1), .... 

Let /(*) « A Q f(0) + A x A/(0) + 4 2 A a /(- 1) + . . . . 
Then, since 

A '/(~l) = 1^/(0); A'/(-l)= r ^L/(0); etc, 

Multiplying up by (1 + A) 7 - 1 , and equating coefficients of A""" 1 , 

. _ (r+x-l)(r + a ;-2)...(x-r + l) 
■ **■* (Sr-l)l • 

And, multiplying up by (1 + A) r , and equating coefficients of A* 

A j.A _ (r + *)(r + x-l)...(x-r+ l) 
A "~ 1 + A » 2f! ' 

Hence, by subtraction, 

. _ (r+x -l)(r+x-2)...(a;- r ) 

Therefore f(x) =/(0) + «A/(0) 

+ i^ A . /( _ 1) + » + l)«(«-l) A . /( _ 1)+ .... 

The other formulas should be proved, in a similar way, as exercises 
by the student*, 

5. The formulas of central differences, although in a different 
form, are intimately associated with those of advancing differences. 
For example, if an interpolated value is calculated by using the 
first three terms of Stirling's formula, it is obvious that the values 

* See J. I. A, Vol. 50, pp. 28-33. 
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of /(— 1), /(0) and /(l) are brought into the calculation. It is 
easy to show that the result is identical with that obtained by 
using the first three terms of the advancing difference formula 
starting with the term/(- 1). 

It may be observed that the first two terms of Stirling's formula 
also involve three values of the function ; the third term merely 
introduces the correction necessary to make the formula true to 
the order of differences (i.e. the second) implied by the use of three 
terms of the series. Thus, as the (2r)th and the (2r + l)th terms 
of Stirling's formula both involve the use of (2r + 1) values of 
the function, there is ordinarily little advantage in using the extra 
(2r + l)th term in any calculation. 

Similarly in Bessel's formula no material increase in precision is 
gained by using 2r terms rather than 2r — 1 terms. 

Gauss' and Everett's formulas are each true to the order of 
differences involved and for general use they would appear to be 
the best of those propounded. 

6. In view of the remarks at the beginning of the foregoing 
article, it may well be asked what are the advantages of central 
difference formulas, as compared with advancing difference formulas 
so chosen as to make the interpolated term as nearly as possible 
the central term of those employed. It may at once be said that 
the theoretical advantages are small but that the practical ad- 
vantages may be considerable. Thus if it be desired to introduce 
further terms of the original series into the calculation, the 
original calculations relating to the central difference formulas 
hold good, and the values of fresh terms of the formula can be 
calculated until the desired degree of approximation is attained. 
If however an advancing difference formula is used, the introduc- 
tion of fresh terms of the original series, while retaining the 
interpolated term in a central position, necessitates the changing 
of the origin and the recalculation of all the terms of the formula. 
An example will make this point clear. 

7. Example 2. Required to interpolate the value of a unit 
accumulated for 17 years with compound interest at 5 per cent, 
per annum, having given the values for 0, 5, 10, ... 30 years. 

For central difference formulas we must take our origin at 15 
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years, and we will take 5 years as the unit. Thus we get the 
following scheme : \ 



No. of 

years 


X 


/W 


A 


A* 


A 1 * 


A 4 


A 5 





-3 


1- 


•27628 


07633 


•02110 


•00580 


•00165 


5 


-2 


1-27628 


•35261 


•09743 


•02690 


•00745 


•00206 


10 


-1 


1*62889 


•45004 


•12433 


•03435 


•00951 




15 





2-07893 


•57437 


-15868 


•04386 






20 


1 


265330 


•73305 


•20254 








25 


2 


3-38635 


•93559 










30 


3 


4-32194 













As the unit of time is 5 years, the required value is repre- 
sented by/('4). 

Central Differences. We will use Gauss' formula, i.e. 

/(•W(0)+*A/<(>) 



+ 



x{x — 1) 



A s /(-l) + 



(x + l)x(x-l) 



A«/(-l)+.... 



2! JK ' ' 3! 

The successive terms, with the corresponding values of /('4), are 

/(0) = 2-07893, 1st approx. = 2-07893 
•4A/(0)= -22975,2nd „ =2-30868 

'^=^A 3 /(-l)= --01492,3rd „ =2-29376 

l-4x-4x--6 frf^y = _ . 00192> 4th n = 2-29184 



1-4 x 4 x - -6 x - 1-6 

24 



A 4 /(-2)= -00017,5th „ =2-29201 



2-4xl-4x-4x--6x- l-6 
120 



A"/(-2)= -00002,6th 



= 2-29203 



It is clear that no further terms would affect the calculated 
value. The true value is 2-29202. 



H. T.B.I. 
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Advancing Differences. 
1st approximation /(0) = 2*07893, 

2nd approximation /(0) + '4 A/(0) = 2*30868, 

3rd approximation 

/(- 1) + 1-4 A/(- 1) + !2ii! A«/(- 1) = 2-29376, 
4th approximation 

/(-l) + l-4A/(-l)+i^A'/(-l) 

+ r4X 'g X "' 6 A'/(- 1) = 2-29184, 



5th approximation 



/(- 2) + 24 A/(- 2) -H^ii!! 4 A . /( _ 2) 



+ ^A>/(- 2 ) 

+ 2-4xl-4 2 X-4x-6 A4/( _ 2) = 2 . 29200> 



6th approximation 

. , 2*4xl-4x 

120 



2*4xl-4x*4x--6x-16 AK „, mN rtrtrtrtrtrt 
(5th approximation) + ^ A 5 /(- 2) = 2-29202. 



It will be observed that in proceeding to the 3rd and 5th ap- 
proximations using advancing differences every term in the formula 
has to be recalculated, whereas, in the application of the central 
difference formula, terms already calculated hold good whatever 
be the degree of approximation. 

It should be noted, however, that both formulas give mathe- 
matically the same results, the difference of a unit in the final 
figure being due to the use of only five places of decimals 
throughout. 

8. As regards other practical points, it may be observed that 
the numerical coefficients in the central difference formulas are 
smaller than those in the advancing difference formulas (see 
Example 2). 

Other advantages arise in special cases. Thus Bessel's formula 
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can conveniently be applied for the bisection of an interval, since 
the alternate terms vanish, giving 

/m - /(0)+/a) l A'/(-l) + A'/(0) 



8 + I | ^/(-*HA-/(-l> + „, (5 , 

Everett's formula gives the same value. 

9. It should be noted as regards Everett's formula, that in cal- 
culating a series of values the work is nearly halved since it will 
be found that terms in the formula can be made to do duty twice, 
tf #" terms reappearing in the calculation as "y" terms. 

This will be seen at once, for, 

+ y /(0) + ^f^A»/(-l) + ^y 3 -W-^) Ay( _ 2) + ,„ 



3! 



5! 



and 

/(i+y)=y/(2)+^f L) A'/a)+ y(y3 ~ 1) 5 ( f~ 4) ^/(0)+ ... 
+ . /( i )+ g^)AV(0) + ^- 1 ^- 4 > Ay(-iH.,. 

the last line being identical with the first. Thus, if we are inserting 
terms in a series by subdividing the interval into five equal parts, 
x = % *4, ... and y = '8, *6, .... Therefore half of the terms used 
in the calculation of /(*2) can be made to do duty in the calculation 
of /(1'8), and similarly for the other terms. 

10. An example will indicate the method of working. 

Example 3. Using Everett's formula, interpolate the missing 
terms in the following series, between /(40) and/(50). 



X 


/(*) 


A 


A 3 


A 3 


A 4 


30 


771 


91 


48 


36 


5 


35 


862 


139 


84 


41 


37 


40 


1001 


223 


125 


78 


28 


45 


1224 


348 


203 


106 




50 


1572 


551 


309 






55 


2123 


860 








60 


2983 











3—2 
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The coefficients of the several terms in Everett's formula are 
•2 - -032 -006336 

•4 - -056 010752 

•6 - -064 -011648 

•8 - -048 -008064 

The work may be arranged in tabular form: 



xf(l) 



x(x*-l) 
3! 



A 2 /(0) 



x{x*-l)(x*-4) 



51 



A 4 /(-l) 



Sum of first 
three terms 
(2) + (3) + (4) 



Sum of 
second 
three 
terms 



Interpolated 

result 

(5) + (6) 



(1) 
•2 
•4 
•6 
•8 

•2 
•4 
■6 
•8 

'2 
•4 
•6 

'8 



(2) 

200*2 
400-4 
600-6 
800-8 

244-8 
489-6 
734-4 
979-2 

314-4 

628-8 

943-2 

1257*6 



(3) 
• 2-6 

- 4-7 

- 6-4 

- 4-0 

. 40 

■ 7*0 

■ 8-0 
. 60 

■ 6*5 
-11-4 
■130 

■ 9-7 



o-o, 

01 
01 

o-o 

0-2 
0-4 
0-4 
0-3 

0-2 
0-3 
0-3 

0-2 



(5) 

197-6 
395-8 
595-3 
7968 

2410 
4830 
726-8 
973-5 

308*1 

617-7 

930-5 

1248-1 



(6) 



796-8 
595-3 
395-8 
197-6 

973-5 
726-8 
4830 
241-0 



(7) 



1037*8 
1078-3 
1122-6 
1171-1 

1281-6 
1344-5 
1413*5 
1489-1 



Columns (2), (3) and (4), which represent the first three terms of 
the formula, are obtained by ordinary multiplication. Column (5) 
gives the sum of these terms. From what has been said above, it 
is clear that column (6), which represents the second set of three 
terms of the formula to fourth central differences, is obtained by 
writing down, in reverse order, the values of column (5) applicable 
to the previous group of terms. The addition of columns (5) and 
(6) then gives the desired result. 

The given values of f(x) have been taken from the tabulated 
values of the probability of dying in a given year of age according 
to the H M mortality table, multiplied by 10 5 . 

The tabular values for the interpolated terms are 1038, 1081, 
1122, 1172, 1281, 1345, 1415, 1490. The small differences between 
these values and the interpolated values are due to the fact that 
the H M table was constructed by means of a mathematical formula 
which is only approximately represented by Everett's formula. 



CENTRAL DIFFERENCES 37 

11. Another method of applying the principles of central 
differences is to express the required function in terms of known 
values of the function among which it occupies a central position. 
This can conveniently be done by Lagrange's formula. The for- 
mulas are of two types according as the number of terms involved 
is odd or even. Thus we have by Lagrange : 

Number of terms 2n + 1. 

3-term formula, 

/(«) = /(-i) /(Q), /(i) ,« 

■*(x*-l) 2(0+1) ' r 2(a;-l) w " 

5-term formula, 

fix) _ /(-2) /(-l) /(0) /(l) /(2) 

0(0 s -l)(0 a -4) 24(a;+2) 6(x+ly 4>x 6 (x-iy 24, (x-2) 

(7). ' 

7 -term formula, 

f{x) _ /(-3) /(-2) /(-l) 

(0 s -l)(rf- 4) (0 s -9) 720(0 + 3) 120 (0 + 2) "^ 48 (0+I) 

360" r 48(0-l) 120(0-2) "•" 720(0-3) v J ' 

Number of terms 2n. 
4i-term formula, 

/(«) = /(-I) , /(-*) /(*) , /(§) 

(* -*)(**-*) 6(0 + f) + 2(0+i) 2(0-0 + 6(0-f) 

(9). 

/(«) ^ /(-i) , /(-i) 

(0 s -i)(0 a -|)(0 s --^) 120 (0 + f)" 1 " 24(0 + f) 

/(-*) , /(*) /(f) , /(#) a0 x 

12(0 + i)"'"l2(0-^) 24 (0-f)" 1 " 12O(0-f) "* ' 

12. These formulas, of course, yield identically the same results 
as other central difference formulas embracing the same terms. 
To illustrate this we will recalculate the value of /('4) in the 
example given in § 7. 
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Example 4. See Example 2. Seven terms are given, the formula 
will therefore be 





/(•4) 






/(- 
720 x 

/(I) 
48 x- 


•3) 

:3"4 

•6 


/(-2) 
120 x 2-4 

/(2) 
120 x- 16 


"*" 48 x 1-4 
+ 720 x- 2-6 


•4(-16- 


-1)(-16- 


4)(16- 

/(0) 
36 x -4 


9) 

i + 



or /(-4) = - -0046592/(- 3) + -0396032/(- 2) - -169728/(- 1) 
+ -792064/(0) + -396032/(1) - -0594048/(2) + -0060928/(3) 
-+ -05055 - -00466 
1-64665 -27647 
1-05079 -20117 
•02633 



+ 2-77432- -48230 
= 2-29202 as before. 
Note, as a check, that the algebraic sum of the coefficients of the 
terms on the right-hand side of the above equation is unity. 

13. For the sake of completeness it is necessary to refer to a 
system of notation in connection with central differences which 
was introduced by W. S. B. Woolhouse and is still in use to some 
extent. This system of notation is compared with that used in the 
previous chapters in the following scheme : 

Ordinary Notation Woolkoustfs Notation 

/(-2) /(-a) 

A/(-3) 
/(-l) A*/(-2) 

AA-1) A 3 /(-2) 

/(0) Ai/(-l) A</(-2) 

A/(0) A¥(-l) 

/(I) ' AS/(0) 

A/(l 
/(2) 

where A/(- 2), A*/(— 2), etc. are denoted by a_j, 6_i, etc. and 
a,, = £,(a_, + a +I ), 
c„ = |(c_! + c +1 ). 
Under this notation Stirling's formula to fourth differences is 

/(g)-/(0) + SB. + gfo + v 31 / Ci+- L 4| — 'd> 



/(-I) 


0_2 


6-1 


«-i 




/(0) 


(«o) 


&o 




* 


/(I) 


«ii 


*i 






/(a) 
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Similarly Gauss' formula can be written 

+ (* + l)*fa-l)(*-2)^ (12) 

14. Another system of notation, which is extensively used, is 
that due to W. F. Sheppard. Two operators 8 and jju are used, such 
that 

*/(- 1) =/(<>) r/(- 1). M/(i) = \ [/(0) +/(l)J 
«/<i) =/(i)- /(OX /*«/(0) = J [a/(« + s/(- 1)1 

etc. etc. 

This notation, although somewjiat complicated, gives the usual 
central difference formulas in very convenient forms. 



CHAPTER VI 

FINITE DIFFERENCES. INVERSE INTERPOLATION 

1. In direct interpolation a series of values of the function is 
given and the problem is to find the value of the function corre- 
sponding to some intermediate value of the argument. 

In Inverse Interpolation the problem is reversed and it is required 
to find the value of the argument which corresponds to some value 
of the function, intermediate between two tabulated values. 

2. In certain cases of mathematical functions the desired result 
can be obtained by direct calculation. Thus if 

J . = 1 °gy 

log a' 

and the value of x can be found equivalent to any given value 
ofy. 

Where this is not the case various methods can be adopted. 
These will be examined in order. 

3. Let y =/(#) be the given value. Then 

y=/(*)=/(0) + *A/(0) + ^=^A«/(0) + .... 

If it be assumed that the higher orders of differences vanish, and 
that the values of A, A 2 , etc. are obtained from the given terms of 
the series, then we have an equation in x which can be solved by 
the usual methods. 

The disadvantages of this plan are firstly that an equation of 
higher degree than the second is troublesome to solve, and secondly 
that for certain functions the degree of approximation may not 
be very close. Since a quadratic equation employs only three terms 
of the series, it often happens that no close approximation can be 
obtained. In all cases the intervals between terms should be as 
narrow as possible, so that accuracy may be increased and the use 
of higher orders of differences obviated as far as possible. 
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4. This difficulty of solving an equation in x of higher degree 
than the second can be overcome in two ways. Assume, for purposes 
of illustration, that four values of the function are given, viz./(0), 
/(l),/(2)and/(3). 

Then 

/(*) =/(0) + * A/(0) + X -^> Ay (0) + gfr-^K*- 2 ) A ./ (0). 

No further differences can be calculated and therefore, since fix) 
is known, the corresponding value of x is found by the solution of 
a cubic equation in x. The solution of the cubic can however be 
avoided by proceeding as follows : 
Taking three terms only at a time 

/(*) =/(<>) + * A/(0) + ^p- } A*/(0), 
and /(^)=/(l) + (^-l)A/(l) + (a! " 1) 2 < j a; ~ 2) A'/(l). 
The third difference error in the first equation is 

«(,-iK«-2) At/(0X 

and, in the second equation, 

If now both sides of these equations be multiplied respectively by 
(3 — x') and x (where x is a rough approximation to the required 
value, obtained by inspection) and the equations so weighted be 
added together, a new quadratic equation in x will be formed from 
which the third difference error will be practically eliminated. The 
work of solving a cubic equation has been avoided, but all terms 
have been used without sensible loss in accuracy. 

If the mere arithmetic mean of the equations were taken, with- 
out weighting as above, it is possible that, in certain cases, a worse 
result would be obtained by taking four terms instead of three. 

5. Alternatively the solution of the equation may be obtained 
by successive degrees of approximation. 

Thus, taking the above equation and neglecting differences of 
the second and higher orders, we obtain as a first approximation 
the value x u where 

. _ /(«)-/«» (1) . 

1 A/(0) ^ h 
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A further approximation is obtained by taking second differences 
into account and writing x x in place of x in the equation, thus giving 



x 2 = 



.(2). 



A/(0) + J<«i-l)A»/<0) 

When third differences are taken into account x 2 is written for x, 
giving 



^s = - 



.(3). 



A/(0) + * (a* - 1) A 2 /(0) + * (a* - 1) (* 2 - 2) A*/(0) 
These processes can be repeated until the desired degree of 
approximation is reached. The method has the disadvantage of 
being somewhat laborious. On the other hand it has the advantage 
that an error of calculation at an early stage does not vitiate the 
result, being rectified by the further approximations. 

6. A different method of procedure is to treat # as a function of 
f(x). Thus since 

we may write x = <j> (y). 

We therefore treat x as a function of y and, since the given 
values of y (i.e. /(0), /(l), etc.) will usually represent unequal 
intervals of the variable y, we must resort to interpolation by such 
a method as Lagrange's or Divided Differences, in order to obtain 
our value of x (i.e. </> (y)) corresponding to the given value of /(#). 

7. The following example is worked out in each of the above 
ways. 

Example. Find the number of which the log is 2 J, having given 

log 200 = 2-30103 

210 = 2*32222 

220 = 234242 

230 = 2-36173. 
Method 7. 



/(*) 


A 


A 3 


A 3 


2-30103 
2-32222 
2-34242 
2-36173 


•02119 
•02020 
•01931 


-•00099 
-•00089 


•00010 



INVERSE INTERPOLATION 43 

f(x) = 2-33333 = 2-30103 + -02119a; - -00099 g ,( g ~ 1 ) 

+ -00010 C(<, " 1 > (< '" 8) . 

6 

Or, by reduction, 

a? 3 - 32*7a? a + 13031 - 1938 = 0. 
Whence, solving the cubic, 

x = 1-5443. 

Since the initial value of x is 200 and the unit of measurement 
is 10, the result of the calculation is to give 215*443 as the required 
value of x. 

8. Method II. 

f(x) = 2-30103 + -02119a? - 00099 ?l£zD , 

also, f{x) = 2-32222 + -02020 (x - 1) - -00089 ( a? ~ 1 K a; " 2 ) , 

The first approximation to the value of x is x' = 1*5, so that 
3 — x =1*5. Since the values of 3 — x and a?' are approximately 
equal, we may take for our "weighted" equation the arithmetic 
mean of the above equations, giving 

/(*) = 2-33333 = 2-30108 + -02161a?- -00047^. 
• Whence, solving the quadratic, # = 1*5443 as before. 



9. Method III. 
1st approximation 



2-33333-2*30103 , „ nyfo 
* 1= -02119 =1 '° 243 ' 



233333 - 2*30103 _ 
^"•02119-i x -5243 x -00099" 6 > 



2nd approximation 

3rd approximation 

= _ 2*33333-2*30103 

Xz -02119 - 4 x -5432 x -00099 + 4 x -5432 x (-5432 - 1) x *00010 

= 1*5442, 

which differs only slightly from the value obtained by Methods I 
and II. 
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10. Method IV. Under this method we may consider the data 
to be as follows : 



X 


/(*) 


2*30103 
2*32222 
2*34242 
2*36173 


200 
210 
220 
230 



It is required to find the value of/(2'33333). By Lagrange, 

flVKUCm ( ' 01111) ( ~ '° 0909) ( " '° 2840) x200 
/ (2 66336) - ( _.Q2119K- -04189) (- '06070) X ZW 

(-03280) (- '00909) (- '02840) 
+ ('02119) (- -02020) (- -03951) 

(-03230) (-01111) (--02840) 
(-04139) (-02020) (- -01931) 
(-03230) (-01111) (--00909) 
"*" (06070) (-03951) (-01931) 

= - 10-7748 + 103-5416 + 1388764 - 16-2006 

= 215-443. 
This result agrees with the values previously found. The true 
value is 215*442 and it will be seen, therefore, that all of these 
methods give results in this case which are closely approximate to 
the true value. 

11. Methods I— III differ from Method IV, in that the first 
three methods assume that/(#) is a parabolic function of x> whereas 
Method IV assumes that x can be expressed as a parabolic function 
of /(#). Both these assumptions may be sufficiently accurate in 
many cases, but in other cases an inspection of the trend of the 
given values of the function may indicate which assumption is to 
be preferred. 

12. By the use of the formulas of Divided Differences given in 
Chapter VIII, Methods I — III can be applied to cases where the 
given values of the function are at unequal intervals. 

The application of Method IV is, of course, perfectly general. 



CHAPTER VII 

FINITE DIFFERENCES. SUMMATION OR INTEGRATION 

1. Summation is the process of finding the sum of any number 
of terms of a given series. This can be accomplished either if the 
law of the series is known or if a sufficient number of terms is given 
to enable the law to be ascertained. As will be shown in Chapter XX, 
if no mathematical law is apparent, methods can be applied by 
which the approximate sum of a series can be obtained. 

2. Consider a function F(x) whose 1st difference is f(x). Then 
we have 

t(\) -F(0) =/(0) 

F(2) -F(l) »/(l) 

F(a) -J»(o-l)-/(o-l) 
F(a + 1)-F(a) -/(a) 

F(n - 1) - F(n - 2) =/(n - 2) 
F(n) -F(n-l)=f(n-l). 
Summing both sides, we obtain 

F(n)-F(0)-f(Q)+f(l) + ...+/<• -1) 
or F(a) - F(0) =/(0) +/(1) + ... +/(« - 1) 

or F(n)-F(a)=f(a)+f(a + l)+...+f(n-l) ...(1). 

It is clear, therefore, that the sum of any number of terms of a 
series of values off(x) can be represented by the difference between 
two values of another function F {x) whose 1st difference is /(#). 
By analogy with the system of notation already adopted for ex- 
pressing orders of differences, the process of finding the function 
whose 1st difference is fix) may be denoted by A -1 f(x). It is 

n-l 

customary to express /(0) +/(1) + ... +/(n — 1) as 2 f(x), the 

o 

terms at which the summation is commenced and terminated 
(designated respectively the inferior and superior limits of sum- 
mation) being indicated iu the manner shown. 
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3. The process of finding the value of F{x) is known as Finite 
Integration and F(x) is called the Finite Integral of /(#). Where 
the limits of summation are known we obtain by summation of 
f{x) the Definite Integral of /(#); if the limits of summation are 
not expressed we obtain merely the Indefinite Integral oif(x). 

4. As stated above, in obtaining the indefinite integral of f(x) 
no point is specified at which the summation is to commence and, 
since an unknown number of terms of the series is included, it is 
necessary to include in the value of F(x) a constant term which is 
of unknown value. 

This constant vanishes in the case of definite integrals, since if 

$f(x) = F(x) + c, 

then V /(*) = \F (») + c] - [F (0) + c] 

o 

= F(n)-F(0). 

5. It is obviously always possible to find the first difference of 
any function, but it does not follow that every function can be 
integrated. The functions which can be integrated are limited in 
number and the process of integration rests largely on the ingenuity 
of the solver aided by such analogous forms as may be obtained by 
the formulas of finite differences. 

Thus we have A a x = (a - 1) a x , 

whence it is easily seen that 



a- 



X _ 



A a* 
a-1' 



and therefore, since the result of differencing — — is to give a x , 
we have 

%a x = =- + c, 

a — 1 

where c is the constant introduced by integration. 

The sum of the series a r + a r+l 4- ..I 4- a r+7lr ~ l is at once obtained 

r+n-l 

by finding the value of the definite integral 2 a x , which by § 2 

r 

a r * n a r 
is equal to = =- , which agrees with the familiar result for 

the sum of a geometrical progression. 
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6. Similarly since A#< m > = mx^-^, 
by analogy 
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^(m+l) 
2^ m) = r- + C. 



ra + 1 

7. The following table gives the values of some of the simpler 
integrals. They should be verified as an exercise by the student. 



Function 



Indefinite Integral 



x 
a* 

xi m ) 
xi~ m ) 



x(x-l) 

2 

a* 

a-l 



m + 1 
xK-™-*) 



+c 
+c 



+ c 



(ax+b)(ax-l + b)...(ax-m + l + b) 



(ax+b) (ax-l + b) ... (ax-m + b) 



+c 



(ax+b)(ax+l + b)...(ax+m-l + b) 
1 



(ax-l+b)(ax+b)...(ax+m-l + b) 
a (m+1) 
1 



+c 



(ax+b)(ax + l + b)...(ax + m-l + b) 
1 



(a# + &)(a#- 1 + &)...(«# - rn + 1 + 6) 



-a(m-l)(a#+&)(a#+l + &)...(aa7+m-2 + &) 
1 



+ c 



-a(m-l)(a#-l-f&)...(a#-wi + l4-&) 



8. If the form of the function is unknown, a general formula for 
the sum of a series of values may be obtained as follows, since 

/(•) -/(0) + * A/(0) + ^p) A'/ (0) 

+ ,(,-l)(,-2) At/(0) + 

Integrating both sides, we have 

V(c)-o + >/(0) + ^>A/(0) + '<*-y- a > Ay(0) + ..., 

or, integrating between limits, we have 

^) = n/(0 ) + W i^A/(0) + W(ra - 1 3 ) ! (w - 2 > Ay(0) + ... 
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or, more generally, 

,. . n(n — 1) . ,. . n(n— l)(n- 2) ..,. . 
= «/(«) + -^-jj — A/(a) + -i ^ '- A 2 /(a) + . . 

(2). 



It is instructive to obtain this result by the method of separation 
of symbols. For 

/(a + l) = (l + A)/(a), 
/(a+2) = (l + A)'/(a), 



Adding 

a+n-l 



/(a + n-l) = (l + A)-»/(a)L 



2 /(a) = {l + (l + A) + (l + A)»+...+(l + A)»T'}/(a) 
(14-A)"-1 



(1 + A)-1 



■/(«) 



1 h ^ 

-V(«)+'" j! ^A/(a)4. "< >t - 1 ><"- 2 > Ay(a) + .. 



9. An example will illustrate the use of the formula. 
Example 1. Find the sum of the first n terms of the series whose 
initial terms are 1, 8, 27, 64, 125. 
We have 



/w 


A 


A 2 


A 8 


A 4 


1 


7 


12 


6 





8 


19 


18 


6 




27 


37 


24 






64 


61 








125 
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Whence 

tt ; l , M - n(n-l) h n(n-l)(n-2) _ 
2 f(x) = nxl+ v 2j 7 xU- ^ -xl2 

n(n-l)(n-2)(n-3) 

4! 
_ n (n 8 + 2n 2 + rc ) __ n* (n + 1 ) a 
4 4 

which agrees with the formula for the sum of the cubes of the 
natural numbers. 

10. Where it is desired to integrate a function which is the 
product of two factors, the following device may often be utilised 
with advantage. 

Let the function be y = u x v x . 

Then ku x v x = u x+1 v x+1 - u x v x 

= %+l OWi - Vx) + % (tt*+i - u x ) 
= u x+1 Av x +v x &u x . 
Integrating both sides of the equation, we obtain 
u x v x = %u x+1 Av x + £v,Atf» 
or 2tfeAti»-ttptfe -£««+! Alfe (3). 

Thus, if the original function can be put in the form v x Au X) its 
integral can be made to depend upon that of u x+1 &v X) and, if the 
latter is in a form which can be readily integrated, the value of an 
apparently intractable integral may often be obtained in this way. 

Example 2. To find the value of Xxa x . 



Since 2a* = 


a® 
= , we may write 








%xa x = 2# 


Aa x 
a-1' 






Using the above formula (3) we 


get 








^ Aa x xa x 

zx - = = 

a— 1 a— 1 


a-1 


Ax 






xa x 
~a-l 


^ a x+l 
a-1 


, since 


A# = 1, 




xa x 
a-1 


a*+ 2 
(a-1) 2 


+ c. 





H. T. B, I. 
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11. Sometimes it may be necessary to apply the formula more 
than once in order to reduce the integral by stages to a standard 
form. The process is illustrated in the following example. 

Example 3. To find the value of %2 x x*. 

Remembering that A2* = 2 X > we may write 
22*# 8 = 2# 8 A2* 

= x*2 z -22 x + 1 Ax* 

= x*2 x - 22* +1 (3a? + 3x + 1). 

It will be observed that in applying formula (3), the degree of 
x y in the terms within the integral, has been reduced by unity. 
Proceeding as before we obtain 

22*«» = x*2 x - 2 (3x* + 3x + 1) A2 aj + 1 

= x*2 x - [2*+ l (3s 3 + Sx + 1) - 22*+ 2 A (3# 2 + 3x + 1)] 
= 2 X (x* - 6^ - 6x - 2) + 22*+ 2 (6x + 6) 
= 2*(^-6tf 8 -6tf-2)+2(6tf + 6)A2*+ 2 
= 2 X (x* - 6x> - 6x - 2) + [2*+ 2 (6x + 6) - 22 a! + 8 A (6x + 6)] 
= 2 X (x* - 6^ + 18a? 4- 22) - 22*+* x 6 
==2*(^-6# 2 + 18a + 22)-6x2*+ s + c 
-2*(rf-ft* + l$r-36) + a 
The above process is analogous to that of "Integration by Parts," 
which is dealt with in the Integral Calculus, Chapter XVIII, § 6. 



CHAPTER VIII 

FINITE DIFFERENCES. DIVIDED DIFFERENCES 

1. A simple method of interpolation is available, where the in- 
tervals between the given terms are unequal, by the method of 
Divided Differences. 

2. The application of the method rests upon the assumption, 
which, as has been shown, is the basis of all theorems for interpola- 
tion by means of Finite Differences, that f(x) is a rational integral 
function of x of the wth degree. 

On this assumption, it can be shown that f(x) can be expressed 
in the form 

A + A l (x - a^ + A 2 (x - a^)(x — a 2 ) + . . . 

+ A n (x - a x )(x- a 2 ) ... (x - a n ), 
where A , A 1} ... A n , o^, a 2 , ... a n are constants. 

3. In order to apply this formula in practice, it is convenient to 
introduce a scheme of notation on the following lines, where the 
symbol of operation is denoted by A' in order to distinguish it from 
the ordinary A. 



Value of x 


Value of 
Function 


A' 


A' 2 


A's 





/(0) 

/(«0 


/(«i>-/(0) 

/(«2)-/(«l) 

a 2 — a x 

/(%)-/ (« 2 ) 
c&3 — a 2 

/(«4)-/(08) 
«4 — a 3 


A'/(ai)-A'/(0) 

c& 2 

A'/(a2)-A'/(a,) 

a 3 -a x 

A'/^-A'/^) 

a 4"~ a 2 


A' 2 /(«i)"A'V(0) 
a 3 

A' 2 /(a 2 )-A'*/(ai) 
ct4 — ai 



etc. 



Generally, we have 

A V(Or) « A,n " 1 /(^)- A,n ' 1 /(^) < 



„ a n+r — a r 



4—2 
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Hence /(a,) =/(0) + a, A7(0), 
and /(oa) =/(c0 + (a* - <h) ^f^) 

=/(0) + a, A'/(0) + (a 2 - a,) { A'/(0) + a 2 A V (0)} 
=/(0) + a, A7(0) + a 2 (a, - a,) A V (0). 
By proceeding similarly for further terms, we find that we can 
write generally 
f{x) =/(0) + x A7(0) + « (* - oO A V(0) 

+ «(*- Oi) (* - a.) A '7(0) + . . . (1). 
This general form can be readily established by the method of 
induction. 

By giving appropriate values to a x , a 2 , ... the ordinary formulas 
applicable to equal intervals can be at once deduced. 



4. The general method of working will be shown more simply 
by an example. 

Example. Find the value of log 4*0180, having given the fol- 
lowing data: 



Number 


Logarithm 


4-0000 
4-0127 
4-0233 
4*0298 
4-0369 


•6020600 
•6034367 
•6045824 
•6052835 
•6060480 



Transposing the origin, we have the following scheme: 



X 


/to 


A' 


A' 2 


o-oooo 


•6020600 


*ssr , -«- 


- -000317 
•0233 


•0136 


0-0127 


•6034367 


*sr-™» 


- -000223 
•0171 


•0130 


0-0233 


•6045824 


•°T"-™« 


- -000186 
•0136 


•0137 


0-0298 


•6052835 


T,r-»'«™ 






0-0369 


•6060480 
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We have to find/(*0180), which is, by the above formula, 
/(0) + -0180A'/(0) + -0180 x '0053 A' 3 / (0) 
[the further terms will not affect the seventh place of decimals], 
where /(0) = '6020600, A'/(0) = '108402, A /2 /(0) = - 0136. 

Thus log 4-0180 = -6020600 + '00195124 - '00000130 = '6040099 
to seven decimal places, which agrees exactly with the true result. 



CHAPTER IX 

FINITE DIFFERENCES. FUNCTIONS OF TWO VARIABLES 

1. Questions involving functions of two variables arise frequently 
in actuarial practice. Thus the tabulated values of functions (e.g. 
annuities) dependent upon two lives may be given only for com- 
binations of quinquennial ages in order to economise space. If the 
value corresponding to any other combination of ages is required, 
resort must be had to methods of interpolation. 

2. In considering the problem of the changes induced in the 
value of /(#, y) by finite changes in the values of x and y we must 
consider x and y as being independent of each other. Clearly, if 
y were a function of x the expression f(x, y) could be made to 
assume the form of a function of x alone and it could be dealt 
with by the methods already developed in previous chapters. 

Thus x may vary while y remains constant, so that, if x changes 
to x+ h t the value of the function becomes f(x + h, y); or y can 
vary while x remains constant, giving a value f(x y y + k); or both 
x and y can vary independently, giving a value for the function of 
/(# + *,* + *). 

3. We shall proceed first to discuss the problem where the 
values of the function are given for combinations of successive 
equidistant values of x and y. 

Thus we may have 

/(#> V) A* + K y) fix + 2h,y) ... f(x + mh, y) 

f(x,y + k) /(x + h.y + k) f(x+2h ) y + k)^ ...f(x+mh,y + k) 

f(x,y + nk) f(x+h,y+nk) f(x+2h,y + nk) ...f(x + mh,y+nk) 

As has already been seen in the case of functions of one variable 
(Chapter III, § 6), this scheme can be simplified, for the origin can 
be placed at the point (x, y), and the unit of measurement can be 
taken as h in the case of the variable x and h in the case of the 
variable y. 
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The scheme then becomes 

/(0,0) /(1,0) /(2,0) /(m,0) 

/(0,1) /(1,1) /(2,1) f(m,l) 

/(0,n) /(1,») f(2,n) f(m,n) 

4. Since x and y may vary independently, a fresh scheme of 
notation must be introduced to express the variations which may 
arise. Thus A^ will be used to denote the operation of differencing 
with respect to x, y remaining constant, a corresponding significance 
attaching to A y , so that 

A./(0,0)«/(l,0)-/(<>,0), 
A,/(0,0)-/(0,l)-/(0,0), 

or, using the method of separation of symbols, 

/(l,0)«(l + A„)/(0,0). 

Accordingly, we have 
f(m, h) - (1 + AJ» (1 + A,)»/(0, 0) 

= (l + mA. + (*)A«.+ ...)(l+»A I , + (;)AV+...)/(0,0) 
=f(0,0) + mA x f(0,0)+(f)A\ /(0, 0) + (*) AV(0,0)+ ... 

+nA J ,/(0,0) + mn A a A y /(0 ) 0)+ (~) nA».A,/(0,0)+... 
+ (g) A» y /(0, 0) + m g) A. AV/(0,0) + . . . 

+ (")A'„/(0 ) 0) + ... 

(I)- 

Here A 2 ^, A 3 ^, ... can be written down by differencing the rows of 
the table of the function; similarly A 2 y , A 3 y , ... are the differences 
of the columns of the table. 
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To find AsAj,, A'^Aj,, etc. we have 

A* A,/(0, 0) - A, [/(0, 1) -/(0, 0)] 

=/(l, l)-/(0, 1)-/(1, 0)+/(0, 0), 
A\A y /(0, 0) - A»,[/(0, 1) -/(0, 0)] 

=/(2, 1) - 2/(1, 1) +/(0, 1) -/(2, 0) 

+ 2/(l,0)-/(0,0), 
and so on. 

Example 1. Table XVI of the "Short Collection of Actuarial 
Tables." To find -4 3 } :63l having given 

^30:60 = -H669. ^sJ !6 o = -13190, AJ !60 = -15494, 

4 3 U = -09809, ^=65 = -H039, 

A^, n = -07812. 
Here ra = -J, n = f , and 

A a! = -01521, AV*" -00783, A^A^- -00291, 

A„ = - -01860, AV = - -00137. 
Whence jig},,, = -10776, the correct value being -10773. 

5. An obvious method of procedure involving only first differences 
is as follows. Obtain the value of f(0, n) by interpolation between 
the values of /(0, 0) and /(0, 1). Similarly, obtain the value of 
/(l, n) from the values of/(l, 0) and /(1, 1). Finally find/(m, n) 
by interpolation between/(0, n) and /(l, n). Thus 

/(0, n) = r^T/(0, 0) + nf(0, 1), 

/(l, n) = 1^/(1,0) + n/(l,l), 



f(m, n)= l-m/(0, n)+ m/(l, «) 

=/(0, 0) + m A x /(0, 0) + nA y / (0, 0) + m?iA * A y /(0, 0) 

(2) 

= 1 - m 1 - n/(0, 0) + n . 1 - mf (0, 1) 

+ m.T^n/(l, 0) + mn/(l, 1) ...(3). 
Employing this formula in the example given above we find 

^63 = -10822. 
The method is suitable if only a rough approximation is required, 
but cannot be depended upon to give an accurate value. 
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6. Obviously the method can be extended by taking higher 
orders of differences. The disadvantage of this procedure is that 
it involves the calculation of further values of the function corre- 
sponding to a given value of x as a preliminary to applying the 
interpolation formula to find the value of f(%, y). The process 
thus becomes laborious and moreover we do not necessarily obtain 
identical values for f(x f y) if we interpolate first with regard to x 
and then with regard to y, or vice versa. 

7. As in the case of functions of one variable, we shall expect to 
obtain the best results when the principles of central differences 
are applied, i.e. when the required term occupies as tiearly as pos- 
sible a central position among the terms employed in the formula. 
The difficulty is that, in dealing with functions of two variables, 
we cannot adapt our formulas to any system of values which may 
be given. Thus an inspection of the advancing difference formula 
(1) shows that it involves points whose coordinates form a triangu- 
lar plan which may be illustrated thus: 



(0,2) 

O 

(0, 1) (1, 1) 
o o 

(0, 0) (I, 0) (2, 0) 

o o o 

This illustrates the formula where two orders of differences are 
taken into account, the black dot representing the interpolated 
term. It will be seen that the scheme is hardly satisfactory from 
the point of view of central differences. For most practical pur- 
poses, however, where ordinary actuarial functions are involved, 
formula (1) will give satisfactory results. 

8. Formulas embodying the principles of central differences can 
conveniently be obtained by an adaptation of Lagrange's for- 
mula. This formula applied to functions of two variables has not 
the same wide application as the ordinary formula of Lagrange 
previously given in Chapter IV, but, as will be seen below, it gives 
expressions for f(x, y) in terms of the neighbouring values. 
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9. General formula for 4 points. 

Taking all combinations of two terms except those which give 
rise to x 2 and y\ let 

f(x,y) = A(x-/B)(y-b) + B(x-j3)(y-a) + C(x-cL)(y-a) 

+ D(x-*)<y-b) t 
then /(a, a) = A (a - /3) (a - 6), 

/(«,&)-*(t-jB)(*-aX 

/(Aa) = D(^-a)(a-J), 
f(0 } b) = O(!3-a)(b-a), 
whence, substituting for J., B, G and D in the original formula, 

j \ > yj j \ > / ( a _ ^ ( a _ j) ./ \ ' (a - £) (6 - a) 

. /V£ 6 x (*-")(y-a) . *,* , (^-«)(y-6) 
+/W 6 ) (i 9- a)( 6^ a ) + /^ a >>( i 8_ a )( a _ & )-W- 

10. General formula for 6 points. 
Taking all combinations of two terms, let 

f(x >y ) = A(x-/3)(y-b) + B(x-l3)(y-a)+C(x-a)(y-a) 
+ D(x-a)(y-b) + E(x-a)(x-P) + F(y-a)(y-b). 

Taking the points a : a, a : 6, a : c, # : a, : 6 and 7 : a, and pro- 
ceeding as before, we arrive at the result 

fir ,A - f(a a\ {<£zBMzh + («-«)(«-/8) , (y-a)(y-6) | 
•' t * ,y, ~ /W ° ; t(/8-a)(6-o) (/9-a) ( 7 - a ) + (b-a) (c - a) J 

f(a Ay-odiy-b) 



/^ ,A f («~«)(y-ft) , («-a)(«-/8) | 
"/ W> a > l0 g _ a)(6 - a) + 08 - a) ( 7 - 8)] 



+f(B n (*- g )(y- o > 

4.// w „\ («- g )( g -fl) /« 

+/(%a) (7-«)(7-^) ^ 
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11. General formula for 9 points. 

Taking all combinations of two terms, each involving x, with 
two terms each involving y, let 

f{<c,y) = A{x-P){x-y){ y -b){y-c) 
+ B(x-P)(x-v)(y-a)(y-c) + C(x-/3)(x-y)(y-a)(y-b) 
+ B(x-a)(x-y)(y-b)(y-c) + E(x-a)(x-y)(y~a)(y-c) 
+ F(x-a)(x-y)(y-a)(y-b) + G(x-a)(x-/3)(y-b)(y-c) 
+ H(a:-a)(a:-0)(y-a,)(y-c)+I(x-a)(x-P)(y-a)(y-b). 

Whence, proceeding as before, we have 
fix v) - f(a a) g -ffX*-7)(y-&)(y-c) 

, f(n M (a?-j8)(a;-7)(y-a)(y-c) 
+/1* °> ( a _ £) (« - y) (6 - a) (6 - c) 

. y-/„ ^ 0-ft)Qc-7)(y-a)(y-6) 
+/ W C; («-/3)(a- 7 )(c-a)(c-6) 

, ,vo „>, Qc-«)Qc-7)(y-o)(y-c) 

+/ w«; (/S _ a) (/S _ 7) (a _ b)(a _ c) 

, f(R h , (a>-«)(x-y)(y-a)(y-c) 
/KP ' ' 09 - «) 08 - 7)(b -a)(b- c) 

4. vya ^ («-«)(«-7)(y-a)(y-6) 

/W J (^-«)(j3-7)(«-i)(c4) 

. *, _% («-«)(«-fl)(y-6)(y-c) 
A7> ; (7-«)(7-£)(a- )( a -c) 

4. /y~ R\ ( g - g K g -^)(y- a )(y-c) 

+/Wff; (7-«)(7-/3)(6-a)(6-c) 

+ /■/„ c \ («-«)(«-£)(y-<»)(y-&) / 6 \ 

+ M7 '^(7-a)(7- / S)(c-a)( C -6) W " 

12. Formula (4) is the general formula corresponding to the 
method of § 5; by altering the notation the identity of the two 
formulas (3) and (4) is apparent. 

Formula (5) includes six values of the function, the co-ordinates 
being related in the manner shown. It will be seen that the formula 
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can be applied to any of the following groups of values, the black 
dot representing the interpolated value: 



(0,2) 
o 








(0,1) 

O 


(1,1) 

O 


(0,1) 
o 


(0,1) 

O 


(i,i) 

O 




(-1,0) 

O 


(0,0) 

O 


(1,0) 
o 


(-1,0) (0,0) (1,0) 
o o o 


(0,0) 

O 


(1,0) 
o 

(i) 


(2, 0) l 

O 


P,-l) 
o 

(ii) 




(-1,-1) (0,-1) 
o o 

(iii) 




(-l ; 


,0) 

O 


(0,1) (1,1) 

O O 

(0,0) (1,0) 

O O 

(i,-i) 

O 




(-1,0) 
o 

(-1,-1) 

o 


(1,1) 

O 

(0,0) (1,0) 

O O 

("i, -i) 

o 



(iv) (V) 

Other systems could be written down, but the above are suf- 
ficient for purposes of illustration. 

System (i) is obviously the same as the advancing difference 
formula. It is obtained by writing a = 0, £=1,7 = 2, a = 0, 6 = 1, 
c = 2 in formula (5). 

Example 2. Find -4 3 i :63 from the data of Example 1, using 
the Lagrange formula applicable to system (i). It will be found 
that the same result is obtained. 

System (ii) should be expected to give a formula which will be 
more accurate than the advancing difference formula, since the 
interpolated term will occupy a more central position in relation 
to the terms employed. The formula is obtained by putting a = 0, 
£ = 1, ry = — 1, a = 0, 6 = 1, c = — 1 in formula (5). Working on 
the same example as before and taking the origin at the point 
(30, 60), the values of the function entering into the formula are 

^25:60 = '10080, Aj> :65 = -09809, 4 3 J !fl0 = -11669, 

A£ sK = -14006, A£ t 65 = '11039, ^60 = '13190. 

We obtain as a result A 3 \ :e3 = t 10770. The degree of approxi- 
mation, though close, has not, in this instance, been improved. 
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System (iii) is an inversion of system (ii) and should be useful 
for interpolation where x and y have negative values. 

The lack of symmetry of systems (iv) and (v) suggests that they 
are not likely to yield good results in practice. 

13. When nine points are used, as in formula (6), the system is 
represented by the following diagram : 



<- 


-1, ]) (0, 1) 

O O 


(1,1) 


(- 


■1,0) (0,0). 

O O 


(1,0) 
o 


-] 


, -1) (0,-1) (1,-1) 
o o o 



It will be seen that this scheme embodies all the principles of 
central differences and should therefore give good results. 

Taking the previous example with the origin at the point (30, 60) 
the six values entering into the formula for system (ii) are used 
together with the following additional values : 

A£ i65 = '08435, ^ 2 J :66 = -12132, ^ :66 == "15972. 

Making use of formula (6) the interpolated value is found to be 
•10771, a slightly better approximation to the true value than 
those obtained previously. 

On general reasoning we should expect a somewhat better result 
by taking the origin at the point (30, 65) so that the interpolated 
value would occupy a more central position. The values A^.^, 
-4gJ :65 , A gj.55 entering into the immediately preceding calculation 
are excluded, and the following values introduced: 

^2J:to= '06642, ^3J :70 = -07812, ^ :70 = -08756. 
On working out the result, however, we arrive at the value -10848, 
which is a worse approximation than the value obtained Jby the 
rough method of § 5. 

14. This apparent inconsistency illustrates one of the chief 
difficulties of interpolating between functions of two variables, 
namely, that one does not necessarily obtain a better degree of 
approximation by proceeding to a higher order of differences or by 
employing more terms in a formula. Changes in the value of f(x, y) 
occasioned by alterations in the values of x and y may be so con- 
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siderable that distant terms may have such a disturbing effect 
upon the formula used as to upset the agreement between the 
approximate interpolation surface and the true surface which 
represents /(#, y). 

It is thus difficult to say what will be the degree of approxima- 
tion of a given formula, but an inspection of the course of the 
differences will be some guide as to the advisability of introducing 
further terms into the calculation. 

15. Other devices may sometimes be adopted which enable the 
interpolation to be reduced to the work of a single variable inter- 
polation. 

Thus, if the sum of % and y is a multiple of 5, by suitably 
selecting the origin we may write 

f(x, - x) =/(0, 0) + x [/(l, - 1) -/(0, 0)] 



+ (*)[/(2,-2)-2/(l,-l)+/(0,0)] 



x— 1 x—2 



f(0,0)-x.x-2f(l,-l) 



+ ^FV(2,-2) (7). 

By referring to the point diagrams on previous pages it will be 
seen that the process is equivalent to interpolating along a diagonal 
line running through the various points. The formula is of the 
advancing difference type; the corresponding central difference 
formula would preferably be employed in practice. 

J. Spencer has given (J.I.A. Vol. 40, pp. 296-301) examples of 
the use of several ingenious methods of this character. 



CHAPTER X 

DIFFERENTIAL CALCULUS. ELEMENTARY CON- 
CEPTIONS AND DEFINITIONS 

1. In the subject of Finite Differences we were concerned with 
the changes in the value of a function consequent upon finite changes 
in the value of the independent variable. In the Differential Calculus 
we consider the relation of Ay to Ax when the value of Ax is made 
indefinitely small. 

The application of the Differential Calculus is largely limited to 
such values of a function as are finite and continuous, and, unless 
otherwise stated, this limitation is to be implied in the following 
demonstrations. In practice these conditions are almost universally 
fulfilled by functions entering into actuarial calculations. 

2. Let y — f{x) and let x receive an increment h. Then the 

change in the value of y is measured by f(x + h)-f(x) and the 

f(x + h)—f(x) 
rate of change of y is — ^ — J-±J. , The limit of this expression 

when h -*- is called the Differential Coefficient or First Derived 
Function of f(x) with respect to x. 

The operation of obtaining this limit is called differentiating /(#). 

Using the notation of Finite Differences the differential co- 
efficient becomes 

Lt $L 

and is variously denoted by -r- , f (x\ f(x), ■ ,— , Df{x). 

The symbol ~- or its equivalent represents an operation of the 

character described; the elements dy 
and dx must not be regarded as separate 
small quantities. 

3. The geometrical representation of 
the differential coefficient is illustrated 
in the accompanying diagram. T M N X 

Let the curve shown represent the 



p y 



«/y 




R 
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equation y =/(#). Let OM=x and ON=w-\-h t and let PM and 
QN be the corresponding ordinates. Let PR be the perpendicular 
from P on QN and let QP be produced to cut OX at T. 

Then /(^^)-/^) ^ ^-^ = || = ^ =tan p m 
A MN PR TM 

When the point Q moves up to, and ultimately coincides with the 
point P, the line QPT becomes the tangent to the curve at the 

point P. The limiting value of — ^ — ^_i_jr j s therefore the 

tangent of, the angle which the tangent to the curve at the point 
(#, y) makes with the axis of oc. 



CHAPTER XI 

DIFFERENTIAL CALCULUS. STANDARD FORMS. 
PARTIAL DIFFERENTIATION 

1. The differential coefficient of any particular function can, of 
course, be obtained by direct calculation, but the process can 
usually be simplified by the application of the following general 
rules. The general similarity to the propositions already demon- 
strated for Finite Differences will be apparent. 

I. The differential coefficient of any constant term is zero. 

This is evident since a constant is a quantity which does not 
change in value in any mathematical operation. 

II. The differential coefficient of the. product of a constant and 
a function of x is equal to the product of the constant and of the 
differential coefficient of the function. 

m, d r £ , x , T , cf(x+h)-cf(x) 

Thus -r-[c. A (#)]== Lt — r ' L ^ J 

dx L J v /J a-a.0 h 

= Ltc /w)-/(«) 

-^ w 

III. The differential coefficient of the algebraic sum of a number 
of functions of x is the sum of the differential coefficients of the 
several functions. 

Let y*=u + v + w-\-..., where u, v,w, ... are functions of x, 

then Ay = Au + Av -f Aw -f . . . 

, Ay _ Au Av Aw 

Ax ~~ Ax Ax Ax '"' 

which, by proceeding to the limit, becomes 

dy _ du dv dw 9 . 

dx dx dx dx '" ^ 

IV. The differential coefficient of the product of two functions is 
the sum of the products of each function and the differential coefficient 
of the other. 

H. T.B.I. 5 
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Let y = uv } 

where u and v are both functions of a?, . J, 
Then Ay = (u + Aw) (u + Av) — wv 

= uAv + *yAtt + AuAv 

= w Av + (v + Av) Au 
l Ay Aw , x Aw 

and a-— s +(,+A,) S' 

whence, taking the limit, when v + Av -* v, 

dy_ dv c?w . 

da? da; di» ^ '' 

which may be written 

y da; w dx v dx ^ '* 

This result may be extended to include the product of any number 
of functions. 

For if y = uvw ; let vw = #, then y = uz. 

Whence 1^ = 1^ + 1^. 

-r, A 1 ^ 1 (in 1 dw 

But ^T-=-T- + -^-- 

Therefore 1^ = 1^ + 1^ + 1^ ( 5 ). 

y dx u dx v dx w dx 7 

Multiplying by uvw y we obtain 

dy du , do dw /m . 

-r- = vw-r- +wu-j-+ uv j— (6). 

dx dx dx dx 

Similarly for the product of any number of functions. 

V. The differential coefficient of the quotient of two functions is 
(Diff. Goeff. ofNumr.)(Denr.)-(Diff. Coeff. of Denr.) (Numr.) 
' Square of Denominator 

Let y — -- 

u + Aw u vAu — uAv 
Then . Ay = 



and 



v + Av v v (v + Av) 

Au , Av 
Ay.^. Aa; . Aa? 
Aa; ~~ v {v + Av) ' 
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whence, taking the limit, 

du dv 

dy __ dx dx .,.. 

S^~ v> Uh 

which may be written 

1 dy 1 du 1 dv /QX 

ss-saria ^: (s> 

VI. TAe differential coefficient of y with respect to x, where y is 
a function of u and u is a function of x, is the product of the dif- 
ferential coefficients of y with respect to u and u with respect to x. 

For ^/-^ Aw 

Ax Au' Ax' 

whence, taking the limit, 

dy^dydu 

dx du' dx ^ 

Si »«r 2-&S-5 ' ; <10X 

and so for any number of functions. 

2. Various standard forms can now be developed, mainly from 
first principles. It is instructive to note the points of analogy with 
the corresponding forms for Finite Differences. 

(i) y = x n , 

ty _ t (x+h) n -x n 
dx h ^ h ' 



/!" 



x" 
= Lt — 



h\ n 



(v3 



Expanding by the Binomial we have 

.U^,[l + ^ + ...] 

= nx n r 1 . v 



5—2 
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(ii) y = a x , 

dy _ T a x + h - a* 
dx k-*>o A 

= a*Lt — j — 
A-*.o n 

= a*Lt ^ fl + h log, a + | 2 (log, a) 3 + ... - ll 

^a x Lt\\og e a+^(log e af^..A 

= a x log, a. 
If ^T^' ^ = #log,e = *. 

(iii) y = log a a;. 

Then a^ = a;. 

But d(gy)^d(a?) dy ^ 

dx dy * dx * 

Hence, using the result established in (ii) and remembering that 

a? = x,we have 

dx „ , dv 

or 1 = a; log* a . -f- • 



Whence 



dy 



dx x log« a ' 

re , dy 1 1 

If y = log e *. ^ = ^ i - r -. 

(iv) 9"i/(")T\ y = e M ; f-log./(«). 

If *-[/(•)]•. 

dy = d[f(*)TT <¥(x) 

dx df{x) ' da; 
-*{/(*)TV'(#). 

Similarly if y = e /**, ^ = «/« . /' (a), 
and if *-**/«. |=^g. 
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(v) 



y = sin x, 
dy _ T sin (x + h) — sin x 

dx ~~ h + { 



h-*-0 



= Lt- 



2 sin ^ . cos 



(* + s) 



T. Sm 2 / *\ 



(vi) 
Then 

whence 
Therefore 



= cos x, 
y = sin"* 1 x. 



sin y = x, 
dx 



dy 



= cos y = Vl — sin 2 y = Vi — #*. 

d# da; Vl - a? 2 * 
3y 



3. The values of differential coefficients for the other trigono- 
metrical functions can be found by methods similar to those employed 
in (v) and (vi). The results are given in the table below and should 
be verified as an exercise by the student. 



Function 


Differential 
Coefficient 


Function 


Differential 
Coefficient 


X* 

a* 

e* 

log a x 

log e x 

sin x 
cos a? 
tana? 
cot X 
secx 
cosec x 


nx n ~ x 
a x \og e a 

1 


sin ~ l x 

COS"" 1 J? 

tan -1 # 

cot -1 X 

sec" 1 a; 

cosec* 1 x 


1 


VI- X 2 

1 


Vi-* 2 

1 
l+# 2 

1 

1 


x log e a 
1 

X 

cos X 

-sin# 

sec 2 .r 

— cosec 2 x 

sec x . tan x 

- cosec x . cot x 


W^-i 

1 


xjx 2 -! 



70 



DIFFERENTIAL CALCULUS 



4. Logarithmic Differentiation. This method is of special value 
in two cases. Thus if y — uvw . . . , where u, v f w, ... are functions 
of x, then 

log y = log u -r log v + log w + . . . , 



and 



1 dy _ 1 du 1 dv 1 dw 
y dx~ udx v dx w dx 



or 



dy Tl du 1 dv 

-f: = uvw |- — +- — + — — + 



■••]■ 



1 dw 
dx [u ax v ax w dx 

a result which agrees with that already obtained in § 1. 
Secondly if y = u v , u and v both being functions of x, 

\ogy = vlogu, 



and 



or 



1 dy _v du , 
y dx u dx ° 



da; 



dy ^du . dv 

ax dx ° dx 



5. We will now give some miscellaneous examples of differen- 
tiation. 

r , a + x 2 

(1) . *~tt- 

By the ordinary rule for a quotient 



dy 
dx 



(&+{*) 



dx 



-(a + a?) 



d{b + x) 
dx 



(ii) 



(6 + a;)' 

_ (b + x)2x-{a + x' i ) _ a?+2bx-a 
(b + x)* ~ (b + xf ' 

il-i 7 



y- 



l-5x + 6a?' . 



This can best be treated by resolving the expression into partial 
fractions. Then 



y = 



4 



1 - 3a; 1 - 2x ' 



dy _ 1 d (1 ~ ac) ~" 1 ^C 1 - 8 ^) 5 <*( 1 - 2g ) ri d C 1 - 2a; ) 

da; d (1 — 3a;) " dx d (1 — 2a;) " da; 



12 



10 



.(1-3*0'.,.. (l.-2a;) s 
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(iii) y = Va + x. . . 

1 
dy = d(a+.x) n d(a + x) _l ±-i 

dx d(a + x) ' dx rA } 

(iv) y = log (log x). 

dy _ d log (log x) d (log x) _ 1 
dx d (logic) * dx ~~#log#" 

(v) y = tan -1 



V^-l' 



^ = V Va; 2 -!/ Va 2 -1 d (a? 2 - 1) 

dx~ 1 ~'d(# 2 -l)* cfo 

* — .-i(^_i)-S.2^ 



1 + 



# 2 -l 
1 



^ Va? - 1 ' 



logy = - log « - 2 -- log - 3, 

*l> * . X 



l^_ld(log£-2) £(0-2) .. — -»■ 1 
yds d(0-2) * H + « , «*- z - „- 

I d (log 0^3) d(g-S) w — -o 1 
~0 d(0-3) * <fe iog*-d.-^ 

1 _1 J. . 0-2 

~x(x-2) -0(0-3) * J °^0-3* 

! I 

Whence ^ = -(^— rJ ~7 — =7>w 5\ + ^ lo S —~5 • ' * 

dx \x — 3/ |_0 (0 - 2) (x - 3) s ° - 3J 
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(vii) y^af.b**. 

log y = x log a 4- c x log 6, 

ydl = 1 ° ga + Ca:l0gcl0g6, 

-^ = a* . b** (log a + c* log c log 6). 

l 

(viii) y^aP + x*. 

In this case logarithmic differentiation must be used, but for this 
purpose the two terms must be taken separately. 

i 

Let af = u and x* = v. 

Then dy = du^d» 

dx dx dx' 

log U = X log Xy 

1 du 1 , 

- -5- = # . - + log a?, 

^ = tf*(l + l0gtf). 

Also log v — - log x, 

X 

v dx x'x g * as 2 ' 

dv ^l/ii \ 
^ = **.-(l-log*). 

Therefore -^ = #* (1 + log x) + x*~* (1 - log x). 

(ix) Differentiate log e x with regard to x 2 . 
Let 2/ = logg x and # = x\ 

Then ^^.^a^.-L. 

cte dx' dz dx' dz* 

dx 

Therefore ^=- — =— . 

xnereiore ^ ^.^ 2 ^. 
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6. In dealing with cases where a function of two variables is 

involved it is convenient to adopt methods similar to those used 

in Finite Differences (see Chapter IX, § 4). Thus we define Partial 

Differentiation as the process of differentiating a function of several 

variables with reference to any one of them, treating the other 

variables as constants. 

7) 7) 
This process is denoted by the symbols ^- , ^- , etc. We will also 

use the symbol Bx to denote a small change in the value of x. 

Let u=f{x,y). 

Then u + Bu=f(x + h, y + k) 

and 

&*«/(* + A, y + k)-f{m t y) 

__ / (x + h, y + k)-f Q, y + k) t f( x,y + k)-f(x,y) 
~ I ^ + I •* 

Proceeding to the limit when h and k successively -*- 0, 
f(x + h > y + k)-f(x,y + k) _ d 

h -dx J ^> y + L) 

™faf( m > V)> since 4^0, 

Therefore Bu = =- . Bx + — . By, 

ox ay 

, Bu du du By 

and ^r- = r- + — . c^ 



or, when Bx -*■ 0, 



Bx dx dy'Bx* 

du _du du dy ,,-.>, 

dx dx dy ' dx 



7. If/(«,y) = 0,g = 0and 



Whence 



n _du du dy 
dx dy' dx* 
du 
dy^__dx ,..(12). 



dy 
Example. If x 2 + xy + y* = 0, 

then ^ = 2a! + 2/' ^-* + 2 ^ 

and •/• = -^f. 

da? x+2y 



CHAPTER XII 

DIFFERENTIAL CALCULUS. SUCCESSIVE 
DIFFERENTIATION 

1. When ■—- is differentiated with respect to x, we obtain a 

further function of x which is known as the second differential 
coefficient or second derived function of y. This operation represented 

by -y-i-f-J is generally written -t-j or /"(#). If the function is 

differentiated n times with respect to x, the result is called the nth 

d n v 
differential coefficient or nth derived function and is written -r~ . 

Other symbols for the nth differential coefficient are/^ (x), D n y, y n . 

2. In many cases the value of the nth differential coefficient can 
be found readily by inductive reasoning. 

Example 1. 

y = log(x + a). 

and by analogy 

Example 2. 

y = a + bx+ cx*+ dx* + ...+ kx n , 
y x = b + 2cx+ 3dx 2 +...+n.kx n ~ 1 i 

y f - 2c + 2.3cta + ...+ri(rc-l) &#"-*, 



2/n = 

Example 3. 



nl*. 



= ?g-l 
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Making use of the method of partial fractions, 
4 5 






l-3a; l-2#' 
4.3 5.2 



(l-Zxf (l.-2#) 2 ' 

_ 4.3 n .n! 5.2 n .n! 
yn ""(l-3^ +1 (l-2^+ 1< 

3. Leibnitz's Theorem. 

d n v 
This theorem gives an expression for -^ , where y is the product 

of two functions of x t say u and v. 

It should be noted first that, as shown in Chapter XI, § 1, the 

operator -j- or D obeys the distributive • and index laws and is 

commutative in regard to constants. In these respects it is similar 
to A (see Chapter II, § 6). * 

tsj d , x du dv 

Now _ (wu) = „_ +tt _, 

and we may therefore write 

D (uv) = A (uv) + A (uv) 

= (A + A) w*>> 

where A operates only on u and differential coefficients of u and 
A operates only on v and differential coefficients of v. 
Therefore 

D n (uv)^(D 1 + D^) n uv 

= [A n +rc A n ~ 1 A + (2) A n ~ 2 A 2 + ... + A n ] ™. 

Now 



AV**-.^, 



dv\ dv d n ' 1 u 



A-A(»)-A-(»$)-|,^„ 

etc. 
Hence 

d n (uv)__ dPu dv_ d n ~ l u ( n \d?v d n ~ 2 u d n v 

(1). 



76 DIFFERENTIAL CALCULUS 

This formula can also be established by induction, but the proof 
by this method is left as an exercise to the student. 

4. Examples of the application of the formula are given below. 
Example 1. y = x*eP x . 

Since 7 n can be written at once as a"^, it is convenient to 
take x* as the factor v, and «■■ as the factor u. 

Then -j- = 3d? 3 , etc. 

ax 

Therefore 

y n = &&#* + w . 3^ . a 71 " 1 ^ + Q) ■ fe ■ a n " 2 ^ + (*) • 6 . a n - 3 e aa? . 

Example 2. If y = # n (log a?) 2 , 

prove that 

^11+2 + ^*1+1 = 2 xn!. 

We have y l = x n . — + (log #)* n . a;"" 1 

x 

or #2^ = 2# n log x + 7i# n (log xf 

= 2x n \ogx + ny (i). 

Differentiating both sides of the equation 

#2/2 + 2/i= \-2logx.nx n ~ 1 + ny 1 

x 

or x*y 2 + xy x = 2# n + 2nx n log a; + nxy x . 

Substituting for x n log x from equation (i) we obtain 

a?y 2 + #2^ = 2x n + nxy l — n 2 y + nxy l 

or a?y 2 — (2n-l)xy 1 + n 2 y = 2x n (ii). 

By differentiating each term of this equation n times, making 
use of Leibnitz's Theorem, we arrive at the result 

# 2 2/n+2+ ny n+i .2x+ n(n-l)y n 

-{2n-l)y n ^.x-n{2n~\)y n 

+n*y n 

^2/n+2 + #2/n+i = 2xnt. 



CHAPTEE XIII 

DIFFERENTIAL CALCULUS. EXPANSIONS. 
TAYLOR'S AND MACLAURIN'S THEOREMS 

1. It is often necessary to expand f(x) in a series of ascending 
powers of x. This has been done by ordinary algebraic or trigono- 
metrical methods in such cases as (x + a) n } eF, log e (1 + x), sin x, etc. 

The various methods which can be employed may be summarised 
as follows: 

I. By algebraic or trigonometrical methods. 

II. By the use of Taylor s or Maclaurin's Theorem. 

III. By the use of a differential equation. 

IV. By differentiating or integrating a known series. [For 
integration, see Chapters XVI-XX.] 

2. I. An example will make this method sufficiently clear. 

1 1 







t 2^ 12 720^' 
by actual division. 


** 


3. II. Taylor's and Maclaurin's Theorems. 


1 


Assume that f(x + h) can be expanded in j 
integral powers of h. 

Let f(x + h) = a + bh + ch? + dh?+ . 


a, series of positive 


Differentiating with regard to h, 

f(x + h) = b + 2ch + 3cZA a + 


... , 


f"(x + h)= 2c +3.2dh + , 
and so on. 
Put h = 0, then we have 

«-/(*). b=f(x), oJ^&, 


... , 
etc. 


Hence 


h n 



/(«+*)-/(«)+ A/(«)+i 1 /'(*)+j- f r w+... +^/ w («>+.... 

This result is known as Taylor's Theorem. 
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4. If in the above result we put a? = Owe obtain 

/(A)=/(0) + A/'(0) + | 2 ! /"(0)+...+^/(»)(0) + ... > 
which, by altering the notation, may be written 

/(*)-/<o)+«/(o) +|r (0)+,.. +jj/w (o)+ .-.. . 

This result is known as Stirling's or Maclaurin's Theorem. 

The student should verify the common algebraic and trigono- 
metrical expansions by means of the above theorems. 

Example 2. 

/(«)-4og(l+'n /(0)-log2, 

■ ^' ( *> = (TT*r /"(0) = i, 

/"'W-(TT^' /'"(0) = o, 

etc. 
Whence, by using Maclaurin's Theorem, 

log(l + e *) = /(a;) = log2 + |+J-^.... 

5. The above results are not universal in their application. Thus 

it can be shown that the method fails if/(#) or one of its derivatives 

becomes infinite or discontinuous between the specified range of 

values of x. Further, the series obtained must be a convergent series. 

Lagrange has shown that the remainder after the first n terms have 

h n 
been taken from Taylor's series can be expressed as —f^ n) (x + Oh), 

where is a positive proper fraction. The corresponding value of 

x n 
the remainder in the case of Maclaurin's Theorem is — / (n) (Ox), 

Unless therefore those expressions tend to vanish when n becomes 
infinite, the series is divergent and the method fails. 
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6. III. The use of a differential equation. 

The following are examples of the use of this method. 

Example 3. 

y = (l + x) n m a + a^x + a 2 a? -f a 3 x* + .... 
Then y x = n (1 + a) 71 " 1 or (1 -f x) y x = tiy. 

Differentiating the first equation, 

y 2 = a x + 2a 2 # + 3a 3 ^ + . . . . 
Therefore 

(1 +#)(a!+ 2a 2 # + 3a 3 # 2 + ...) = n(a + OiX + a^x 2 -^ ...). 
Whence, by comparing coefficients, 

O! = na , • 
2a> 2 + a l = waj, 
3a 3 + 2a 2 =na 2 ,- 
etc. 

Putting # = in the original equation we have 

a =l, 
and successively 

a x = wa = n, 

_ (n-2) „ j(nrl)(n-2) 
a 3 - — — a,- ^ , 

etc. 
Hence (1 + co) n = 1 + nx + (*) ^ + (3) & + ... . 

Example 4. 3/ = sin -1 # = a + a^ + a 2 x 2 + . . . , 
1 

and (1-^)^=1. 

Differentiating (1 - x 2 ) 2y x y 2 - txy* = 
or (l-x 2 )y 2 = xy 1 . 

Now sin -1 x = — sin -1 (— #). 

Therefore 

a + c^a? + OzX* + a,^ + ... = — (a,,— a^+'aatf 2 — 030 s + ...). 
Whence a = a$ = a 4 = . . . = 0. 
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Thus we may write 

y = ai x + a z a?+ ... +a» l+1 d aB+1 + ..., 

x 



y l -a 1 + 3a 3 # 2 +... + 2n+lo 2n+1 # 2n + ..,, 



2/ 2 = 3.2a 3 a;+...+2w + 1.2n.a 2n+1 aJ 2n - 1 +.... 

Therefore 



(l-a? 2 )(3.2.c^a?H-...+2n4-1.2n.a 2n+1 a^^ 1 4-..0=^K + 3a 3^+--- 



+ 2n + la2n+i a 27 * +...)• 
Whence a 1 =3.2a s , 

3a 3 = 5.4a 5 — 3.2os, 

(2n + 1) ^ = (2w + 3) (2ra + 2) a^ - (2?i + 1) 2n . 0^+, 
or (2n + l) 3 a„ +l - (2/i + 3) (2n + 2) a^. 

But in the limit when x -*■ 0, - = 1. 

Therefore «i — 1» 

and accordingly a 3 = ^— « , 



a B = 



1.3 



* 2.4.5' 
etc. 

, 1 ^1.3 a\ 

and sm i«-« + -.- + — •.- + ... . 

7. IV. Differentiation or integration of a known series. 

The method is sufficiently indicated by the following examples. 

Example 5. sin x = x — ;- + -r. — . . . , 

dsina? - a? X* 

___ ==C0SiC=1 __ + __.... 

We have therefore obtained the expansion of cos x in terms of x 
from that of sin x. 
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Example 6. Let 

y = l°g«( 1 - x) = a + a^ + a^x* + ..., 

^ = -= = a 1 + 2a 2 x + 3a z x 2 +.... 

ax L — x 

But = l+x + x 2 + x*+.... (#<1) 

J. ~~ it? 

Therefore, comparing coefficients, 

Oi = — 1, a 2 = — £, a 3 = — J, etc. 
Also putting x = in the original equation, we find a Q = 0. 

Hence log*(l -a)~- &- ^ - - - ... . (x<l) 

Thus our knowledge of the expansion of in terms of x enables 

X "~" x 

us to find the expansion of log e (1 — x). 



H. T.B.I. 



CHAPTER XIV 



DIFFERENTIAL CALCULUS. MISCELLANEOUS 
APPLICATIONS 

Limiting Values or Undetermined Forms 

1. In certain cases the value of a function cannot be at once 

x n — 1 
ascertained. A common form is illustrated by Lt =- , the value 

of which is not apparent since it takes the form ^ . The function 

in this case is said to take an Undetermined or Indeterminate Form, 
and since, in the limit, its value is represented by the ratio of two 
indefinitely small quantities, it may be expected that the methods 
of the Differential Calculus will enable the Limiting Value of the 
function to be found. 

2. The process to be employed can best be illustrated by examples. 

|~# n — 1~| 
Example 1. To find Lt . 

Let x = 1 + h, then 

J5i[Trr]"is[ ( a+*Hi] 

-Lt[«+(;)A+(i).*+...] 

= tt. 



Example 2. To find Lt H° g f * + x) ] . 
Lt Pog(l + «>1, Lt 



X 2 X s 

X "2 + S"" 

x* X s 
* + 2! + 3! + '" 



= Lt 1 - x + ^ x 2 - . . . (by actual division) 



= 1, 



since the series obtained by division is obviously convergent. 
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Example 3. To find Lt ^- — ^ . 






T^I> a;,08C + ^ (1O8C)2+ -"-l] 

= Lit 5 — = 

5T I J* log C + |j (log C) s + . . . J log 

-Lt 



4logo + g(logc)' + ...} i( -M 



T fffflogfflogc+termsinvolvingar'andhigherpowersof x 
xhgg + £(loggy+... 

T ff log g log c+ terms involving x and higher powers of a? 

log ^r + ljOog <?)*+... 
= flrlogc. 



3. The formulas of the Differential Calculus can be used for the 

solution of problems of this class, in the following way. 

fix) 

Let , \ \ take the form -^ when x-^a. 
<f>(x) 

Then if x = a + h, by Taylor's Theorem, 

/(*) = /(<*)+*/'(<*) + -•- 
<£(#) <j){a) + h<p'(a)+...' 

Therefore, when x-*~ a, and consequently A, /(a) and <f>(a) all 
-*-0, we have 

If this is still an undetermined form, the process can be repeated. 
Taking Example 3 above, we have 

f(x) = g*-g, <£(#) = #*- 1, 

/' (#) = g^^ogg log c, <£' (x) = tf logg. 

6—2 
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Therefore Lt fM = Lt^J°^ lo g C 

= (grlogc, as before. 

4. Other Undetermined Forms are found when a function in 
the limit assumes one of the following forms : 

Oxoo, — , oo-oo, 0°, oo°, or l 00 . 

00 

These can all be made to depend upon the form -r . 

Thus if /(#) and <f> (x) be the two functions involved, in the first 

f(x) 

case by writing f{x) x <j> (x) as ^^^ we have the form jr and can 

proceed as before. 

Similarly for the second form. 

The third form can be written as oo \l J- ; the second term 

within the bracket can then be evaluated, and the quantity within 
the bracket will then be found to be zero or otherwise. In the 
former event the product is of the form oo x and can be evaluated ; 
in the latter event the value of the product is clearly infinite. 

The last three cases are solved by taking logs, when the logarithm 
of the function will be found to fall within one of the cases already 
discussed. 

5. In connection with the form l 00 , it is useful to remember 
that, from the development of the Exponential Theorem in algebra, 
it follows that 

and Lt(l + -y = e». 

Maxima and Minima 

6. If /(#) continuously increases with x until it attains a certain 
value a and subsequently decreases, the value a is said to be a 
maximum value of the function. Conversely if the value of the 
function continuously decreases, while x increases, until it attains 
a value b and subsequently increases, the value b is said to be a 
minimum value. 
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By this definition there may be several "maximum" values and 
several "minimum" values; further, a "minimum" value may con- 
ceivably be greater than one of the "maximum" values. 

These conceptions are illustrated in the following diagram. 

Thus the points A and D are maxima, whilst the points B and 
- E are minima. At the point G 
the value of the function ceases 
momentarily to increase but is in- 
creasing both before and after the 
point is reached. Such points, 
% where the rate of increase or de- 
crease in the value of the function tends to become constant as the 
point is reached, are called points of inflexion. 

7. In obtaining criteria for ascertaining the maximum or 
minimum values the following considerations are of importance. 

If f(x) increases in value with x, it is clear from the definition 

that -j- is positive. Conversely if f(x) decreases with x t -^ is 

negative. 

Further, if f(x) is continuous and assumes two values f(a) and 
/(&), one of which is positive and the other negative, it follows that 
for some value of x between the values a and b,f{x) must be either 
zero or infinite. 

8. Now for a maximum value the function increases up to that 
value and then begins to decrease. Therefore at a maximum value 

-j- must change from positive to negative ; similarly at a minimum 

value -j- is changing from negative to positive. It follows from the 
preceding article that at a maximum or minimum point the 
value of -~ is zero or infinite. 

9. Further, at a maximum point -j- is decreasing from positive 
to negative, and therefore -^ must be negative or zero. At a 
minimum point the converse is the case and -r* is positive or zero. 
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Hence in order to find the maximum and minimum values of 
f(x) it is necessary 

(i) to find the values of x which satisfy the equation /' (x) = 
or/'(#) = oo, 

(ii) to ascertain for each of these values whether/" (x) is negative 
or positive. 

If /" (x) = 0, a further test is necessary. For the investigation 
of these cases, which include values which give rise to points of 
inflexion, the student is referred to more advanced treatises on the 
subject. 

Example 1. What fraction exceeds its pth power by the greatest 
number possible ? 

Let x be the fraction ; we have to find the maximum value of 

y = x — x p , 

dy 



dx 



= 1 -pxP~\ 



For a maximum or minimum value -^- = 0, therefore 





1- 


-p X P-i = 0, 


and 




x*- 1 = - 
P 


or 




i_ 

x=p p~K 


Now 


d*y 
dx* 


-p(p-l)x*>-* 




rs 


P-2 

-p(p-l)p p -\ 



and is therefore negative. 

Hence the above value of x is a maximum value. 

Example 2. What are the dimensions of the largest rectangular 
box on a square base the area of whose surface does not exceed 
12 square feet? 

Let x be the volume of the box, a the length of the base and b 
the height. 

Then the surface is 2a 3 + 4a6 = 12, whence 

7 6-a* 
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Also x — a'6 

2^ - da- 2' 

dx 3a a 

<fc* 2 ' 

d 2 « „ 

For a maximum or minimum value -y- = 0, therefore 

or a = V2, whence 6 = V2, 

and since -t-j is negative when a = V 2, this value gives a maximum 

value for x. 



CHAPTER XV 

RELATION OF DIFFERENTIAL CALCULUS TO 
FINITE DIFFERENCES 

1. By Taylor's Theorem we have 
But/(# + 1) = Ef{x\ therefore if the operation -y~ be denoted by 



D we may 


write 


symbolically 








Ef(x) = (l+D + D>+. 


••)/(«) 






-•»/<*) 




or 




E = e D , 




whence 




A = e D -l 




and 




D = log e (l + A). 




Therefore 




<Z A J A' 
dx 2 + 3 


... • 



(1). 

If the interval of differencing be A, we have 

/(■+*)-/(*)+fc/'(«)+|5r(c)+... 

and therefore A = e hD — 1. 

Further, D 5 = [log, (1 + A)] 2 , 

, d 2 y T A A a A 3 "I 3 

whence »"[ ¥ + «~'"J 8 ' 

= A»y-A»y + ttA<y-$A»y + (2). 

Similarly g = A«y-f A'y + f A'y - (3). 

2. The above formula (1) gives a convenient expression for the 
differential coefficient in terms of advancing differences. But from 
the nature of the differential coefficient, it may be expected that a 
better result will be obtained by applying the formulas of central 
differences. 
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Thus, taking Stirling's formula, we have 



2 " 2 



«(«*-!) A»/(-l) + A'/(- 2) 
+ 3! 2 



A/(0) + A/(-l) A'/(-l) + A'/(-2) 
/ W- g 12 + - 



Differentiating with regard to x, 

f± 

+ terms involving x and its powers. 
Putting x = we obtain 

//(0) ,A/(0) +2 A/(-l)_Ay(-lHA>/(-2) + (4) 
Taking the first term only, we have 

f > (0) _ a/(0) + Afl- 1) = /d) -/(-i) (5)> 

a very useful approximation, which, by altering the origin and the 
unit of measurement, may be expressed more generally as 

f > ix)= x*+*)-f(«-Q ( 6). 

Alb 

If h = \ we arrive at the approximate result 

/'<«)« A/(«-|) (7) 

-A(l + A)-*/<*) 

= [A -JA« + §A»-AA« + dfeA» -...]/(*)• 
Comparing this with the value given by formula (1) we see 
that the error involved in taking f {x) = A/(# — \) is 

[- A A-» + A ^ - ^V A 5 + • • •]/(*)> 
which is approximately equal to — ^A 3 f(x — §). 
Hence we may write 

A») = A/(«-i) -£*/(*-» (8)- 

3. The introduction of the second term in equation (4) gives 
the useful approximation 

f , (Q) 8 {/(l) -/(- 1)} - 1/(2) -/(- 2)} (9) 

The other formulas of central differences, when treated as above, 
give these and other expressions for f (x). The above, however, 
are the most accurate and useful for general use. 
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4. Using the central difference formula, we shall obtain for the 
differential coefficients of higher orders, 

/"(0) = A»/(-l)-^A'/(-2) (10) 

/'"(Q)- A ' /( ~ 1)+A ' /( ~ 2) 



and 



= A 3 /(-f)approx (11). 



5. An example of the working of the above formulas is now given. 
Example. Find the first three differential coefficients of log a?, 
when x = 2. having given the values of log* 1'80, log* 1*90, .... 



Number 


Natural 
log 


A 


A* 


A 3 


A 4 


A 5 


1-80 


•587787 


•054067 


- -002774 


•000271 


-•000038 


•000007 


1-90 


•641854 


•051293 


- -002503 


•000233 


-•000031 


•000005 


2-00 


•693147 


•048790 


- -002270 


•000202 


- -000026 


•000004 


2*10 


•741937 


•046520 


- -002068 


•000176 


- -000022 


•000004 


2-20 


•788457 


• -044452 


- -001892 


•000154 


- -000018 




2-30 


•832909 


•042560 


- -001738 


•000136 






2*40 


•875469 


•040822 


- -001602 








2-50 


•916291 


•039220 










2-60 


•955511 













Using the advancing difference formula (1), starting at the term 
2*00, and remembering that the interval of differencing is not unity 
but '1, we get 

f « = 10 (-048790 + ^ + « + ^26 + -000004) 

= -50000. 

Similarly 

/" (x) = 100 (- -002270 - -000202 - {£ x -000026 - f x -000004) 

= - -2499, 

/'" (x) = 1000 (-000202 + f x -000026 + J x -000004) 

= -248. 

The true values are, of course, given by the differential coefficients 

1 12 
of log a;, which are respectively -, — — , — . When x = 2 these 

X Xr X* 

become £, —J, \\ the closeness of the above approximations is 
apparent. 

Formulas (6) and (9) give respectively for/' (x) values of '50042 
and -50000. 



OHAPTEE XVI 

INTEGRAL CALCULUS. DEFINITIONS 
AND ILLUSTRATIONS 

1. In Finite Differences we discussed under the heading of Finite 
Integration the problem of finding the value of /(0) -f/(l) + ... 
+f{n— 1), or, changing the origin and the unit of measurement, 
the more general series 

h [/(^) +f(a + h) +f(a + 2h)+... +/(a + n^Th)]. 
The Integral Calculus deals with the value which this summa- 
tion assumes when h becomes indefinitely small. 

2. Now let F{x) be a function such that f(x) is its differential 
coefficient. Then by definition 

where a x is a quantity that vanishes when A-^0. 
Then we have 

hf(a) =F(a + h) -F(a) +ha lt 

hf(a + h) = F(a + 2h)-F(a + h) +ha 2 , 

hf(a + 2h) = F(a + 3h) -F(a + 2h) + hct z , 



hf(a + n-lh) = F(a + nh)--F(a + n-lh) + hct n . 
By addition, 



A lf(fl) +f( a + A) +/(a + 2fc) + . . . +/(a + n - 1A)] 

= F(a + nh)-F(a) + h[a 1 + a 2 +...+a n l 

Now if a denote the greatest of the quantities a u c^, ... a n > the 
last term is clearly less than nha, or (b — a) a, if a + tiA be put equal 
to 6. This term therefore vanishes in the limit and we are left with 
the relation 

Uh[f(a)+f(a + h)+...+f(b-h)] = F(b)-F(a). 
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This result is usually expressed in the form 
rb f(.x)dx = F(b)-F(a), 



I 

J a 



where F(x) is a quantity such that/(#) is its differential coefficient 
and dx represents, in the limit, the indefinitely small interval 
between the terms which are summed. 

The expression I f(x) dx is called the definite integral of f(x) 

J a 

with regard to x, and b and a are called respectively the superior 
and inferior limits of integration. 

Where no limits to the summation are expressed and we are 
concerned merely with the form of the function, we obtain the in- 
definite integral. 

3. A geometrical illustration of the process of integration may 
be obtained as follows. 

Let the equation y=f(x) be represented by the curve in the 
diagram. Let OA = a and OB = b t so that 
MA**f(a) and NB=f(b). 

Further, let AB be divided into n parts 
each equal to A. Then the sum of the rect- 

angles shown is clearly equal to 2 hf(x) f and, 



x=a+h 




m 



A 



'« 



B 



in the limit, when h-^0 the value of the 

integral becomes equal to the area between the curve and the axis 

of x bounded by the ordinates MA and NB. 



CHAPTER XVII 
INTEGRAL CALCULUS. STANDARD FORMS 

1. In the preceding chapter we found that the process of integra- 
tion is the converse of that of differentiation. In other words, that 
given /'(#) we have to find a function f(x) such that f'(x) is its 
first differential coefficient. The analogy with the process of finite 
integration is apparent and the remarks made in Chapter VII, 
§ 5, apply equally to the present case. 

It may, therefore, be said again that the process of integration 
cannot be carried out for every function, that definite rules cannot 
be laid down to apply in every case, and that the student must look 
to applying the results of the differential calculus as a guide in 
solving any problem presented by the function under consideration. 

2. In the first place it is necessary to point out that the ordinary 
algebraic laws apply to the integrating symbol fdx in the same way 
as they have been shown to apply to the symbolic operations A and 

£. Thus: 
ax 

(1) The operation is distributive for, if u, v, w, ... be any func- 
tions of x, 



dx 
and therefore 



j- \ \udx+ lvdx + jwdx + ...> = u + v + w ... 



wd#+ lvdx+ \wdx + ...= l(u + v + w + ...)dx ...(1). 

(2) The operation is commutative with regard to constants, for, 
.„ dv d(cv) __cdv __ 

dx~ * dx dx 



therefore 



lcudx = cv = c \udx (2). 



3. It should be added that the process of integration introduces 

dv 
a constant into every indefinite integral, for, if u = -p , 

d (v + c) dv 



dx dx 



= u, 
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and therefore ludx = v + c, 

which result is, indeed, obvious from the consideration that a con- 
stant term disappears on differentiation. 

4. In accordance with the symbolic notation which has been 
previously developed (Chapter XV, § 1), we may write 

//<«>*-(£)"/<•> < 3 > 

-[ui^+55>« < 4 >- 

5. Standard Forms. 

Jm-i 



d( — ) 
U+l/ _ (n + l)d» _ 



(a) *\ Since yi/ a yi»Ti^ = 

f x n+1 

therefore \x n dx = + c. 

J n + 1 

(6 ) ^ Since ^ -1, 

dx x 

therefore J — = log a; + c. 

J & 

(c) e* Since ^ = e* 

eta; 

therefore I e^da: = e* + c. 



(d) «-. S i nce _Uog^ = a^a = 

therefore I a^da? = ; h c. 

J log„a 

Other integrals maybe obtained from a consideration of the table 
of standard forms of differential coefficients given in Chapter XI, 
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§ 3. The results are embodied in the following table. Their verifi- 
cation is left as an exercise to the student. 



Function 


Integral 


lx n dx 


n + 1 


\ -dx 

J X 


log e x 


J a x dx 


a x 
log e a 


i4*dx 


<* 


J cosxdx 


sin x 


1 sinxdx 


-cos X 


1 sec 2 xdx 


tsmx 


1 cosec 2 x dx 


— cot# 



In each case the constant c has been omitted for the sake of 
simplicity, but its importance in the result must not be forgotten. 

6. Sums of the above functions can be integrated by the use of 
the distributive property of the operation of integration. Thus 



\{x + 2 X ) dx = J J dx + ]2 X dx 



I 2 X 
log* 2 



+ c. 



Similarly Uax 4, + bx" 4 ) dx = a \x*dx + b jx"* dx 



ax 1 bx" 3 



+ c. 



7. The problem of integration may be presented in another way. 
Thus to find the sum of 



1 _J_ 1 

x x + m x + 2m 



+ ...+- 



x + xm 



when x is indefinitely increased. 
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The series may be written as 

1 r 1 11 



m 



i+2 i+^ 

x x 



+ ... + • 



m 
l + #.— 
x 



1 m 



The terms in the denominator consist of values of a quantity, which 



m 



can be represented by y, increasing by equal increments of — . Also 

x 

the sum of the series is multiplied by — , which is the increment 

x 

in the value of the denominator. 

Since the initial value of the denominator y is unity and the 
final value is 1 -f m, from our definition of an integral we may write, 



when x-^qo and therefore — ^0, 

x 






- + —J— + 1 

X X 



1 -I m _ 1 f 1+m dy 

, t __ m x ~~ mJi y 



l+#.- 
x 



= — losr 1 + m. 



CHAPTER XVIII 

INTEGRAL CALCULUS. METHODS OF INTEGRATION 

I. The Method of Substitution 

1. Various methods are available for treating the integration of 
functions which are not among the standard forms. Of these methods, 
one of the most important and most easily applied is the method of 
substitution. By this method the independent variable x is changed 
to y where y is a function of x. Thus an integral in terms of y can be 
obtained which can frequently be made to assume a standard form. 

For suppose y = F(x). 

Then dy = F'(x)dx. 

So that if the original integral is of the form lf[F(x)] F'(x)dx, 

it can be expressed as lf(y)dy, which may be a standard form. 

No general rules can be given, but the following examples will 
indicate methods which can frequently be employed with advantage. 

2. (1) j(a + bx) n dx. 

Put a + bx = y. 

Then bdx = dy. 

Therefore l(a + bx) n dx = f^dy = , *** . +c = ( a +^T + c. 
J J b * (n + l)b (n + l)b 

f a?dx 
KZ) }(a + bx) n ' 

Put a + bx = y. 

Therefore 



f a?dx _ f (y - a) 2 dy 



6 3 y n 
_1[" 1_ 2a 



(n-2)j^ T (»-l)y 
where y = a + bx. 

H. T. B. I. 



+ c 
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(3) j(b + 2cx) (a + bx + c# a ) n dx. 

Put a + bx + czP^y. 

Then (b + 2cx) dx = dy. 

Therefore 

I (b + 2c#) (a + bx + c^) n cfcc = \y n dy = -^ — - + A 



(a + bx + c# 2 ) n+1 , 
= — h fc> 

71+1 



(4) f (log^) n <fe > 

. J # 

Put log a: s= y. 

Then -dx = dy. 



Therefore 



r(iog£)^ = /- j^ (jogxy^ 

J ar .F y n+1 n + 1 ' 

(5) Jtc^t- 



Put /(*)=y- 

Then f (x)dx = dy. 

Therefore 



(6) f^fL. 

Put &=y- 

Then 2xdx = dy. 

-*[Ml%] 

=i[y-iogi+y] + c 



a) 

Put 

Then 
Therefore 

sin xdx 
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sin xdx 
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f sin 

JaTi 



■bcosx* 
a + bcosx = y. 
— b sin xdx = dy. 



J a -f 6 cos # 
(8) 



1 fdy 1, . log(a + 6cosa;) 



dx 
sin# 



dx 



2 sin -= cos ~ 



sec 2 ^ a# 



Put 

Then 

Therefore 

Corollary. Since 



2 tan^ 
tan| = y. 



| sec 2 ^ dx = dy. 

J sin# J y oa 

x 
= log tan h + c. 

cos# = sin(-^ + x\, 

/^ =1 °s tan S + i) 



+ c 



II. Rationalisation by Trigonometrical Transformation 

3. The above methods can be extended by the use of trigono- 
metrical functions and the relations which exist between them. 
The following are examples : 



(1) 



Let 

Then 

and 



f dx 
J Va 2 -^ 2 ' 

x = a sin y. 
dx — acosydy t 
Va 2 — x*= a cos y. 



7-2 
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Therefore . = \dy = y + c 



x 
= sin" 1 - + c. 



v } J *Jx*-a? 

Let x = a sec y. 

Then dx= a sec y tan 3/ dy, 



and Va 3 — a 2 = a tan y. 

Therefore 

rd^ = rjy_ = logtan(j + fW c 
JV^-a 3 J cosy 5 V4 2; 

l+tan| 

-log - + c 

l-tan| 

. 1 + sin y 
cosy 



1 + 



y-i 



= log +c 

6 a 



x 



'x + Vtf 2 - a 2 \ 

+ c. 



(3) 



=Io n — — j 

[ __dx__ 
J vV + a 2 ' 



Let # = a tan y. 

Then <fcc = a sec 3 y dy, 



and V^Ta 2 = a sec y. 

Therefore 

r^^ = fA = logtan(j4-|)-fc 
JV^Ta 3 J cosy & \4 2; 

= log (sec y + tan y) + c [See Ex. (2)] 
; + *Jx* + a 2 \ 



=H (^f^y c . 
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From the above integral, two other\ itaj&r'fanV, Jfatfris' £&n be 
obtained. 

dx 



(4 ) J d * 



This can be written as 



' dx f d (ex + b) 

/ (aK + 6) 2 + ac-6 a ~J Vc {(ex + 6) 3 + ac - b 2 } ' 

Disregarding for the moment the constant multiplier, this can 
clearly be transformed into one of the forms given in Examples 
(1)— (3) above. 

Corollary. Since 

\ =21og(\^^ + \/^6) + c, 

J n(x — a)(x — b) 



f dx . Ix — a , 
I = 2 sin * a / f h < 



< s > /^ 



(p + £#)<&£ 

+ Ibx -t- c# 2 



( n.«, )+ (g=gE] 



Va + 26# + c# 2 



da? 



_ q f d(a + 2bx + cx*) (pc-q b) f cfo 

~ 2c J Va + %bx + c# 3 c J »Ja + 2bx + cx 2 ' 

The value of the first integral is 2 (a + Ibx + as 2 )* and the second 
integral can be dealt with as in Example (4). 

4. General remarks in regard to the method of substitution. 
When the precise substitution which will enable the integral to 
be solved is not apparent the device illustrated in the following 
examples is often of assistance. 

f (x + a)dx 
K } J(x + b)\/x+c' 

The appropriate substitution is not apparent, so assume 

V# -f c = y T y 

where r may have any value, to be determined subsequently. 
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and dx = 2ry* r - 1 dy. 

The integral thus becomes 

f (y ar + a-o)2ry ar ~ 1 rfy 

Clearly this integral can be evaluated if ^ = y 2r—1 , i.e. if r= 1, 
when we get 

J Vc - 6 Ly - ^c - 6 y + Vc-6j ^ 
a a — 6 . y-\/c-b . 



cfa? 



where y = V# + a • 

(2) f 

J(x-p) v (a + 26a; + as 2 ) 

Let us take the linear function as likely to lead to simpler results 
and, as before, put 

x — p = y r 

so that cfa = ry r - 1 dy. 

Then the transformed integral is 

f ri/ r ~ 1 dy f rdy __ 

Jy r 's/a + 2b(y r +p) + c (y r + p) 2 Jy Va + 26 (y r + jp) + c (y r +p) 2 ' 

Obviously if r m — 1, the denominator reduces to the form 
Va' + 26'y + c'y* 
and the expression can be immediately integrated. 

(3) f-Ar- 

This example is a little more difficult than the two previous 
ones. As before, let x — y r , so that dx = ry*""" 1 dy. Then we have 



f rfa? _ f ry r ^ 1 dy 
J(a + cx*f J{a + cy^f 
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The expression within the brackets in the denominator is of an 
even degree in y, so we will express the integral in the form 

f ry^djtf) f ry*-*d(tf) 

h(a + cy*f J 2y* r {ay~^ + cf ' 

This can clearly be integrated if r — 2 = 3r, when r = — 1. The 
integral then becomes 

1 f df = 1 f d(ay* + c) _ 1 |/: 

2J( a2/2 + c )5 2a J (ay* + cf a(atf + c)* 

+ fc 



a(a + ca?f 

5. Definite Integrals. The above examples of the application of 
the method of substitution all relate to indefinite integrals and the 
final form of the indefinite integral is expressed, in each case, in 
terms of x. 

In the case of a definite integral, a step may be saved in the 
operation by avoiding the final substitution of x in terms of y, 
provided appropriate changes are made in the values of the limits 
of integration. An example will make the point clear. 



Example. \ 

Jo 



xdx 



'oVl +# 2 ' 
Put x 2 = y. 

Then 2xdx=dy. 

Also when x = 0, y = and when x = 2, y = 4. 
Thus we get 

r2 xdx r* dy r . "I* .— 

III. Integration by Parts 
6. From the Differential Calculus we have 

d (uv) __ dv du 
dx dx dx' 

Whence, by integration, 

uv = ju -T- dx + p-T- dx, 
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f dv , f du 7 

or 1^-7- dx — uv — \v-T- dx. 

J dx J dx 

Consequently an integral of a function of the form u -7- can 

OLX 

always be made to depend upon that of a function of the form v -7-. 

The latter may prove to be a standard form and the use of the 
formula given above may enable an apparently intractable ex- 
pression to be integrated at once. 

The advantage of the method is best exhibited by applying it to 
a few elementary cases. 

(1) Jxlogxdx=^jlogx-^ dx 

-s[*"*-H*FH 

= 2# 2 log# — ^ jxdx 

(2) j x ^dx=±fx^dx 

= -\xe? x - \^ x -^dx\ 

1 T #*~\ 
= - \xe ax +c. 

(3) I sin" 1 xdx=\ sin" 1 x -7- dx 

. , f dfsin^x) 7 
= # sin" 1 x — \x , dx 

. , f xdx 
= x sin" 1 x — f , 



= a; sin" 1 x + - )— ^-/ 



= #sin" 1 # + Vl — a^-hc. 
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dx 



(4) [Va 2 + x 2 dx = jV^ + a? 2 ^^ 

= #Va 2 + # 2 -L- 

r ^ 

= Wa 2 + # 2 - -7=; 



ds/a? + a? , 



a?dx 



+ a? 



J L va 2 + a 2 J 

r r a 2 cfo 

= #Va 2 + # 2 - \>s/a 2 + a?dx+ . 

J JNtf + x* 

Therefore 2 / ^tf + x 2 cfo = # Va 2 + a; 2 + a 2 log [# + Va 2 + # 2 ] + c, 
and I'SaFVaPdx = | Va 2 + # 2 + |- log [^+Va 2 + ^] + 1 . 

(5) /log (x -f V# 2 + a 2 ) cfo? = flog (x + V# 2 + a?)-^-dx 

= x\og(x + ^¥T^)-jx dlo ^ X ^^^dx 

= a; log (a; + \/^ + a 2 ) - V^ + a 2 + c. 

7. The above process can be expressed in a general form as 
follows: 



ju t v t dt= U 



Ut —dt~ dt 



= u t jv t dt - I ~ (j v t dt) dt. 
To apply this formula to the definite integral I u t v t dt t we can 

J a 

write \v t dt as I Vjbdi, the upper limit being taken as t since the 



integral must itself be a function of t 
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Thus we get 

/ u t v t dt = \ u A v k dk — J -3-' ( I v k dk) dt 

= u b \ v k dk — u a \ v k dk—j -rri] v k dk\dt 

= u b I v k dk — -=-' ( I v k dk) dt, 

since J v k dk = 0. 
J a 

Alternatively, if we express \v t dt as the definite integral 
— I v k dk, we obtain 

u t v t dt — u b \ v k dk+\ -rpij v k dk)dt*. 

8. Reduction Formulas. Where integration by parts is not im- 
mediately successful, continuation of the operation may ultimately 
lead to the evaluation of the integral. 

Thus L*(?dx = fa? % dx = ^-\^% dx 

= x 2 0*-2 \<?xdx 

= a*e*-2Jx^dx 

= x i e*-2(xe e -je?dx\ 
= x 2 e?-2xeP + 2<? + c. 

9. The above case is an example of a Reduction Formula, this 
designation being used since the application of the method of 
Integration by Parts effects the reduction by successive steps of 
one term in the integral. 

Thus 

x m (l-x) n dx=- (1 **' + — -n(l-xy-*dx. 

J v f m+1 Jm+l v 

* See J.LA. Vol. 44, pp. 402-409. 



METHODS OF INTEGRATION 107 

Thus the degree of the term (1 — x) n has been reduced by unity. 
Successive applications of the formula will reduce the degree 
further, until, if n is an integer, the term ultimately equals (1 — x)° 
or unity. 

In the special case where the value of this integral is taken 
between the limits and 1, we have 

x m (l-x) n dx=\ -(l-#) n + = n (l-x^dx 

z m+1 (l-#) n ~ 1 *^ 



n f 1 <n 
= 7 x 

m + ljo 



whence by successive applications of the formula we find ulti- 
mately if n be an integer that 

I x m (l-x) n dx = r x m + n dx = — . 

Jo ri + ml J o n + m+ 11 

Many other important integrals can be dealt with in a similar 
way. 

IV. Integration by the use of Partial Fractions 

10. Any expression of the form ^ ^ , where F(x) and/ (x) are 

both rational integral algebraic functions of x, can be expressed as 
the sum of a number of terms of which the general forms are a r x r 

T 

and - — T -— . For if the degree of F(x) is equal to or greater than 
yx — c r ) 

that of f(x), by division we can obtain a quotient, together with a 

new fraction in which the numerator is of lower degree than the 

denominator. The quotient provides the terms of the first form and 

the new fraction can be split up by partial fractions into terms of 

the second form. The integrals of these general forms are known 

and consequently the whole expression can be integrated. 

The following are examples : 

(i) f (p+g;)^ ' i \\«j±p Ji±E\ dx 

' J(x-a)(x — b) a — b)\_x—a x — b] 

= ^T-b ^- aq +P ^ log (^ ~ a ) ~ ( b Q + P) lo S ( x ~ W + c - 



(2) / 



dx 1 



ax 2 -+- 2bx + c a 



dx 



ac-6 3 



/ 6\ ac- 
\ a) a 2 
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The form of the integral depends upon whether ac — b 2 is positive 
or negative. 

In the former case the integral clearly is equal to 

1 , , ax + b 
-r== tan -1 _ . 

v ac - 6 2 V ac - 6 2 

In the latter case the integral becomes 

1 . ax + b- V6 2 - ac 
-log- 



2 V& 2 - ac ax + b + */b 2 -ac' 

(3) / \«+f)t+2)-/oT*" 8 "^ + ^] <fa 

= 1^ — 3a? — log(a? + 1) + 8 log(# + 2) l 1 
= 81og3-91og2-|. 



CHAPTER XIX 

INTEGRAL CALCULUS. DEFINITE INTEGRALS. 
MISCELLANEOUS APPLICATIONS 

General Propositions 

1. It is desirable to place on record several general propositions, 
in regard to change of limits, which are in the nature of being 
self-evident. 

In Chapter XVI it has been shown that 

f b f(x)dx = F(b)^F(a) > 

J a 

where f{x) is the differential coefficient of F(x). 
It follows that 

I. \ b f{x)dx=\ b f(z)dz (1), 

since neither x nor z occurs in the result. 



II. 



f b f(x)dx=-f*f(x)dx (2). 

J a J b 



Thus the interchange of the limits results in a change of sign of 
the definite integral. 

On the left-hand side we have regarded the increment of dx as 
positive, so that, while x increases from a to b t the value of the 
integral is F(b) - F(a). 

On the right-hand side the increment dx is negative and x 
decreases from b to a, giving a value for the integral of 

F(a)-F(b). 
III. \ C f{x)dx=\ C f(x)dx+\ b f{x)dx (3). 

J a J b J a 

For the left-hand side is F(c) — F(a) and the right-hand side is 
F(c)-F(b) + F(b)-F(a). 

IV - f/(«) dx =/ /(<*-<") dx (4). 
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For if we substitute a — z for x, then when x = a, z = 0, and when 
x = 0, z = a. 

Also — dz = dx. 

Therefore 

f /(a?) d# = l°-f(a -z)dz= f /(a - *) ds (by II) 

Jo J a JO 

= \ a f(a-x)dx (by I). 

Jo 

2. Differentiation of Definite Integrals. 
Let w = I /(a?, c) (£a? 

Ja 

be a definite integral where the quantity c is independent of x> 

and the limits a and b are independent of c. 

du 
To find -j- let At* be the change in u corresponding to a change 

Ac in c. Then, since the limits are unaltered, 

fb 

Au = I {f(x, c + Ac) —fix, c)] dx. 

J a 

Therefore Ay = fM/(,, c + Ac)-/(,, c)}^ 
Ac J « Ac 

Proceeding to the limit, we have 

s-W*- < 6 >- 

Thus the differential of the definite integral is reached by a 
process of differentiating under the sign of integration. 

3. An important use of this theorem is that of finding the values 
of other integrals from those of known form. 

For example, if the equation 



fee I 

e~ a *dx = - 
Jo a 



* The student is referred to more advanced treatises for exceptions to this 
general result, 
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be differentiated n times with respect to a, we get 
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Jo 



x n q-mc fa - 



nn+l * 



/. 




4. Areas of Curves. 

It has been shown that geometrically the definite integral 

f(x) dx represents the area enclosed between the curve y =f(x) 
and the axis of x bounded by the ordinates x = a and x =■ 6. 

Example 1. Prove that the area of the parabola 
j/ 2 = 4<ax bounded by the curve, the axis of x and 
any ordinate is two-thirds of the rectangle con- 
tained by the ordinate and the intercept on the 
axis of x. 

Let 0^=6. 

Then area 0PN= j ydx = j *J~iaxdx = 2a* U x$\ 

= %a*b% 

«f(2a*6*)6 

=§PN.0N 
= f rectangle. 

5. Mean Value and Probability. 

Definite Integrals can be used to find the mean value of a 
function whose value is changing continuously by indefinitely 
small increments. 

Thus to find the mean value of f(x) for all values of x from a to 6. 
If we divide b — a into n portions each equal to A, the mean value 
of the functions /(a), f(a + h), . 

/(&)+/(<* + &) + 



•/(a + n- \h) is 



+ /(a + ft - lh) 



^ h[f(a)+f(a + h) + ...+f(a + n-lh)] 
6 — a 
since nh = b — a. 

If now we make h indefinitely small, we shall have the mean 
value of all values of f(x) from f(a) to f(b). Consequently the 



mean value is 



1 f b 



dx. 
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Example 2. A number n is divided at random into two parts; find 
the mean value of their product. 

Let x be one part, then the product is x (n — x). We have to 
find the mean value of this product for all values of x from to n. 
Thus we get for the mean value 

If* v , 1 /n 3 n 3 \ n a 

- x(n-x)dx = - hr-o" )=-*' • 
nJo * n, \2 3/ 6 

6. The class of problem in Probability that admits of being 
treated by the methods of the Integral Calculus can best be 
illustrated by an example. 

Example 3. Three events A, B, G are known to have happened 
in the same century. What is the chance that the events happened 
in the order A, B, G\ B happening within n years of the middle of 
the century? 

The chance that the events happened in the order A, B, G is 
clearly £, if there is no limitation as to when B can happen, since 
there are six possible orders in which the events A, B, C can 
happen. 

To ascertain the further chance that, with the events occurring 
in that order, B happened within n years of the middle of the century, 
let x be the number of years from the beginning of the century 
to the event jB, so that B occurs during the interval of time between 
x and x + dx. Then the A cases may have occurred in any of these 
x years, and the C cases may have occurred in any of the following 
100 — a? years. The total number of possible cases is therefore 
x (100 — x) in respect of the above interval of time dx. 

But B must have occurred within n years of the middle of the 
century, therefore x can have any value between 50 — n and 50 + n. 
Hence the required chance, being equal to the number of favour- 
able cases divided by the whole number, is 

JO 

and the answer to the question is 

6L 3 \^/ 4 vWJ"2llOoJ"3\IOoJ ' 
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7. An example of a somewhat different type is given in the 
following problem. 

Example 4. An event has always happened on an average once 
a year. Find the chance that it did not happen in a given year and 
prove it equals e" 1 . 

The event has happened on an average once a year. The chance 

that it occurred in any given nth part of a year is, therefore, -, and 

n 

the chance that it did not occur in that period is 1 — . 

By the ordinary rules of probability, the chance that it failed to 
occur in n consecutive periods each of -th part of a year is 

(l\ n 
1 — j . To obtain the solution to the problem, we proceed to the 

(l\ n 
1 — ) , which by algebra equals e~\ 



H. T.B.I. 



CHAPTER XX 

APPROXIMATE INTEGRATION 

1. In many cases integration, or continuous summation of the 
values of a function, cannot be accomplished, either because the 
quantity to be integrated cannot be expressed as a mathematical 
function, or because the function itself is not capable of being 
integrated directly. 

In these cases formulas of approximation can be used, which may 
conveniently be divided into two classes, viz. formulas expressing 
the value of the definite integral in terms of 

(i) the sum of the successive values of the function and of its 
derivatives, 

or (ii) the values of isolated values of the function, not necessarily 
successive. 

We shall now proceed to consider formulas of the first of these 
classes. 

2. The Euler-Maclaurin Expansion. 
Let 2f(x) = F(x), 

so that f(x) = &F(x). 

Then "% f(x)=f(0)+f(l) + ...+f(n-l) 

- F(n) - F (0). [See Chapter VII, § 4] 
Now F(x) = A- 1 f(x) 

= ( e D-l)-*f(x), since A = 6^-1, 
which, by actual division, 

(1 ID D' \, /N 

- D J{X} 2 /W + 12 dx 720 d* 3 + - 
- J/(c) dx - */<*) + hf (*) - 7k/"' (*) + -• 

Therefore 

F(n)-F(0)=jy(a ! )da ; -i{f(n)-f(0)}+^{f(n)-f(0)} 

-7h{f'"(n)-f'"(0)} + .... 
But F(n)-F(0)=/(0)+f(l)+...+f(n-l). 
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Hence 

/ o 7(*)«te-i/(0) +/(!)+•..+/(« - 1)+ */(n)- A {/' (n)-/'(0)} 

+ Tk{/'"(«)-/" / (0)}-... (1). 
This result is not limited to the case where the ordinates are at 
unit distance apart, for, as has been remarked in connection with 
Finite Differences, by changing the origin and unit of measurement 
the formula can be given a more general form. Thus 

± /(*) dx = if (a) +/(a + r) + . . . +/(o + » - lr) + J/(a + nr) 

r J a t 

- fr {/' (a + nr) -/' (a)} + ^ {/'" (a + «r) -/'" (a)} 

-3^0^ <V)(a + Wr)_/(V)(a)} + < 2 >" 

f 25 cfcc 
Example. Calculate — . 
r J20 # 

Taking values at unit intervals we have, remembering that 

dx ■ x 2 dx 3 x 4 ' 

J 20 a: 2"20^21^22^23^24^2 25 12\ 25» + 20V 

L20'v 25 4 + 20V 



+ 120' 



025000 
047619 
045455 
043478 
041667 
020000 



•223219 - ^ (-0009) + ^ (-00000369) 
= -223144. 
The value of the integral is clearly equal to log e f§, which to 
five places of decimals is -22314. 

3. Woolhouse's Formula. 

This formula gives a relationship between the sum of consecutive 
ordinates and the sum of equidistant ordinates at greater intervals. 
It can be derived from the Euler-Maclaurin formula as follows. 

8—2 
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Formula (1) gives 
f/(x)dx~if(P) + f(l) + ... + f(n-l) + lf(n)-&lf(n)-f(fi)} 

+tU/'»-/"(<>)}-.... 

If the interval is — the formula becomes 
m 

«/>)A»-V(0)+/£)+/© + .-.. 

+ i/(«)- I ^{/'(«)-/'(0)} + 72L,{/"(n)-/'"(0)}-..., 
whence from the first of the above equations 
m f "/(*) dx - m {/(0) +/(1) + . . . +/(»)} 

JO 

-f{/(0)+/(«)}-5{/(^)-/'(0)}+^{/ / "(n)-/'"(0)}-... 

and from the second equation 

m jj(x) dx = |/(0) + / (4) + • • • +/(»)} " * {/(°) + /(»» 

whence, by subtraction, 

/(0) +/© + - +/(») = * {/(°) +/(!) + • • • +/(«)} 

+m£\r(*)-r(o)}---- (3) - 

4. By putting a = and r= 1 in formula (2) and proceeding as 
above, we obtain the following result: 

/(0) +/(1) +/(2) + ... +/(mn) « » {/(0) +/(«) +/( 2w ) + • •• 

+ f(r n n)\-^{AO)+f(mn)}- , ^{f'(mn)-f(0)} 

+^tr(«»)-/ w (o)}- w. 

This formula enables the sum of consecutive terms of a series 
to be expressed in terms of those at greater intervals. Thus the 
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work may often be shortened, particularly where the values of in- 
dividual terms of the series require to be obtained by calculation 
in the first instance. 

5. Lubbock's Formula. 

Lubbock's formula is similar to that of Wool house, but it sub- 
stitutes finite differences for the differential coefficients. 

Now Z> = A-£A 2 + £A 8 -;JA 4 +... 

and i) 3 = A 8 -fA 4 +.... 

Therefore, substituting in the above formula (3), we get 

'2> 



/<°> + /® + /(S+ •"+.«»> 



ra — 1 



-*> l/(0)+/(i)+/(S)+ ... +/(»)} -:-_-t(/(o) +/(,)( 

+ ( -=^S^(AyW-A'/(0)! + ... (5). 

6. By an alteration in the unit of measurement the formula may 
be expressed as 

/(0)+/(l)+/(2)+ . . . +/(i«n)-»{/(0)+/(»)+/(2«)+ . • • +/M| 
- ^ {/(0) +/(««)} - ^=i lA/(«») - A/(0)} 

+ ^T^"/(^)-^/(0)} 

+ (W2 " 1 480 9 J" 1) ^(^)-Ay(0)} + ... (6), 

wfore A, A 2 , ... express the values of the differences taken over the 
interval n, and not those taken over unit intervals. 

In this form the formula corresponds to formula (4) given above. 
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7. In cases where the terms of a series decrease steadily and the 
term corresponding to the upper limit tends to the value zero, the 
above formulas can be simplified somewhat, since the final term of 
the series, its differences and derived functions all vanish. Also it 
is frequently the case that the end terms of the series are unim- 
portant and thus m and n may be taken at any suitable figures so 
as to correspond with tabulated values of the function. 

To illustrate this we will calculate, first by Woolhouse's formula 

£C = 00 

and then by Lubbock's, the value of 2 (1*1)"*. 

a?=lO 

Example 1. By Woolhouse's formula. 

Let us use formula (4) and take n = 15. The values of (l'l)"* 
tend to become unimportant when # = 100, therefore we will take 
m = 6 and ignore the values of the differential coefficients at the 
upper limit. 

Then 

/(•)-(l-l)-, / , (^) = -(ri) a! l0g e l-l, /-(#) = -(ll)Ml0gel-l) S . 

Also , log. 1-1= -09531. 

Thus we get 



— o - 1 (M)- 10 =2-6988 

15 2 — 1 

X(l'l)- 10 x -09531= -6859 



X 


(1-1)- 


10 


•38554 


25 


•09230 


40 


•02209 


55 


•00529 


70 


•00127 


85 


•00030 


100 


•00007 



12 

15*- 1 
720 



x(l'l)- 1( >x(-09531) 3 = -0234 



•50686x15 = 7*6029 

And the final result is 

7-6029 - 2-6988 - '6859 + -0234 = 4-2416. 
The true value is, of course, 



( 1-1 ) '° = 4-2410, 
11 



I- 1 



which shows that the error involved by disregarding the final terms 
of the series is small. 
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Example 2. By Lubbock's formula. 

In this example we will take n = 10, using formula (6). 



X 


(1-1)-* 


A 


A 2 


A 3 


A 4 


10 


•38554 


-•23690 


+ •14557 


- -08946 


+ 05500 


20 


•14864 


•09133 


•05611 


•03446 




30 


•05731 


•03522 


•02165 


-|/(10) = -1-7349 


40 


•02209 


•01357 




50 


•00852 






qq 


60 


•00328 






120 A= - ' 1954 


70 


•00127 






99 


80 


•00049 






-^ A2=a - *° 600 


90 


•00019 






240 


100 


•00007 


10 = 6*2740 




99x1899 3 


— - -0233 


•62740 x 


720,000 












-S^ 










480,000 




-2-0238 



The result is 6*2740 - 2*0238 =* 4*2502 as compared with the true 
value of 4*2410. 

8. The above examples illustrate points of disadvantage which 
may arise in connection with the formulas. 

Woolhouse's formula can only be applied to a mathematical 
function where the differential coefficients can be obtained. In 
other cases, Lubbock's formula must be used, and then it may often 
happen that the terms in the formula do not converge sufficiently 
rapidly to give a good result, unless a large number of terms is 
used. This accounts for the relatively poor result shown in the 
example above. If intervals of 15 had been used instead of 10, a 
worse result would have been shown. 

Some of these disadvantages can be met by the use of the 
formulas given in the following Articles. 

Other formulas of Approximate Integration 

9. The definite integral of any function can be expressed in terms 
of the individual values of any number of ordinates by assuming 
that the function can be represented, to a sufficient degree of 
approximation, by a parabolic function of the requisite degree in x. 

f 2 

10. Thus to express the value of f(x) dx in terms of/(0),/(l) 

Jo 

and /(2), let f(x) = a + bx + cx\ 
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[N.B. Since we are to express the result in terms of three values 
of the function, we introduce three terms involving three unknown 
constants a, b and c in the expression representing/^).] 



r 2 r 2 

Then f(x)dx=\ (a + bx + ca?) dx 

Jo .' o 

f bx 2 ex*!? 



Also /(0) = a, 

/(l) = a + 6-fc, 
/(2) = a + 26 + 4c. 

These equations could be solved for a, 6 and c, and the resulting 
values substituted in the value of the integral found above. 
Alternatively, we may proceed as follows : 

Let (7<*)<fe=P/(0) + ff/(l) + »/(2). 

Jo 

8c 
Then 2a + 26 + -k =pa + ? (a + b + c) + r (a + 26 + 4c), 

so that jo + <? + r = 2, 

g + 2r = 2, 

<7 + 4r = f . 
Whence r-J, g = f, p = £, 

and r/ W &-/<°> + ^ 1 > + > f < 2 > (7). 

.' o o 

By analogy we may write 

n 



/: 



| o /(*)cte-g{/(0) + 4/(»)+/(2ii)} (8). 

This result is known as Simpson's Rule. 

By dividing the whole area from to 2n into n consecutive 
equal spaces and applying Simpson's rule to each, we obtain 

J*/(«) d*nm\ {/(0) + 4/(1) + 2/(2) + 4/(3) + 2/(4) + . . . 

+/(2n)}...(9). 
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This formula should normally yield better results than a single 
application of the formula (8) over the whole range of integration. 
It can, of course, only be applied where a sufficient number of 
values of the function to be integrated is available. 

11. Where the number of terms is four, by proceeding in a 
similar manner to the above, we arrive at the result 

j o f(x)dx = f {/(0) + 3/(l) + 3/(2)+/(3)}...(10). 

12. The process of obtaining the desired formula may be sim- 
plified somewhat, when the given terms are arranged symmetrically 
about the central point, by adopting the central point as origin. 

This may be illustrated by proving the well-known Weddle's 
Rule. 

In this case / f(x) dx is to be expressed in terms of 

/(0),/(l),.../(6). 

Let f(x) = a + bx + ex* + dx 3 + ex* +fx 5 + gx 6 . 

mu f" 4 " 3 ^/ si n -10 486 4374 

Then I f(x)dx= 6a + 18c +-p-e + -y- g. 

Also /(0) = a, 

/(l)+/(-l) = 2a + 2c + 2e + 2< 7 , 
/(2) +/(- 2) = 2a + 8c + 32e + 128g, 
/(3) +/(- 3) = 2a + 18c + 162e + U5Sg. 

Solving these equations for a, c, e and g, and substituting the 
resulting values of these constants in the expression for the in- 
tegral, we obtain 

jv- * 

J -3 



/(*) dx = ^ (272/(0) + 27/(1) +/(- 1) 



+ 216/(2) +/(- 2) + 41/(3) +/(- 3)}. 

This is not in a very convenient form for numerical work, so 
we add 



i5o A, -f<- 3) = i4o i- 2 °f (0) + 15 ^ (1) + f(~ x) 

- 6/(2)+/(-2) +/(3)+/(-3)}, 
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thus giving 

fj(x) dx + jL A «/(- 3) = ^ (252/(0) + 42/(1) +/(-!) 



+ 210/(2) + /(- 2) + 42/(3) +/(- 3)}. 
Neglecting the term r— ^ A«/(— 3) which will usually be very 
small, we arrive at the final result 
f + *f(x) dx = A {/( 3 ) + 5/(2) +/(1) + 6/(0) +/(- 1) 

" . . +5/(-2)+/(-3)}...(ll), 

or, changing the origin, 

/'/(*) dx = A {/(0) + 5/(1) +/(2) + 6/(3) +/(4) 

+ 5/(5)+/(6)}...(lt). 
13. Another powerful formula can be obtained by expressing 

jj(x)dx in terms of /(0), /(l), /(3), /(5) and /(6). 

Having five values of the function, let 

f{x) = a + bx + cx 2 + dx 9 + ex*. 

f(x)dx=Ga + 18c + ?Y-e. 

Also /(0) = a, 

/(2) + / (- 2) = 2a + 8c + 32e, 
/(3) +/(- 3) = 2a + 18c + 162e. 

Solving for a, c and e, and substituting in the equation for the 
integral, we find 

jj(x) dx = 2-2/(0) + 1-62 {/(2) +/(- 2)} + -28 {/(3) +/(- 3)}, 

or, altering the origin and limits, 

rQn 

J f(x) dx = n {"28/(0) +/(6») + 1*62/ (n) +/(5w) + 2*2/(3n)} 

* . (13). 

Similarly 

run 

/(*)(&= n {-28/ (6n)+/(12n) + r62/(7n) + /(llw) 
J en 

+ 2-2/(9n)}. 
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Then \f(x)dx=\ f(x)dx + f(x)dx+... 

Jo JO J Qn 

- n [-28 {/(0) + 2/(6n) + 2/(12n) + . . .} 

+ 1-62 {/(n) +/<*») +/(7n) + ...} 
+ 2-2 {/(3»)+/(9») + ...}]. 

Now in a series of values decreasing to zero, if In be so chosen 
as to fall just within or just without the limits of the table of the 
function to be integrated, we obtain the convenient formula 

{" f(x)dx = n {'28/(0) + l-62/(n) + 2*2/(3w) + l-62/(5n) 

+ -56/(8n)+ V62/(7n)}.'..(U). 

The formulas in this Article are due to G. F. Hardy and the 
last is usually known among actuaries as 39 (a), since that is the 
number assigned to it in the original Text Book, Part II. 

14. A useful formula, involving only a simple summation of 
certain terms of the given series, is as follows : 



/, 



'°7(«) <** = ^ [/(«) +/(*») +/(6») +/(9»)]- • -(15). 

* 



By expanding /(a?) in a series of ascending powers of x by 
Maclaurin's Theorem, it is easy to show that the formula involves 
a small second difference error. It can, however, be used con- 
veniently where only a rough result is required. 

15. An alternative method of obtaining formulas of this char- 
acter is to express f(x) in terms of the given values of the function. 
This can be done by Lagrange's formula and the resulting ex- 
pression is then integrated between the desired limits. 

As an example, we will develop formula (13) in this way. 

The given values of the function are/(3),/(2),/(0),/(- 2) and 
/(-3). Then 

( g +2)*(s-2)(*-3) „ (x + 3)x(x-2)(x-3) 

/(*) = -lx-3x-5x-6 /( ~ 3) + lx-2x-4x-5 ' ( "" 2) 

(a? +3) (a? + 2) (a? -2) (a? -3) , , m (a? + 3) (a? + 2) a? (a? - 3) ,, , 
+ 3x2x-2x-3 /W + 5x4x2x-l y v; 

(a? + 3)(a;+2)a;(a;-2) 
" + " 6x5x3x1 J w ' 
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Therefore 
f + /(.)d* = /(- 3) f + ^ + 2)*(*-2)0E-3) ^ + s . m . lar temg 

= -28/(-3) + l-62/(- 2) + 2-2/(0) + 1*62/(2) + -28/(3) 
as before. 

16. The above methods are perfectly general and afford means 
of expressing the value of an integral in terms of any given values 
of the function. 

It will be of value to indicate the method of application and to 
show the degree of accuracy in each case by applying them to an 
example. 

f 8 dr 

Thus — = 2 loe. 2 = 1-38630. 



0.21„ g .2-l-: 



(i) Simpson's rule applied once gives 

*(* + * + *) = 1-42500. 

(ii) Simpson's rule applied three times over the values 2-4, 4-6, 
6-8 gives 

Hi + * + ! + f + * + * + *)- 1-88770. 

(iii) The "three-eighths" rule gives 

§ (t + t + 1 + 4) -1-40635. 

(iv) Weddle's rule gives 

&(t + t + i + « + * + * + i)-l-88679. 
(v) Formula (13) gives 

{•28 (1 + *) + 1-62 (i + 1) + ~-l = T38643. 
(vi) If formula (15) is used, n = "6 and we get 

i[^ + H + ^ + ^]- lfl8889 ' 

As the values of the function are changing rapidly over the 
period used, the above is a somewhat severe test of the formulas. 
It will be noticed, however, that Weddle's rule and formula (13) 
differ from the true value only in the fourth place of decimals. 



CHAPTER XXI 

PROBABILITY 

1. In considering the subject of Probability a clear appreciation 
must be obtained of the distinction between its mathematical or 
theoretical treatment and its arithmetical or practical development 
which is the basis of actuarial science. The difference between 
these two points of view will become clear by examining a few 
simple illustrations. 

Probability or chance is merely an expression of relative degrees 
of uncertainty in relation to an event in respect of which our 
knowledge is incomplete. In the ordinary phenomena of life prob- 
abilities may not have a numerical value although distinctions 
are drawn. Thus we should say, in the absence of more precise 
knowledge, that a man of 30 was as likely as not to die before 
another man of 30, would most probably outlive a man of 60 and 
would almost certainly outlive a man of 80. The mind has clearly 
formed definite opinions on these points but the respective prob- 
abilities cannot be expressed in numerical form without further 
analysis. 

2. To obtain a measure of probability some unit must be taken, 
and for reasons of convenience it has been the universal custom to 
take absolute certainty as having unit probability. Starting from 
this standard, certain probabilities can be found from general 
reasoning, whereas others, such as the chance of dying in a year, 
or of becoming the father of twins, must be deduced from the ob- 
servation of suitable statistics. 

3. As an example of the former class, one may say that, on 
tossing a coin, it is certain that either "heads" or "tails" will appear, 
and as these possibilities are equally likely, the chance that heads 
will appear is clearly ^. 

We cannot however assume from this that, if we toss a coin four 
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times, two heads and two tails will necessarily appear. Indeed, as 
will be shown later, the respective probabilities are 

heads 4 tails T ^ 

1 head 3 tails J 

2 heads 2 tails § 

3 heads 1 tail J 

4 heads tails T ^ 

T_ 

What we understand is that, if the number of trials were 
sufficiently extended, the proportion of cases in which heads 
appeared would approximate more and more closely to £. It is in 
this sense that the mathematical definition of probability is to be 
understood in connection with actuarial work. 

The following Articles deal with the treatment of the subject 
from its mathematical aspect. 

4. From the foregoing considerations it will be evident that if an 
event can happen in 'a* ways and fail in l V ways, and each of these 

ways is equally likely, the probability of its happening is y , and 

that of its failing is , . For the sum of these two probabilities 

gives the chance of its either happening or failing, which is neces- 
sarily unity. 

Alternatively, we may say that the probability of an event is equal 
to the ratio of the number of cases favourable to the event, to the total 
number of cases. 

It follows that if p is the probability that an event will happen, 
the probability of its not happening is 1 — p. 

5. Another method of statement is to say that the odds are < a i 
to*V in favour of the event or 'b' to 'a' against the event 

Example. The chance of throwing a four at a single throw of a 
die is £, for there is one favourable result (a 4) and five unfavourable 
results (1, 2, 3, 5 or 6), all of w T hich are equally likely. 

The chance of not throwing a four is # , or, in other words, the 
odds are 5 to 1 against the event. 
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6. The following proposition is also self-evident : 

If an event can happen in more than one way (all ways being, 
however, mutually exclusive), the probability of its happening at all 
is the sum of the several probabilities of its happening in the several 
ways. 

Thus if the chance of scoring a "bull" be -fa, that of scoring an 
"inner" be £, that of scoring a "magpie" be £, and that of 
scoring an "outer" be ^, the total chance of hitting the target at 
all (all events being mutually exclusive) must be -^ + £ + £ + \ = £#. 

7. The solution of elementary questions in probability depends 
therefore upon general reasoning, but calculation is aided in some 
cases by the theorems of permutations and combinations. 

Some elementary examples will now be given : 

(1) The odds in a given race against three horses are 11 to 4, 
13 to 3, and 7 to 2 respectively. Find the chance that one of them 
will win the race, a dead-heat being assumed to be impossible. 

The chance that the first horse should win is ^, 
„ „ second „ „ -$$, 

third „ „ f. 

Thus the chance that one of them should win is 

T5 + tS + 9 = 7S17* 

(2) A has three shares in a lottery where there are three prizes 
and six blanks. B has one share in another, where there is one 
prize and two blanks. Show that A has a better chance of winning 
a prize than B in the ratio of 16 to 7. 

To be successful A may draw either 3, 2 or 1 prizes. 

He may draw 3 prizes in 1 way. 

He may draw 2 prizes and 1 blank in( 2 JxLj = 18 ways. 

He may draw 1 prize and 2 blanks in L) x (A = 45 ways. 

The total number of ways in which he can win at least one prize 
is therefore 1 + 18 + 45 = 64. 

Now three tickets can be selected in L) = 84 ways. 

Therefore A's chance of success is §£ = •£?-. 

B's chance is clearly J. 

Therefore A's chance is to B's chance in the ratio 16 : 7. 

Alternatively A may draw all blanks in (A = 20 ways. His 



128 PROBABILITY 

chance of non-success is therefore §£ and his chance of success 
1 — f£ = £f- as before. 

(3) If four cards be drawn from a pack, what is the chance that 
there will be one from each suit ? 

Four cards can be selected from the pack in ( J = 270725 ways. 

Four cards can be selected so as to be one from each suit in 
13 4 = 28561 ways. 

Therefore the required chance is ^n^V^T&Ws* 

(4) Out of a bag containing 12 balls, 5 are drawn and replaced, 
and afterwards 6 are drawn. Find the chance that exactly 3 balls 
were common to the two drawings. 

The total number of ways of making the second drawing is 

To comply with the conditions, it must contain 3 balls out of 
the first 5 chosen, and 3 balls out of the 7 left on the first choice. 

The respective ways of making these selections are UJ and (^j . 
The total number of selections favourable to the event 'is therefore 
(3) x (3) = ^®> an( * ^ e re( l u i re d chance is fff = |f . 

(5) Twelve persons take their places at a round table. What is 
the chance of two particular persons sitting together ? 

Let the two persons be A and B. Then, since we are only con- 
cerned with the relative positions of the persons, we may regard 
A's place as fixed. There are then 11 other seats in all, 2 of which 
are adjacent to A. B's chance of occupying one of these is there- 
fore &. 

(6) What is the chance of throwing more than 10 in a single 
throw with two dice ? 

A score of more than 10 can be made by the following throws : 

5 and 6, 6 and 5, 6 and 6. 
The total possible number of combinations is 6 x 6 = 36 and, 
as 3 of these are favourable to the event, the required chance is 



*■ 



It should be noted that exactly 11 can be thrown in two ways, 
since (regarding the throws as consecutive) either a 6 or a 5 may 
appear first, but 12 can only be thrown in one way. 



PROBABILITY 129 

8. If the chance of an event is p and the measure of the 
quantity dependent on the event be x, the product px is called the 
expectation. 

Thus we may have the expectation of a person in a lottery, or the 
expected value of a prize, or the expected number of deaths among 
a given number of persons, and so on. 

Examples, 

(1) What is the expectation of a person who is to draw one 
envelope from a bag which contains one £1 note, two 105. notes and 
three blank pieces of paper, each placed in an envelope of uniform 
size? 

The chance of drawing the £1 note is £. 
The chance of drawing a 10s. note is £. 
The value of the expectation is therefore 

t x 205. + £ x 105. = 65. Sd. 

(2) If the chance of dying in a year is ^, the expected number 
of deaths among 100 people is ■£$ x 100 = 5. 

9. In the foregoing articles we have considered what are, in 
principle, single events. We have now to determine the appropriate 
formulas for combinations of two or more events. 

In considering problems of this kind, a close watch must be made 
to see if the events are dependent or independent. 

Thus, a bag contains 6 white and 4 black balls, and it is desired 
to estimate the combined chance of drawing at the first draw 
3 white balls and at the second draw 3 black balls. If the balls are 
replaced after the first draw, the second event is clearly independent 
of the first. But, if they are not replaced, the drawing of 3 white 
balls at the first draw will obviously affect the chance of drawing 
3 black balls at the second draw, and the two events will be 
dependent. 

Dependent events are also called contingent. 

10. We are thus led to the following proposition: 

The chance of two independent events happening is the product of 
the chances of their happening severally. 

For if the first event can happen in a ways and fail in b ways, 
and the corresponding figures for the second event are a! and b\ the 
total possible combinations of events are (a + b) (a + &'). Of these 

H. t. b. i. 9 
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Both events may happen in aa! ways. 
The first event may happen and the second event 

fail in ab' ways. 
The first event may fail and the second event 

happen in ba ways. 

Both events may fail in bb' ways. 

The chance of their both happening is therefore 7 TN , . — rrr 

- 5 (a +6) (a' + 6') 

or — -r . , , , i.e. the product of the respective chances of their 

happening. 

If the respective chances are p and q, 

the chance that both happen is pq, 

the chance that the first happens and the 

second fails is p (1 — q), 

the chance that the first fails and the 

second happens is (1 — p) q, 

the chance that bpth fail is (1 — p) (1 — q), 

the chance that at least one happens is 1 — (1 —p) (1 — q) 

= p + q-pq, 
the chance that one and only one happens 

is 1>(1 -?) + (! -1>)S 

= p + q-2pq. 

11. The above results can be applied by a slight modification 
of reasoning to events which are not independent. 

Thus if p be the chance of the event happening, and q be the 
chance of a second event happening when the first has happened, 
then the chance that both should happen is pq. 

12. By successive applications of the above reasoning we can 
arrive at a formula for any number of events. 

Thus if p be the probability of an event A ; and, when A has 
happened, q be the probability of another event B; and when A and 
B have happened, r be the probability of another event 0; and so on 
for any number of events; the chance that all the events will happen 
is the product pqr 

It follows that if p be the chance that an event will happen in 
one trial, the chance that it will happen in each of a series of r 
trials is p r , and the chance that it will happen at least once in a 
series of r trials is 1 — (1 — p) r . 
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13. This result can be developed more generally in the following 
way: 

Ifpbe the probability of an event happening in one trial, what is 
the probability of its happening once, twice, three times, . . . exactly 
in n trials? 

Let q be the probability that the event does not happen, so that 
p + q = 1. Then if the event is to happen exactly r times, it must 
happen in each of a given combination of r of the trials and fail in 
each of the remaining (n — r) trials. The chance of this occurring 
is, as seen above, p r q n ~ r . But the particular set of r trials can be 

chosen in y\ ways, each of which is equally likely. Therefore the 

total chance of its happening exactly r times is ( n ) p r q n ~ r , or the 

term containing p r in the expansion of (p + q) n . 

Thus the successive terms of this expansion represent the 
probabilities of the event occurring respectively n, n — 1, n — 2, . . . 
times in n trials. 

14. The foregoing propositions can best be grasped by applying 
them to certain examples. 

(i) If four cards be drawn from a pack, what is the chance that 
there will be one from each suit? 

Let one card be drawn, which may be of any suit. The chance 
that a second card is of a different suit is ff , for there are 51 cards 
remaining, 39 of which will be of different suits from the first. 
Similarly the chance of drawing a third card of a different suit from 
the first two is §§, and that of a fourth different card is ^f. 

The combined chance is therefore ff x ft x if = sVWs- 

This result should be compared with that obtained in Example 3 
of § 7, where the problem was treated as that of a simple prob- 
ability. It is clear that the same result must be obtained whether 
the cards are treated as being drawn simultaneously or successively. 

(ii) A man throws a six-faced die until he gets an ace; he is to 
receive £1 if he succeeds at the first throw, ££ if he succeeds at 
the second throw, ££ if he succeeds at the third throw, and so on; 
given that log € 6 = 1*79176, find the value of his expectation. 

The chance of succeeding at the first throw is £. The chance of 
succeeding at the second throw is compounded of the chance of 
failing at the first throw and succeeding at the second throw, etc. 

9—2 
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•} 



The value of the expectation is therefore 
* + #■*•* + (#.*•*+•.• 

-i{- log. (i-Dl 

- -3583 = 75. 2d. 
(iii) Let it be assumed that the probabilities of dying within 
ten years after the ages specified are, on the average, as follows: 



Age 



30 
40 
50 



Probability of dying 
within next 10 years 






What is the chance 

(a) that a person A now aged 30 should die between the 
ages of 50 and 60 ; 

(6) that two persons A and B, aged respectively 30 and 
40, should be alive 10 years hence; 

(c) that of two persons A and B, aged respectively 30 and 
40, A should die between the ages of 40 and 50 and B should 
survive to the age of 60 ? 

(a) The required chance is compounded of the chances that 
A should survive successively to ages 40 and 50 and should then 
die within the next ten years. 

The several chances are (1 — ■&), (1 — $) and £; and the required 
chance is the product of these three factors, namely 

(b) The respective chances of surviving ten years are (1 - ^) 
and (1 — £). The chance that both A and B should survive that 
period is therefore 

tt*ir 

(c) The chance that A should die between the ages of 40 and 50 
is tt x £. 

The chance that B should survive to the age of 60 is § x f . 
The required chance that these two events should both happen 
is, therefore, 

<ttX*)x(tX*)-T«fr 



ttx|-W. 
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(jv) The faces of a die are marked with the consecutive numbers 
1, 2, ... 6. What is the chance that, after seven throws, the sum of 
the numbers exhibited equals 30 exactly ? 

This example is important, as it illustrates a method which can 
frequently be employed. 

The number of ways in which the seven numbers exhibited can 
total 30 is given by the coefficient of x M in the expansion of 

for this coefficient arises from the combination of the indices of x, 
taken together in such a way as to produce a total of 30. 

Writing (x + a? + . . • + afj as x 1 (1 - a«y (1 - x)-* it is seen that 
we require the coefficient of x 23 in the expansion of 

(1 -a«y (I - x)~ 7 = (1 -7a« + 2W 2 -35x ls + ...) (I -x)-**. 

The expression within the first bracket need not be expanded 
further, since no term higher than x 23 is required. It remains to 
combine the given terms with appropriate terms taken from the 
expansion of the second bracket in" such a way that the power of 
x given by the product of the two terms is 23. 

Thus we get 

(i) ^;;;; x x ^ )- 475,020^ 

(- V) x ( ^x.'.-xe 3 x ") = " 706 ' 629 *" 

(210 x ("^j;;; *" «*») = 259,896 *> 

(-SSx^xf 6 ,*-* 1 } a?) =- 16,170a; 23 
v / \lx...x6 / 

12,117 of 3 

The number of ways in which a total of exactly 30 can be made 
is thus 12,117. 

Now, on any of the seven throws, any of the six numbers may 
be exhibited. The total number of possible combinations of num- 
bers is, therefore, 6 7 . The number of combinations giving a total 

12 117 
of 30 being 12,117, the required chance is — — — . 
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15. The following miscellaneous examples are taken from the 
examination papers of the Institute ; they illustrate various devices 
which may be employed with advantage in the solution of questions 
of this character. 

(i) If n whole numbers be multiplied together, find the chance 
that the last digit of the product is a five. [1910, Paper I, Q. 7.] 

The chance is compounded of the separate chances that none of 
the final digits of the n numbers is even and that one at least is a five. 

The required chance is therefore 

Kn _A.1I 

(i) n {i-(*)"} = ^-- 

(ii) Four coins are tossed together and A is to receive £2 if 
exactly 2 heads turn up, and to pay £1 in any other event. Find 

the probability that after four trials A is £1 out of pocket. 

[1913, Paper I, Q. 8.] 

To satisfy the conditions A must win once and lose thrice. 

The chance that exactly 2 heads turn up equals the middle term 

in the expansion of (£ + £) 4 = f *) * (i) 4 = §. 

Therefore the chance that he should win once and lose thrice 
equals the second term in the expansion of (§ 4- § ) 4 , which is 

(iii) If a coin be tossed 15 times, what is the probability of 
getting heads exactly as many times in the first 10 throws as in 
the last 5? [1915, Paper I, Q. 9.] 

If n coins be tossed the chance that exactly r heads turn up 
equals the (r + l)th term in the expansion of (^ -f \) n . 

Now in the last 5 throws we may get 0, 1, ... 5 heads and we have 
to combine the chance of any of these with the chance of getting 
the same number of heads in the first 10 throws. Thus we have 
for the several chances of getting 

heads in both sets of throws = (^-) 5 
lhead „ „ „ = 5(J) 5 



2 heads 

3 „ 

4 „ 

5 „ 



= 10(i) 5 
= 10(i)» 



(4) 10 

10(i) M 
.120(i)» 



= 5<iy.210(i)» 
= (i)».252(i)» 
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Therefore the total chance that the same number will turn up in 
both sets of throws is the sum of each of these distinct probabilities, 
3003 

(iv) A and B throw for a certain stake, each throwing with one 
die; A's die is marked % 3, 4, 5, 6, 7, and &s 1, 2, 3, 4, 5, 6, and 
equal throws divide the stake; prove that A's expectation is ffof 
the stake. 

What will A's expectation be if equal throws go for nothing? 
[1911, Paper I, Q. 8.] 

Let e equal A's expectation. 

The chance of A's throwing any of the numbers marked on the 
die is the same for each number. 

If A throws a 2, B must throw a 1 if A is to win the stake, or 
a 2 if A is to divide the stake. The chance that A throws a 2 is £; 
the chance that he then wins the stake is £ and that he divides the 
stake is also £. His expectation if he throws a 2 is therefore £ -f £ . J. 

If equal throws go for nothing, A's expectation after an equal 
throw clearly remains at e. In that case his expectation if he throws 
a 2 is £ + £e. 

We therefore have the following scheme: 



A'a throw 


Chance 


A's expectation if equal 


A's expectation if equal 


thereof 


throws divide the stake 


throws go for nothing 


2 


* 


*+*•* 


*+*• 


3 


* 


*+*.* 


*+*• 


4 


i 


*+*■* 


*+*« 


5 


i 


m-i 


m« 


6 


i 


#+*■* 


Hie 


7 


i 


i 


1 



The total expectation being the sum of the separate expectations, 
we have in the first case e = fy, being the sum of the figures in the 
third column of the above table, and in the second case e = ^ + -fee. 

Whence e = ||. 

(v) A man tosses 20 pennies and removes all that fall head up; 
he then tosses the remainder and then removes all that fall head 
up, and so on. How many times ought he to be allowed to repeat 
this operation if he is to have an even chance of removing all the 
pennies before he has finished? [1906, Paper I, Q. 11.] 
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The problem is clearly the same if all the pennies are tossed 
each time; we then have to find the chance that all the pennies 
have turned up heads once at least 

If n be the required number of throws, the chance that any- 
particular penny has turned up heads once at least is {1 —(J) 71 }. 

The chance that all the pennies have turned up heads once at 
least is therefore {1 — (J) 71 } 20 and by the terms of the question this 
must equal \. 

Solving this equation we find (£) n = 1 — (J) 20 , 

n = 4*87. 

— > 

(vi) 2 n players of equal skill enter for a tournament; they are 
drawn in pairs, and the winners of each round are drawn again for 
the next. Find the probability that two given competitors will play 
against each other in the course of the tournament. If n = 5, show 
that the probability that a given player will either win or be 
beaten by the actual winner is ^. [1912, Paper I, Q. 6.] 

As regards the first part of the question, the total number of 
games that will be played is 

2 n-i + 2 n - 2 +...+l=2 n -l. 

2 n (2 n — 1) 
Also two players can be selected in ^ ways. The chance 

that two will meet is therefore the quotient of these two values, 
1 

The second part of the question can be proved by an inductive 
process. 

Let u n be the required chance. 

Then in the first round he may either win or lose. If he wins, 
he passes into the next round and his chance of winning or being 
beaten by the ultimate winner then becomes u n ^ 1 . If he loses, his 
opponent passes into the next round where the latter's chance of 

becoming the ultimate winner is increased to -^—^ . 
We thus have 

I 1 _1_ 



PROBABILITY 4 137 



Q- 'I 1 1 1 11 

Similarly g a n _, = ^ ^ 2 + j . — 2 , 

11 11 



etc. etc. \ 



2»-2 ^ 2 ~~ 2 n_1 ^ 2 n_1 ' 2 * 
As ^! obviously equals 1, we have by addition 

1 ?i-l n+l 



' n 2 n_1 + ~~2**~ ' 2 n ' 

It is instructive also to prove the first part of the question 
inductively. 

Let u n be the required chance. This chance is compounded of 
the two chances that they meet in the first round, or that, not 
having met in the first round, they both pass into the second 
round where the corresponding chance becomes iV-i» 

The chance that they meet in the first round is ^ — -; the chance 

2 W — 2 
that they do not meet is ^ — =- . Therefore we have 

1 1 2" -2 

Un " 2nZi +u ^ ' 4 • 2» - 1 ' 

or (2» - 1) Un ~ 1 + n ~ u^. 

(2«— 1 — 1\ * 1 2 n ~ 2 — 1 
^ ) u ^-i = 2 + 4 — Un ~*> 

/2«-2-l\ 1 2"- 8 - 1 

V 4 J Un ~ 2 ~ 4 + 8 Wn " 3 ' 

etc. etc. 



/ 2 3 ~1 \ _ J_ 2 2 -l 
\ 2 W ~ 3 / Ud ~~ 2 n ~ 8 "*" 2 n ~ 3 Ws> 

/ 2 2 -l \ _ _1 

^ 2n-a y^ 2 "2 n " 



2-1 
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By addition 




(2»- 


-l)tfe-l + 2+3+-H 


(since u x obviously equals 1) 




1-1 

2 n 2 n -l 


- 


2 


Whence 


1 

^» — 2^=1 » 


as before. 





16. The application of the Integral Calculus to problems of 
mean value and probability is shown in Chapter XIX, §§ 5-7. 
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CHAPTER III 

1. Given 

/(l) = (*-2)(*-3), /(2) = (*-7)(#-5), /(B)-(*-l6)(*+l), /(4) = 63, 
obtain a value for #, assuming second differences are constant. 

2. Find the nth term of the series 2, 12, 36, 98, 270, 768, etc. and the sum 
of n terms of it. 

3. Given that 

/(0) = 66566, /(1)=65152, /(2) = 63677, /(3) = 62136, 
find/(9). 

4. Find A- 6 * +29 



(3* + l)(3ff+4)(3a?+7)(aaM-10)" 

5. Find the value of 

A»(3ar + 1) (3#+4) (3^7 + 7) ... (3*+3n-2). 

6. Prove that 

(7i-l)^/(^)+77i(^-l)^-V(^ + l) + ^y-\^-l) m " 2 /(^+ 2 ) 
+ ...+/(#+m)=w w »/(;r)-i-mn Tn - 1 A/ (a;) 

+^ff-^ ) « m - 2 Ay(^)t...+A«/W 
Hence find the sum of the series 

l^+m.2H m( ^" 1) 32+... + (^ + l)l 

CHAPTER IV 

1. Given 

/(0)=70795, /(1) = 72444, /(2) = 74131, and /(6) = 81283, 
find /(3),/(4), and /(5). 

2. Find /(35) given 

/(20)='01313, /(30) = *01727, /(40) = -02392, /(50)="03493. 

3. Supply the missing term in the following table : 

/(0) = 72795, 
/(1) = 71651, 
/ (2) = 70458, 

/(4)=67919, 
/(5) = 66566, 
/(6)=65152. 

4. Given 

/(0) = 11, /(3) = 18, /(6) = 74, and/(9) = 522, 

find the intermediate terms. 
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. 5. Given 

/(0)=98203, /(1) = 97843, /(2)=97459, /(3)=97034, 
find/(2'25). 

6. Given 

/(0) = 98023, /(10) = 97651, /(20) = 97246, /(30) = 96802, 
find/(15). 

7. Given 

/(0) = 58*842, /(2) = 55'257, /(4)=51'368, /(10) = 37'977, 
complete the series /(0),/(l), .../(10). 

8. If you were asked at very short notice to obtain approximate values for 
the complete series /(0), /(l), /(2), .../(20), being given that /(0) = '013, 
/(10)='248, /(15) = '578, and /(20) = -983, what methods would you adopt, 
and what value would you obtain for/(9)? 

9. /(1)=1; /(2)+/<3)-6-41; /<4)+/(5)+/(6)- 18-47; 

/(7)+/(8)+/(9)+/(10)+/(ll)+/(12)=90'36. 
Find the value of f(x) for all values of x from 1 to 12 inclusive. 

10. Apply Lagrange's formula to find/ (5) and/(6), given that 

/(1) = 2, /(2)=4, /(3)=8, /(4)=16and/(7) = 128; 
and explain why the results differ from those obtained by completing the series 
of powers of 2. 

11. Find the simplest algebraic expression in x which has the values 5, 3, 
9, 47 and 165 when x has the values 0, 1, 2, 3 and 4 respectively. 

12. Prove the following formulas for approximate interpolation : 

/(l)=/(3)-3|y(5)-/(-3)]+-2|y(-3)-/(-5)] (1), 

/(0)=i[/(l)+/(-l)]-i[H/(3)-/(l)}-i{/(-l)-/(-3)}]...(2), 
and apply them to find the logs of 45, 46, 47, 48, 49, being given 
log 42 = 1-62325, log 50= 1 '69897, 
log 44=1-64345, log 52 = 1'71600. 

13. Given 

10 10 10 

2/00 = 500426, 2/ (x) = 329240, 2/(#) = 175212 and / (10) =40365, 
find/(l). 

CHAPTER V 

1. Use Gauss' interpolation formula to obtain the value of/ (41) given 
/(30)=3678-2, /(35) = 2995-l, /(40)=2400'1, 
/(45)=1876-2, /(50) = 1416-3. 
Verify your result by using Lagrange's formula over the same figures. 
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2. Given the following table find /(28) using Stirling's formula : 
/(20) = 98450, /(25) = 96632, /(30) = 94472, 
/(35) = 91852, /(40) = 88613. 



3. Prove that 

2! " 2 



/W-i[/(-*)+/(4)l+^(-*>+^* • *Zfcf>+*Kfc# + ... 



and apply the formula to find/ (32) given 

/(25)=-2707, /(30) = -3027, /(35) = -3386, /(40J--3794. 

4. From the table of annual net premiums given below find the annual net 
premium at age 25 by means of Bessel's formula : 

Age Annual Net Premiums 
20 -01427 

24 -01581 

28 -01772 

32 -01996 

5. Use Everett's interpolation formula to complete the series /(25) to /(35), 
given that 

/(15) = 305, /(20)=457, /(25)=568, /(30)=671, 

/(35) = 897, /(40) = 1190, /(45) = 1481. 



CHAPTER VI 

1. Given the following table of f(x) : 

/(0) = 217, /(1) = 140, /(2) = 23, /(3)=-6, , 
show how to find approximately the value of x for which the function is zero. 

2. Given that, when 

#=0, /(*)«<), 

*-l, /(#) = 100, 

x=% /(#) = 2000, 
find x when f(x) =1900 by Lagrange's formula of interpolation (applied 
inversely) and explain why the result does not agree with that found by using 
the formula/(#) = (l + A) x /(0) and solving the quadratic, 

3. The following values of /(#) are given : 

/(10) = 1754, /(15)=2648, /(20)=3564. 
Find, correct to one decimal place, the value of x for which f{x) — 3000. 

4. /(30)=-30, /(34)- -13, /(38) = 3, /(42) = 18. 

Apply Lagrange's formula of interpolation inversely to find x, where f(x) =0. 
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CHAPTER VII 



1. Find the value of 



1 



7l(?l + l)(7i + 2)...(?^ + 7•-l) , 

and use the result to find the value of 

* 1 

2- 



2 n(n + l)(n+2)' 

2. Prove that 

2(#2+l). xl=n.(n+l)\. 

3. Show that 

where <j> (x) is any rational integral algebraic function of x. . 

4. Prove that if a„ diminishes as n increases and converges to the limit zero, 
the sum to infinity of the series a\ — a 2 + a 3 ••• is ^ ne same as the sum to infinity 
of the series \a x — £Ao2+|A 2 a 3 - etc. 

Find the sum to infinity of the series ^-^x+ ^ ... true to four decimal 
places. 

5. Prove that if the fourth and higher differences are ignored 'the sum of n 
successive terms of a function of which /(0) is the central term is 



«/(o)+^Ay(-i), 



where n is an odd number. 



CHAPTER VIII 

1. Given the following values of/(#): 

/ (1-41) = 7092, /(1-49)=*6711, 

/(l -62) = -6579, /(l -53) = -6536, 
find/(l*45). 

2. Use Divided Differences to find/ (80) to the nearest integer, given 

/(70) = 235, /(71) = 256, /(79)=436, /(81) = 484. 

3. Given 

/(20) = -342, /(23) = -391, /(31) = -515, /(34) = -559, 
find/ (30) by means of Divided Differences, and check the result by applying 
Lagrange's interpolation formula. 
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CHAPTER IX 

1. The following values of /(#, y) are given : 

/(35, 55)= 10*020, /(35, 50) = 11*196, /(35, 45) = 12-019, 
/(40, 55)= 9*796, /(40, 50) = 10'894, /(40, 45) = 11 -641, 
/(45, 55)= 9*583, / (45, 50) = 10*591, /(45, 45) = 11'243. 

(i) Using only six of the above values, find/(42, 52). 

(ii) Making use of all the data calculate /(44, 51). 

2. Prove that if /(0, 1)=/(1, 0) and/(0, 2)=/(2, 0), then 

/(*, ^^/^OJ + ^+^^+^tf^A^/CO, 0) + sy[/(l, l)-/(2, 0)]. 
Find/ (39, 33), given 

/(35, 35) = 3'151, /(35, 45)=3'912, 

/(35, 40) =3-471, /(40, 40) = 3-766. 



CHAPTER XI 

i 
1. Find the differential coefficients with regard to x of log#* and 

V^3 



(#-1)^-7 

2. Differentiate with respect to x 

, Z / a*+a?x+ax*+a* 
g V at + ax+x* • 

3. Find the differential coefficient with regard to x 2 of 

V# 2 -2 



4. Differentiate 



©X ff 

log -, with respecVto log - , 

XCf" 

(2) log ——j , with respect to x, 

(3) x&°«*)\ with respect to x. 



5. If(l-^+(l-^)*=0, nnd^. 

6. Obtain the differential coefficients of 



(1) lo g^> ( 2 ) «**> witn respect to x. 
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CHAPTER XII 

1. Find (1) the nth differential coefficient of — — — -^ with respect to 

\#~r *■) [X — 1) 
x, and (2) the second differential coefficient of log a (1 +x) with respect to log« x, 

2. Find the nth differential coefficient of eP* {x - 2) 2 . 

3. Having given that 

X dx*+ X dx+ y -^ 
prove that ^^ +( s» + i)*. ^ +( ^ +1) g, a 

4. Obtain the second differential coefficient of log (a + bx + ex 2 ) with respect 
to six, 

5. Find the nth differential coefficient with respect to x of 

x+1 
2x 2 -5x+2* 

CHAPTER XIII 

1. Prove that ^-r=l- % + ^ - ^ ... 

e?-l 2 12 720 

and show that no odd power of x beyond the first can occur. 

o d it, r log(#+Vl+# 2 ) _ ,,2.4 . 

2. Prove that &v , * =x-&a?+— -x 5 - .... 

n/1+# 2 s 3.5 

i 

3. Prove that the first three terms of the expansion of (1 +x) x in powers of 

x are e - — 4- i} ex 2 , 

4. Expand log (x + "/a 2 + x 2 ) in ascending powers of x. 

5. Expand — ^ — - — ' in ascending powers of a? as far as the term in- 
volving x 6 , 

CHAPTER XIV 

2. Find the maximum value of x (x - 1) (x - 2) between the limits and 1. 

3. Find the values of x at the points where the graph of the function 
(1 +x 2 ) e 1 "* 2 has its greatest slope. 
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4. A man in a boat at sea, 5 miles distant from the nearest point of a 
straight shore, wishes to reach a place 12 miles distant along the shore, 
measuring from this nearest point. At what point should he land to reach this 
place in the minimum time, if he can row at 3 miles an hour and walk at 
4 miles an hour? 

5. Given /(0) = 1876, /(1)= 777, /(3) = 19, and/(6) = - 218, interpolate the 
values of /(2), /(4), and /(5) and find the values of x for which f(x) is a 
maximum or minimum. 

6. Find the minimum and maximum values of 

#4_4#3_8#2 + 4&r-48. 

7. A window is in shape a rectangle with a semicircle covering the top. If 
the perimeter of the window be a fixed length £, find what is its maximum area. 



8. Find (1) Lt 7 , 

(2) LtCL+a,- 3 )* 2 . 



CHAPTER XV 

d z f(x) 

1. Show that —^-^=A 3 /(^-f) approximately. 

By considering the function f(x)=a+bx+d B and using the above relation 
prove that loggC=c*-c~* approximately, where c is a small quantity. 

2. Show that AO" 1 - ~0 m +^- m ...=0, when m > 1. 

3. Prove that, if f(x) be a function the fourth differences of which are 
constant, 

^P=Mf(*- 2)- 8/(* - 1)+8/(*+ 1)-/(*+2)] 

and hence find an approximate value for -r- [log f(x)] where 

f(x- 2) = 42-699, f(x-\) = 40*365, /(#) = 37*977, 
f(x+l) = 35-543, f(x+ 2) = 33*075. 

CHAPTER XVIII 

1. Find the integrals 

m/jra* »/A* mj 11 ?* 1 * 



\x*a x dx. 



2. Evaluate 

H. T. B. I. 10 
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3. Evaluate 



(D/^rn&d. 00/^. (3)/^- 



4. Find the value of J -^-r — ova/ . , N c?#. 



/•27^-l71^+256 
2(#-3) 3 (#+4) 
5. Evaluate 



/ (jy--l) a tg- 7) <to; J l °Z( l -**)<**'> jWog*) 3 ^- 



CHAPTER XIX 

1. If /(#)=* a + for +c# 2 , find expressions for /(J), ^ ■ and / f(x)dx in 

terms of/(0),/(l) and/(2). 

ri ^ 

2. Evaluate | — — r-rr*^- 

3. Find the value of ^2 ^9^-2 

4. Find the average value of x {x— 1) (a? — 2) between the limits and 1. 

5. Find the value of / x 2n - l e?dx. 

CHAPTER XX 

1. From the table 

/(0)=217, /(1)=140, /(8)-88, /(3)=-6, 

find an approximate value for / f(x) dx, and explain why the result differs 
from/(0)+/(l)+/(2)+/(3). 

2. Prove that approximately 

J*|/(*)^=&{/(-l) + 22/(0)+/(l)}, 
and find thereby an approximate value for 

f H 10,000 x 4~* *2 y dx. 

3. Discuss the error in assuming that 

J l 9 /(x) «fc-/(i)-i[/«»+/(i)] 

if (1) f(x)=a+bx+ca?, (2) /(#) = 100x4*. 
4 Given /( - 2), /(0), /(2), find an approximate value for 

f[ 3 f(x)-dx. 
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5. Prove that 



/: 



f(x)dx=& [ 5 /0Q + 8/(0) -/(-!)] approximately. 



If the speed of a train on a non-stop run is as shown in the table below, find 
the approximate mileage travelled between 12.0 and 12.30, using the above 
formula. 



Time 




Speed in 
miles per hour 


11.50 
12.0 
12.10 
12.20 
12.30 . 


242 
35 
41-3 
42-8 
39-2 



6. If/(1) = 4157, /(2) = 4527, /(4)=5435, find approximately the value of 

f{x) dx. 



i: 



CHAPTER XXI 

1. n persons are sitting at a round table, and from them three are selected 
at random ; show that the chance that no two of those selected are sitting next 

one another is 



(n-4)(n-5) 



(n-l)(n- 2)' 

2. A heap of playing cards contains 6 hearts, 5 spades and 4 clubs. A card is 
chosen at random 9 times in succession and is not replaced. Find the chance 
(1) that there are no hearts left, and (2) only hearts are left. 

3. A man has two sovereigns and four shillings in his pocket ; from these 
he selects at random 3 coins, and again chooses from these 3 coins at random, 
the first coin for A, the second for 2?, and the third for C. 

Find the values of the expectations of A, B, and 0. 

4. Find the probability that out of 5 persons aged 45 exactly 3 will die in 
a year. 

Probability that a person aged 45 will die in a year ='01224. 

5. A die with six faces is thrown three times and the sum of the throws is 
twelve : find the chance (1) that the first throw was a four : (2) that four was 
thrown each time. 

6. Two men, A and B 9 each draw a card from a well-shuffled pack of 
playing cards, find that they are of the same value, and replace the cards ; 
they do this four times in succession. Find the chance that this would happen 

1 



and show that it is approximately 



y 



84,000" 
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7. The probability that a man aged 50 will survive one year is -98428. 
Show that the probability that, out of 5 men aged 50, 3 at least will die within 
a year is -0000385. 

8. The 26 letters of the alphabet are placed in a bag. A and B alternately 
draw a letter from the bag, the letters drawn not being replaced. The winner 
is the one who draws most vowels. A starts and draws a vowel with his first 
draw. What is his chance of winning ? 

9. If a number of five figures containing any five of the ten digits once only 
is written down at random, what is the probability that it is divisible by 9? 

10. Given the following table find the probability that one at least of three 
persons aged respectively 20, 30, 40 will die between the 10th and 20th year 
from now : 



Age 


Probability 

of surviving 

10 years 


Probability 

of surviving 

20 years 


20 
30 
40 


•93363 
•91740 

•88476 


•85651 
•81167 
•71517 



11. Two persons, A and B, play for a stake, each throwing alternately two 
dice, A commencing. A wins if he throws six, B if he throws seven, the game 
ceasing as soon as either event happens. What ratio will A's chance of winning 
bear toil's? 

12. The sum of two positive integers (excluding zero) is 100 ; find the chance 
that their product exceeds 1200. 

13. The following table shows the probability that a woman of the age 
specified will marry in a year : 

Age Probability of marriage 

20 -0665 

25 -1033 

30 -0649 

40 -0183 
Find the probability that, out of 4 women aged 20, 25, 30, 40 respectively, only 
one marries within a year. 

14. A bag contains 8 counters, numbered 1 to 8. Four are drawn at 
random. Find the chances that 

(1) The sum of the numbers on the four counters amounts to at least 17. 

(2) The counters numbered 2 and 3 are among the four. 

(3) The four counters contain at least two of the three counters numbered 
3, 5 and 7. 

15. A penny is tossed six times. Find the chance that neither heads nor 
tails have occurred three times in succession. 
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16. An urn contains counters marked with the digits 6, 7, 8 and 9 ; and the 
number of times each digit occurs is equal to the value of the digit. If counters 
are drawn one at a time, each counter being replaced when drawn, what is the 
probability 

(1) that the digit 6 is drawn before the digit 9 ; 

(2) that the sum of the first three digits drawn is exactly 20? 

17. A and B play a set of games, to be won by the player who first wins four 
games, with the condition that if they each win three they are to play the best 
of three to decide the set. A'a chance*oi winning a single game is to B'a as 
2 to 1. Find their respective chances of winning the set. 

18. The probability of any one of 10 men each aged 30 surviving a year is 
•99229. 

Show that the probability that exactly 5 men out. of the 10 survive a year 
is 6-6 x 10 ~ 9 . 

Find also the probability that of the 10 men one particular man will die 
first and another particular man last. 

19. A point is taken at random within the area bounded by the curve 
y= x log x, the x axis, and the ordinates at the points #=1, and x=4. 

Find the probability that the distance of the point from the y axis is less 
than 2. 

20. In a game of whist the dealer found that on turning up the last card he i 
had the Ace, King, Queen, Knave, Ten, and 3 other trumps in his hand. Find 
the chance that this would occur. 

21. A and B cut a pack of cards, the player who wins the cut six times to 
be the winner. A, having won four times to 2?'s once, cuts a five. Find the 
chance that A will be the winner. 

22. In a line AB of length 3a, a point P is taken at random and then in 
AP a point Q is taken at random. What is the probability that P@.exceeds a ? 
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4 (1) aP»(^"^loga.log£ +a »logf (kg2+l) + l], 

00 i-5^. (3) ^T^ 2 - 

5 - - ;g:g; , :g:g:5g:-g - «• l-*« <"***+* 
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Maximum value when #= 
Minimum value when x 
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